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ABSTRACT 


This  report  contains  a summary  ot  intercal location 
measurements  carnea  out  on  four  different  reactors  unuer  various 
irraaiation  conditions.  Tnese  reactors  were  the  Northrop  TRIGA, 
SANDIA  SPR  II  ano  ACPR,  and  the  Aberdeen  APRF  reactor.  The  inter- 
calibration was  carried  out  by  (a)  measuring  neutron  spectra  for  a 
variety  of  oasic  and  modified  ("filtered")  reactor  neutron  fields 
with  a standard  threshold-foil-activation  neutron  spectrometry 
method,  (b)  folding  in  each  spectrum  0(E)  with  the  silicon  carnage 
function  D(E)  to  obtain  ^e(j#  the  number  of  near-1  MeV  neutrons 
required  to  proouce  the  same  damage,  (c)  irradiating  a batch  of 
2N2222A  transistors  simultaneously  with  the  neutron  oosimetry  stack 
and  measuring  the  total  damage  D * K 0 , and  (d)  dividing  D = K <t 
by  0 to  obtain  values  oi  the  carnage  coetticient  K for  the 
diverse  neutron  fields  measured. 

The  values  of  K were  constant  (±5%  for  1 standard  devia- 
tion) for  all  the  spectra  investigated,  indicating  that  ade- 

quately characterizes  the  permanent  radiation  carnage  effectiveness 
from  one  oasic  (and/or  modified)  reactor  spectrum  to  another  for  a 
practical  type  of  transistor. 

These  results  suggest  that  can  quantitatively 

relate  reactor-spectrum  and  threat-spectrum  radiation  damage,  since 
a variety  of  spectral  shapes  was  utilized.  For  the  FBR  leakage 
spectrum  modified  by  iron  and  polyethylene,  the  spectral  shapes 
above  the  0.2  MeV  "threshold"  for  radiation  damage  in  silicon  were 
signiticantly  different  from  all  the  others. 

This  work  is  a continuation  of  a Defense  Nuclear 
Agency  ( DNA)  sponsored  effort,  initiated  in  1974,  in  which 
several  different  spectra  (wSMR  FBR  glory  hole  and  leakage , 
and  TRIGA  leakage)  were  measured  ana  variational  studies 
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carried  out  to  determine  the  sensitivity  of  to  tne  many 

possible  sources  of  error.  To  make  this  report  as  complete  as 
possible  tor  the  user  these  previous  results  have  been  included 
as  Appendix  A. 

The  culmination  of  this  work  was  a set  of  ASTH 
uratt  standards  (Appendix  B)  which  were  followed  in  detail, 
except  tor  an  important  improvement  recently  developed  by  one 
of  the  authors  of  this  report  (VVV).  This  involves  a simple 
and  highly  accurate  method  of  measuring  the  error  in  ganuna-ray 
detector  efficiency  that  arises  tor  some  cases  where  cascade- 
gamma-ray  emissions  occur  (see  presentation  of  this  work  in 
the  text).  Once  measured,  the  error  is  corrected  tor  in  a 
very  simple  way. 

One  of  the  outputs  of  the  intercal ibr at  ion  effort 
is  a value  of  *eq/S  for  each  reactor  facility,  so  that  *Aeq 
can  be  determined  for  each  irradiation  by  the  sulfur-toil  acti- 
vation S (i.e.  the  i2P  disintegration  rate)  when  the  sulfur- 
foil  counter  is  properly  calibrated,  as  described  in  this  report. 


INTRODUCTION 


This  work  is  a continuation  of  an  effort  initiated 
in  1975,  when  a series  of  reactor  spectral  measurements  and 
calculations  were  carried  out  on  the  White  Sands  Missile 
Range  (WSMR)  fast  burst  reactor  (FBR)  and  the  General  Atomic 
(GA)  TRIGA  reactor.  The  neutron  spectral  measurements  were 
made  with  threshold  activation  foils,  and  a sensitivity  study 
was  carried  out  on  the  effect  of  foil-counting  errors,  foil 
activation  cross  section  errors,  choice  of  trial  spectrum,  and 
choice  of  reactor  unfolding  code.  These  results  appear  as 
Appendix  A of  this  report.  The  standard  neutron  threshold  foil 
spectrometry  method  appears  as  Appendix  B,  as  submitted  in  the 
form  of  a set  of  draft  ASTM  Standard  Methods. 

A neutron  field  can  be  uniquely  specified  by  a table 
of  spectral  intensity  .versus  neutron  energy.  An  accurate  spec- 
tral determination  is  prohibitively  expensive  and  in  general  can 
be  carried  out  only  once.  In  addition,  it  involves  the  carrying 
of  a large  number  set. 

Historically,  the  neutron  field  has  been  monitored 
from  irradiation  to  irradiation  by  means  of  sulfur  foil  activa- 
tion. These  foils,  having  an  activation  threshold  of  3 MeV, 
only  sense  the  upper  "tail"  of  the  spectrum.  The  spectral 
flux  can  be  accurately  inferred  only  if  the  spectral  shape  and 
the  sulfur  activation  cross  section  a'(E)  are  known  accurately. 
Further  quantification  had  been  applied  by  utilizing  the  Pu/S 
ratio,  since  Pu  has  an  effective  threshold  at  0.001-0.01  MeV 
when  adequately  shielded  inside  a ^B  shell.  The  S-foil  activa- 
tion and  the  Pu/S  ratio,  although  further  quantifying  the  neutron 
field,  fall  far  short  of  doing  so  to  the  extent  that  work  car- 
ried out  with,  say,  a TRIGA  reactor  can  be  directly  compared  to 


Note:  Manuscript  submitted  December  12,  1978. 
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FBR  irradiations. 


In  the  present  work,  neutron  spectra  (E)  were 
measured  for  a variety  TRIGA  and  FBR  facilities  both  for 
the  free  field  case  and  for  these  fields  modified  by  "filters". 
These  spectra  were  characterized  and  quantified  in  terms  of 
radiation  damage  effectiveness  for  silicon  devices  by  a single 
parameter  $eq-  Here,  is  the  equivalent  near-l-MeV 

fluence  for  producing  the  same  radiation  damage: 


eq 


-/ 

.01 


18  MeV 

$(E)  D(E)dE 


/ 


/ 


0(0.85-1.15  MeV) 


(1) 


Here.  D(E)  is  the  calculated  energy  dependent  neutron  cross 
section  for  producing  permanent  atomic  displacements  in  bulk 
silicon. 

Both  $ and  <J>  /S  are  presented  below  so  chat  4 

eq  eq  r eq 

can  be  determined  from  the  sulfur  activation  S,  following 
proper  (absolute)  calibration  of  the  sulfur-foil  beta  counter 
(see  below) . 

In  order  to  determine  the  validity  of  some 

measurements  of  radiation  damage  D - were  also  made  on 
batches  of  2N2222A  transistors.  These  were  simultaneously 
exposed  with  the  neutron  threshold  foil  spectrometry  stack, 
so  that  both  and  D ■ K4>  could  be  determined  for  the  same 
exposure.  From  the  ratios  D - K4>  to  <>  , the  damage  coeffic- 

ients K were  determined.  The  extent  to  which  the  measured  K 
are  constant  over  a wide  variety  of  neutron  spectra  determines 
the  extent  to  which  characterizes  the  radiation  damage 
from  any  of  the  radiation  fields  measured  in  this  task. 
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The  free-field  FBR  spectrum  was  modified  with  two 
different  filters  (10  cm  of  polyethylene  and  5 cm  of  iron)  in 
order  to  check  <t>e(j  over  a more  diverse  set  of  radiation  fields 
This  yielded  a total  of  5 different  types  of  spectra  (TRIGA, 
FBR  glory  hole,  FBR  leakage,  and  two  modified  leakage  spectra) 


2.  RADIATION  FACILITIES  INTERCALIBRATED 

A list  of  the  different  neutron  spectra  measured 
previously  and  in  this  program  is  presented  in  Table  1.  This 
list  covers  all  the  major  reactor  facilities  extensively  used 
for  radiation  damage  studies  in  the  United  States  at  this  time. 

Table  1.  Reactor  Spectra  Measured. 

The  spectra  in  parentheses  are  those 
measured  earlier,  without  simultaneous 
irradiation  of  2N2I22A  transistors. 


FAC IL ITT 

SPECTRUM 

FBR , WSMR 

(Glory  Hole) 

FBR , WSMR 

(Free  Field,  50  cm) 

TRIGA,  G . A . 

( J-Tube , Leakage) 

TRIGA,  NORTHROP 

Exposure  Room,  V'  Boral 

TRIGA,  NORTHROP 

Exposure  Room,  V Boral  + 2"  Pb 

ACPR , SANDIA 

ACPR  Core 

FBR.  SANDIA 

Glory  Hole 

FBR,  SANDIA 

Leakage,  Free  Field 

FBR,  SANDIA 

Leakage,  Free  Field,  Burst 

FBR,  SANDIA 

Leakage,  10  cm  Polyethylene 

FBR,  SANDIA 

Leakage,  5 cm  Iron 

APRF , ABERDEEN 

Glory  Hole 

APRF , ABERDEEN 

Leakage,  Free  Field 

APRF,  ABERDEEN 

Leakage,  Free  Field,  Burst 

APRF.  ABERDEEN 

Leakage,  HA  Table*,  Free  Field 

3. 


FOIL,  TRANSISTOR,  anu  TLD 
IRRADIATIONS;  PRESENT  PROGRAM 


3.1  FOIL  IRRADIATIONS 

The  toils  were  stacked  iu  oroer  ot  increasing  threshold 
energy,  Ettl,  so  as  to  minimize  self-shielding  ettects.  The  low- 
est threshold  foil  was  nearest  the  neutron  source.  Three  separate 
stacKS  were  used  as  follows:  ( 1 ) Au,  23jU,  (2)  In,  21oU, 

2^2Th;  and  (3)  Fe,  Ni,  rtg , Al , ana  Zr.  Note  that  the  first  two 
stacks  terminate  on  '•^Np  and  2J‘"lh,  respectively,  since  these 
are  the  thickest  foils  (See  Table  I,  Appendix  A);  the  last  five 
toils,  on  the  other  hand,  are  all  relatively  thin. 

Short  descriptions  ot  the  irradiation  runs  arc  given  in 
Table  2 together  with  the  corresponding  calculated  l-MeV(Si)  equi- 
valent neutron  tluences,  some  related  spectrum  characteristics,  and 
the  transistor  damage  coefficients  obtained  lor  transistors  irra- 
diated at  the  same  time  as  the  foils.  (Discussions  of  the  results 
contained  in  Table  2 are  given  in  sections  e,  7,  and  rt.) 

1 r\  O 

All  foil  stacks  were  placed  in  U spheres  (l.b5  g/ern  ) 
except  in  the  FBK  glory  holes.  There  was  no  room  in  the  glory  hole 
tor  the  sphere,  and  the  sphere  was  not  needeu  there  because 

the  reactor  fuel  provided  even  better  low-energy-neutron  shielding. 
In  the  FUR  leakage-spectrum,  Runs  11  ana  12,  toils  were  placed  both 
Inside  anu  outside  the  l^i3  sphere  to  help  determine  the  effect  ot 

the  l0u  shielding  (scattering)  on  foils  with  high  E. . values. 

lo  tn 

B-Scatterinq  corrections  were  the  order  ot  10-15%,  and  were 

applied  to  all  data  where  the  transistors  were  exposed  outside  the 

shell  while  the  toils  were  located  Inside  the  shell. 

The  use  of  the  shell  is  important  tor  provi- 

2 i s 

ding  an  effective  threshold  tor  U and  Au  toils,  and  also  tor 
drastically  reducing  the  cttects  of  impurities  in  the  toils  that 
leao  to  the  same  activation  nucleus  through  the  much  more  effi- 
cient route  ot  thermal-neutron  activation.  Examples  ot  toil- 
impurity  pairs  are  2*Al(n,or),  21Na(n,y)j  >bFe(n,p), 

^**rtn(n,7)j  2^ttl)(n,t),  21>U(n,t).  The  second  foil  in  each 
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Table  2.  Radiation  Damage  Coefficients 


Run  a 

Dascrlptlon 

V 

_5S 

♦ 

*«9/S 

"1 

tOuA 

2 

lOOuA 

3 

1mA 

4 

10mA 

U 

ACPR  Cor* -8 10 

9 95+12 

1.06 

5.86+31 

4.83-6 

2.79-6 

1 64-6 

0 87-6 

IS 

ACPR  Cora 

ll  44+12 

1.06 

5.86+31 

-.84-6 

2.77-6 

1.59-6 

0.84-6 

■* 

SPR  11  Clory  Ho  la 

2.61+13 

1.16 

5.80+31 

4 83-6 

2.93-6 

1.70-6 

0.90-6 

3 

3PR  II  Laak.  Burse 

3 17+13 

1.35 

6.89+31 

4.57-6 

2.74-6 

1 63-6 

0 35-6 

4 

SPR  II  Laakag* 

4.21+13 

1.42 

6.97+31 

4.20-6 

2.68-6 

1.61-6 

0 . 84-6 

j 

SPR  II  - Ch, 

1.82+13 

1.38 

5.08+31 

4 80-6 

2.82-6 

1.65-6 

0 88-6 

6A 

SPR-II  - Fa 

6 26+13 

0 99 

1.20+32 

- 50-6 

2.72-6 

1.74-6 

0.89-6 

68 

SPR  II  - CH, 

1.41+13 

1.38 

5 08+31 

4.39-6 

2.33-6 

1.65-6 

0 89-6 

7A 

TRICA-Boral-S10 

2.45+13 

1.18 

5 04+31 

4.24-6 

2.57-6 

1.48-6 

0.  ’8-6 

78 

TRICA-Boral 

2.45+13 

1.18 

5.04+31 

4.73-6 

2.85-6 

1.67-6 

0 37-6 

3A 

mGA-Pb-B10 

2.37+13 

1.10 

5 59+31 

4.6—6 

2.73-6 

1.58-6 

0.82-6 

38 

TRIGA-Pb 

2.37+13 

1.10 

5 . 59+31 

- . 63-6 

2.76-6 

1 .60-6 

0 35-6 

9 

APRF -Laakag* 

1.00+14 

1.19 

6.53+31 

3.91-6 

2.40-6 

1 . 6+-6 

0 33-6 

10 

APRF-Glory  Hola 

7.79+13 

1. 19 

6.59+31 

4.36-6 

2.63-6 

1.73-6 

0.38-6 

11 

APRF  on  H.A.  Tab  la 

9.64+13 

1.42 

7.02+31 

3.96-6 

2.50-6 

1.72-6 

0 35-6 

12 

APRF  Laakag* . Burse 

3-9+13 

1.41 

7.08+31 

4.50-6 

2.70-6 

1.57-6 

0 81-6 

Maan 

4.55-6 

2.73-6 

1.65-6 

363-6 

0 

.32-6 

. 14-6 

062-6 

.0*2-6 

(7*.) 

(5*.) 

(4*0 

(5*) 

* 

-22.  i 01-18  M*V) 
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case  is  activated  via  thermal  neutrons,  which  have  a dis- 
proportionately large  reaction  cross  section.  (See  Appendix 
A for  method  of  making  accurate  impurity  corrections  when 
thermal  neutrons  are  shielded  out.) 

3.2  TRANSISTOR  AND  TLD  IRRADIATIONS 

Batches  of  7 to  10  transistors  were  tightly 
wrapped  in  aluminum  foil  and  placed  ■vl  cm  outside  of  the 
ball  containing  the  neutron  threshold-activation  foils. 

The  aluminum  foil  and  the  other  transistor  cans  very  likely 
shorted  out  the  transistor  leads  in  all  cases.  Great  care 
was  taken  to  keep  the  transistors  and  the  neutron  spectrometry 
foils  at  the  same  distance  from  the  reactor  for  free-field 
measurements . 

For  the  case  of  abundant  moderated  neutrons,  such 
as  for  the  ACPR  and  the  Northrup  TRIGA  reactors,  transistors 
were  placed  both  inside  and  outside  the  ^B  sphere  to  deter- 
mine the  radiation  damage  due  to  low  energy  neutrons.  In  all 
cases,  the  measured  transistor  damage  inside  the  ^B  shell 
was  less  than  that  outside  by  about  the  same  factor  (10-157.) 
as  high-energy  threshold  foils  differed  inside  and  out.  From 
this  observation,  it  is  concluded  that  in  no  case  was  any  net 
damage  observed  for  low  energy  neutrons  (i.e.,  neutrons  below 
the  'vlO  keV  neutron  "threshold"  for  the  ^B  shield)  . 

A batch  of  4-5  TLD's  was  placed  next  to  the  trans- 
istors during  each  irradiation  to  measure  total  dose.  The 
TLD's  were  the  thin  CaFjiMn  type,  embedded  in  teflon.  They 
were  wrapped  in  aluminum  foil  of  a total  thickness  of  approxi- 
mately 0.025  cm  so  as  to  shield  out  most  of  the  extraneous- 
electron  dosage. 

Some  of  the  unirradiated  TLD's  of  the  same  batch 
were  exposed  to  a calibrated  pure  ^Co  gamma-ray  field  at  NBS 


by  J.  Humphreys  and  read  with  the  same  equipment  as  the 
reactor-irradiated  TLD's  by  Louise  Miles  of  the  Naval 
Surface  Weapons  Center  (Silver  Springs,  Maryland). 

The  TLD  data  were  used  to  provide  "total  dose" 
corrections  to  the  neutron  damage,  as  discussed  in  Section  8 
be  low . 


4. 


COUNTING 


Counting  the  y-rays  activities  induced  by  the  foil 
irradiation  was  accomplished  for  the  most  part  according  to 
the  ASTM  standards  method  outlined  in  Appendix  B.  Ge(Li) 
Detectors  were  used  exclusively,  both  in  the  irradiation  site 
for  the  short-lived  activities  and  at  the  SAI-La  Jolla  Labora- 
tory for  the  long-lived  ones.  Calibration  of  (and  inter- 
calibration between)  individual  detectors  was  accomplished 
using  an  NBS-prepared  mixed  y-ray  source  covering  the  energy 
range  of  interest  (100  to  ^ 2,000  keV) . 

A number  of  different  foil-to-detector  distances 
were  utilized  in  the  measurement,  depending  on  the  exhibited 
total  activity  of  each  individual  foil  and  the  energy  of  the 
photopeak  of  interest,  so  as  to  minimize  excessive  dead-time 
and  background  effects.  The  detector  efficiency  for  each 
geometry  was  determined  using  the  same  NBS  standard  11-line 
source  for  all  detectors  and  detector  geometries.  A sig- 
nificant addition  to  the  ASTM  Standard  Method  (Appendix  B) 
during  this  work,  was  the  study  of  the  effect  of  y-ray  cas- 
cades on  detector  efficiency;  significantly  in  the  case  of 
the  537  keV-  and  1596  keV  lines  observed  in  all  fission  foils 
used  in  this  work.  This  was  accomplished  by  determining  the 
count-rate  of  the  relevant  photopeak  at  different  distances 
D from  the  "face"  of  the  detector  and  plotting  l//fl" versus  D 
where  N is  the  count  rate.  This  yields  a correction  of  the 
count  rate  due  to  cascades  when  the  foil  being  counted  is 
located  close  to  the  effective  detector  center;  it  also 
locates  the  approximate  effective  detector  center  for  a given 
energy  gamma  ray.  Figures  1 and  2 show  the  results  of  this 
study  for  the  NBS  standard  source  and  the  537  keV  and  1596  keV 
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Figure  L.  Cascade -Gamma  Effect  or.  Determining 
Decector  Efficiency:  Calibration 
Sources . 
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lines  observed  from  the  fission  foils.  Note  in  Figure  1 (NBS 
source)  that  the  failure  of  the  ^03Hg  line  to  pass  through 
the  measured  point  at  O distance  shows  that  a cascade  correction 
is  required.  The  correction  is  significant  and  is  due  to  the 
associated  T1  x rays.  The  cascading  b0Co  1173-  and  1332-keV 
lines  (data  for  the  latter  not  shown  in  Figure  1)  do  not  require 
correction  because  the  efficiency  for  detecting  these  high-energy 
cascades  is  low  enough  to  render  the  effect  insignificant  for  the 
detector  used. 

Note  also  that  for  the  cases  where  cascade  corrections 

are  not  significant  all  points  pass  through  a straight  line  as  ex- 

2 

pecteo  from  the  1/D  dependence  of  N. 

In  Fig.  2,  the  correction  factor  required  for  the  537- 
kev  and  1596-keV  lines  in  the  case  of  the  fission  foils,  was 
determined  to  be  1.0b  when  these  were  counted  at  the  face  of  the 
detector . 


5. 


DATA  ANALYSIS 


The  resulting  gamma-ray  peak  areas  were  further 
corrected  for  decay,  irradiation  and  counting  time  effects. 

By  further  taking  into  account  fission  yield,  branching 
ratios,  detector  efficiency  at  each  relevant  energy  and  foil 
mass,  the  number  of  activations  (fissions)  per  atom  was  deter- 
mined for  input  into  the  SAND  II  code.  The  relevant  formulas 
can  be  found  in  Appendix  B. 

It  may  be  mentioned  that  where  possible  (as  for 
instance  in  the  case  of  fission  foils)  more  than  one  y-ray 
line  was  analyzed  to  provide  a consistency  check  on  the 
results . 


6. 


UNFOLDING 


Because  of  the  limited  data  carried  by  the  threshold 
foils,  especially  below  1 MeV , the  unfolding  code  (SAND  II  in 
this  case)  can  only  be  used  as  a mild-perturbation  type  unfold- 
ing code.  A good  trial  spectrum  is  required;  e.g.,  one 
obtained  from  a simple  transport  calculation  utilizing  the  cor- 
rect neutron  source  and  the  approximate  material  composition 
and  geometry  involved.  For  the  work  carried  out  here,  the  trial 
spectra  {$Cr(E)}  were  as  follows: 

(1)  ACPR-Core  calculation  for  ACPR  and 
TRIGA  reactors  (DOT) 

(2)  Glory  hole  calculation  for  FBR  glorv- 
hole  runs  (IDF) 

(3)  Leakage  calculation  for  FBR  free-field 
runs  (IDF) 

(4)  Polvethvlene-perturbed  SPRII  spectrum 
(ANISN) 

(5)  Iron-perturbed  SPRII  spectrum  (ANISN) 

The  code  carries  out  the  necessary  number  of  itera- 
tions in  perturbing  $tr(E)  until  the  solution  4>(E)  is  obtained 
that  agrees  with  the  input  foil  activations. 

The  results  of  the  0^  iteration  were  used  to  detect 
spurious  foil  activation  values.  These  were  rejected  whenever 
they  were  identified,  and  the  SAND  II  unfolding  code  was  rerun 
without  them. 

The  damage  function  D(E)  was  added  to  our  version  of 
the  SAND  II  code  with  an  energy  mesh  matching  the  620-group 
structure  of  SAND  II  over  the  energy  range  10"^-17.9  MeV,  and 
the  quantities 
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® e ./C , ana  ®eq/S  were  numerically  calculated. 

Note  that  the  sulfur  foil  activations  S are  here 
defined  as  the  activation  rate  per  atom  of  "*2S  initial  target 
material  transmuted  via  the  32S(n,p)32P  reaction  to  32P 
radioactive  nuclei 


S £ NX  /N(  *S) 


3 2 

where  NX  is  the  proauct  of  the  number  of  P nuclei  produced 
and  the  '}2P  aecay  constant,  ana  N(32S)  is  the  number  of  32S 
nuclei  in  the  sulfur  pellet. 

Actual  facility  sulfur  monitor  counts  were  not  usea 
since,  in  one  case,  the  counter  calibration  went  through  three 
different  values  over  a time  span  of  only  a few  months  (Itsu 
Arimura,  private  communication,  1978).  Rather,  S was  obtained 
iteratively,  by  choosing  a value  of  sulfur  activation  that 
agreed  with  the  final  SAND  II  solution  for  0(E).  After  per- 
forming this  exercise  several  times,  it  was  establisned  that 
the  sulfur  activation  could  be  obtained  to  within  +2%  by 
utilizing  the  ratio 


2 . 92*N  X (:>0Ni(n,p):>8Co)  - NX  ( J2S(  n ,p)  32P) 


7. 


NEUTRON  SPECTRA 


Figure  3 shows  the  SANDIA  ACPR  core  spectrum,  as 
measured  by  SAI  and  SANDIA  (SAI's  foils,  SANDIA  detectors 
for  short-lived  foils  and  SAI's  detector  for  long-lived  foils), 
and  Figure  4 shows  the  same  spectrum  measured  by  SANDIA  with 
their  detectors  and  foils.  All  detectors  were  calibrated 
against  an  NBS  11-line  mixed  source  standard.  Both  spectra 
were  unfolded  with  SAND  II,  using  the  same  trial  spectrum. 

This  trial  spectrum  was  calculated  for  the  ACPR  with  a DOT-code 
calculation  by  J.  Odom  of  SANDIA. 

The  "spectral  shape  factor"  lt>e^/<P  is  1.06  for  the 
data  of  Figure  3 and  1.03  for  Figure  4.  These  independent 
runs  show  good  reproducibility  in  terms  of  silicon  radiation 
damage  effectiveness  of  the  neutron  field. 

Figure  5 shows  the  results  of  an  attempt  at  unfold- 
ing the  APRF  spectrum  WITHOUT  the  a priori  knowledge  of  a 
physically  meaningful  trial  spectrum.  The  IUNFOLD  code  of 
SANDIA  was  used  to  obtain  the  spectral  shape  with  only  the 
foil  activation  data.  Above  1 MeV,  where  many  foils  have 
thresholds  and  where  the  spectrum  is  dropping  rapidly,  the 
unfolding  appears  to  be  fair  even  though  the  unfolded  spectrum 
drops  off  a bit  faster  than  the  SAND  II  data  with  adequate  a 
priori  data.  But  below  1 MeV,  where  only  one  threshold  exists 
between  0.01  MeV  and  0.50  MeV,  the  unconstrained  solution  is 
physically  meaningless  (i.e.  not  possible  for  the  source 
spectrum,  the  geometry,  and  the  materials  involved) . Other 
attempts  will  be  made  to  use IUNFOLD  with  a good  trial  spectrum, 
and  much  more  reasonable  results  are  clearly  expected  with 
the  constrained  solution. 
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Calculation  (DOT) 


Energy  (HeV) 

Figure  3.  SANDIA  ACPR  Core.  SAND  II  solution  and  trial  spectrum 
(DOT  calculation).  Unfolded  at  SAI  with  SAI  foils  and 
detector,  but  with  SANDIA  detectors  for  short-lived 
activations . 


UlcuUttoa  (DOT) 


Energy  (MeV) 

Figure  U.  SANDIA  ACPR  Core.  SAND  II  solution  and  trial  spectrum 
(DOT  calculation).  Unfolded  at  SAI . Independent  data 
with  SANDIA  foils  and  Ge (Li) detectors . 


Figures  6 and  7 show  the  Northrop  reactor  spec- 
trum just  outside  the  boral  wall  in  the  exposure  room  adjac- 
ent to  the  reactor  core.  The  shape  factors  $ /$  are  1.18 

€Cj 

and  1.10  for  the  data  of  Figure  6 (V  boral  next  to  the 
reactor  core)  and  Figure  7 (V  boral,  2"  lead). 

In  Figure  8,  the  SPRII  free  field  spectrum  was 
passed  through  a 10-cm  thick  slab  of  polyethylene.  This 
results  in  a highly  moderated  component  below  1 MeV  with 
the  spectral  shape  looking  much  like  that  for  a TRIGA  reac- 
tor. However,  careful  examination  shows  that  above  1 MeV 
the  filtering  has  produced  a harder  spectrum  due  to  the 
falling  off  of  the  C and  H cross  sections  with  increasing 
energy.  This  is  seen  by  the  larger  value  of  d>e<^/4> . being 
equal  to  1.38  which  is  much  larger  than  the  values  of  1.03 
to  1.18  for  the  four  ACPR-core  and  TRIGA- leakage  spectra. 

A greatly  distorted  spectrum  is  shown  in  Figure  9 . 
This  was  achieved  with  a 10-cra  thick  iron  "filter".  The 
inelastic  scattering  cross  section  is  quite  large  for  iron 
at  just  below  1 MeV  and  on  up,  which  depletes  the  spectrum 
at  high  eneriges  and  builds  it  up  a bit  just  above  the  0.2 
MeV  "damage  threshold"  level  where  the  silicon-damage  cross 
section  first  becomes  large.  As  a result,  the  neutrons  above 
the  0.2  MeV  "threshold"  produce  less  damage,  on  the  average. 
This  is  seen  by  the  lowest  value  of  <J,eq/4)  observed  for  all 
these  runs:  i.e.,  $eq/$  ■ 0.99  (see  Table  2 below). 

The  SPRII  and  APRF  glory  hole  spectra,  shown  in 
Figures  10  and  11  respectively,  have  by  far  the  hardest  spec- 
trum below  1 MeV.  Yet,  the  respective  4>eq/4>  values  for  these 
two  spectra  are  1.16  and  1.19;  very  much  like  the  TRIGA 
reactor  spectra.  This  is  due  to  a very  soft  spectrum  above 
1 MeV  which  is  due  to  the  flux  peaking  at  the  core  (glory 
hole)  of  the  reactor.  The  neutrons  reaching  the  glory  hole 
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(Modified  ACPR  Core.  DOT) 


Energy  (MeV) 

Northrop  TRIGA,  Jfc-in  Boral.  SAND  II  solution  and  trial 
spectrum  (1  DF  calculation  modified  by  boral  absorption 
below  10" 2 MeV). 


CeUuWtton  (Modified  ACPR  Core.  DOT) 


Northrop  TRIGA,  %-in  Boral  and  2- in  Lead.  SAND  II  so 

and  trial  spectrum  (1  DF  calculation  modified  by  boral 
absorption  below  10"*  MeV) . 


Calculation  (AMISH) 
SAND  II 


SPR  II  Leakage,  Modified  by  10-cm  Polyethylene.  SAND  II 
solution  and  trial  spectrum  (ANISN  calculation  of  incident 
SPR  II  leakage  neutrons,  from  1 DF  calculation  transmitted 
through  polyethylene) . 
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Figure  10.  SPR  II  Glory  Hole.  Steady  state  SAND  II  solution  and 
trial  spectrum  (1  DF  calculation). 
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have,  on  the  average,  passed  through  the  least  material  of 
any  of  the  neutron  fields  encountered  in  this  series. 


f 

The  SPRII  and  APRF  FBR  reactor  leakage  spectra 
shown  in  Figures  3,  13,  14,  15,  and  16,  yielded  respective 
^eq'^  1.42,  1.19,  1.42,  and  1.41.  The  value  for 

SPRII  steady  state  spectrum  of  Figure  13,  $eq/$  ■ 1.42,  is 
only  slightly  higher  than  the  value  of  1.35  for  the  fast 
burst  spectrum  of  Figure  12.  On  the  other  hand,  the  value 
for  the  APRF  steady  state  spectrum  of  Figure  is  only  1.19, 
which  is  lower  than  the  value  of  " 1.41  for  the  fast- 

burst  counterpart.  This  value  is  nearly  the  same  as  for  the 
leakage  spectrum  measured  at  the  hardness  assurance  table 
(Figure  15).  Thus,  the  measurements  of  fast-burst  spectra 
show  no  clear  difference  in  hardness  or  silicon-damage  capa- 
bility from  the  steady  state  counterparts,  either  in  ^gq/41  or 
or  in  detailed  spectral  shape. 

It  should  be  noted  that  the  above  leakage  spectra 
were  used  as  Inputs  to  ANISN  calculations  for  the  spectral 
shape  of  neutrons  emerging  from  the  10  cm  thick  slab  of  poly- 
ethylene and  the  5-cm  slab  of  iron:  the  ANISN  (emergent) 
spectra  were  in  each  case  used  as  the  trial  spectrum  $ for 
the  polyethylene  and  iron  filter. 

The  foil  activation  data  used  to  obtain  the  spectra 
of  Figures  3-16  are  given  in  Table  3.  Here,  the  run  numbers 
correspond  to  those  described  in  Table  2. 

The  unfolding  outputs  corresponding  to  the  spectra 
of  Figures  3,  4,  and  6-16  are  given  in  Table  4,  along  with 
the  trial  spectrum  used  for  each  SAND  II  unfolding  operation. 

Again,  the  run  numbers  correspond  to  those  described  in  Table 
2,  except  that  Run  13  in  Table  4 corresponds  to  the  indepen- 
dent ACPR  core  data  (Figure  4)  obtained  from  Hugh  Scott  of 
Sandia.  Although  not  stated  in  Table  2,  all  FBR  leakage 
spectra  were  with  the  FBR  ^1%  m above  the  floor  except  for 
Run  12  ( 6m). 
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Table  3.  Activations  Used  for  SAND  II  Unfolding 


Table  4A.  Tabulated  SAND  11  values  of  E *$  . . and  E*<6(E).  E is  in  MeV 
and  $(E)  is  in  neutrons  per  car  per  MeV.  cJlE)  is  normalized  to  unit 
fluence  above  .01  MeV.  The  multiplier  given  in  the  last  row  may  be  used 
to  obtain  the  absolute  fluence  for  each  corresponding  run. 


RUN  1 

RUN  13 

RUN  7 

RUN  8 

En(H.V) 

EV 

E*e(E) 

E*e(E) 

E-tr1 2 

E*6(E) 

£*♦(£) 

1.0-4 

4.95-2 

3.60-2 

6.88-2 

1.60-2 

1.50-2 

1.54-2 

1.5-4 

5.00-2 

3.84-2 

6.78-2 

1.77-2 

1.66-2 

1.70-2 

2.0-4 

5.04-2 

4.14-2 

6.73-2 

1.93-2 

1.80-2 

1.85-2 

3.0-4 

4.96-2 

4.20-2 

6.34-2 

2.13-2 

1.99-2 

2.05-2 

4.0-4 

5.15-2 

4.50-2 

6.49-2 

2.31-2 

2.16-2 

2.22-2 

5.0-4 

5.20-2 

4.60-2 

6.09-2 

2.47-2 

2.31-2 

2.37-2 

6.0-4 

5.20-2 

4.57-2 

6.53-2 

2.60-2 

2.43-2 

2.49-2 

8.0-4 

5.23-2 

4.57-2 

6.24-2 

2.81-2 

2.63-2 

2.70-2 

1.0-3 

5.25-2 

4.64-2 

6.143-2 

2.98-2 

2.79-2 

2.86-2 

1.5-3 

5.25-2 

4.67-2 

6.18-2 

3.24-2 

3.03-2 

3.11-2 

2.0-3 

5.41-2 

4.87-2 

6.30-2 

3.47-2 

3.25-2 

3.34-2 

3.0-3 

5.11-2 

4.79-2 

5.86-2 

3.77-2 

3.52-2 

3.62-2 

4.0-3 

5.40-2 

5.06-2 

6.20-2 

3.96-2 

3.72-2 

3.82-2 

5.0-3 

5.48-2 

5.16-2 

6.27-2 

4.16-2 

3.90-2 

4.00-2 

6.0-3 

5.52-2 

5.21-2 

6.31-2 

4.X-2 

4.03-2 

4.14-2 

8.0-3 

5.57-2 

5.31-2 

6.34-2 

4.54-2 

4.25-2 

4.36-2 

1.0-2 

5.67-2 

5.51-2 

6.44-2 

4.72-2 

4.42-2 

4.54-2 

1.5-2 

5.80-2 

5.70-2 

6.57-2 

5.00-2 

4.68-2 

4.81-2 

2.0-2 

6.03-2 

5.98-2 

6.82-2 

5.27-2 

4.93-2 

5.07-2 

3.0-2 

6.38-2 

6.41-2 

7.18-2 

5.70-2 

5.34-2 

5.49-2 

4.0-2 

6.71-2 

6.80-2 

7.53-2 

6.71-2 

6.29-2 

6.47-2 

5.0-2 

7.09-2 

7.22-2 

7.94-2 

7.09-2 

6.65-2 

6.84-2 

6.0-2 

7.39-2 

7.57-2 

8.27-2 

7.39-2 

6.93-2 

7.12-2 

8.0-2 

8.13-2 

8.38-2 

9.08-2 

8.13-2 

7.62-2 

7.83-2 

1.0-1 

8.76-2 

9.12-2 

9.93-2 

8.76-2 

8.22-2 

8.45-2 

1.5-1 

1.02-1 

1.07-1 

1.13-1 

1.02-1 

9.62-2 

9.89-J 

2.0-1 

1.17-1 

1.24-1 

1.30-1 

1.17-1 

1.11-1 

1.14-1 

2.4-1 

1.29-1 

1.37-1 

1.42-1 

1.29-1 

1.22-1 

1.26-1 

3.0-1 

1.46-1 

1.56-1 

1.60-1 

1.46-1 

1.39-1 

1.44-1 

3.4-1 

1.57-1 

1.69-1 

1.72-1 

1.57-1 

1.52-1 

1.57-1 

4.0-1 

1.61-1 

1.75-1 

1.73-1 

1.61-1 

1.60-1 

1.67-1 

4.5-1 

1.69-1 

1.86-1 

1.80-1 

1.69-1 

1.73-1 

1.81-1 

5.0-1 

1.86-1 

2.06-1 

1.97-1 

1.86-1 

1.95-1 

2.05-1 

6.0-1 

2.08-1 

2.34-1 

2.17-1 

2.08-1 

2.26-1 

2.39-1 

7.2-1 

2.31-1 

2.63-1 

2.40-1 

2.31-1 

2.58-1 

2.72-1 

8.0-1 

2.41-1 

2.76-1 

2.54-1 

2.41-1 

2.73-1 

2.87-1 

9.2-1 

2.36-1 

2.71-1 

2.48-1 

2.36-1 

2.71-1 

2.84-1 

1.0 

2.34-1 

2.69-1 

2.48-1 

2.34-1 

2.72-1 

2.82-1 

1.2 

2.51-1 

2.90-1 

2.74-1 

2.51-1 

2.99-1 

3. 0S-1 

1.5 

2.57-1 

3.05-1 

3.14-1 

2.57-1 

3.25-1 

3.22-1 

1.7 

2.58-1 

3.03-1 

3.14-1 

2.58-1 

3.30-1 

3.20-1 

2.0 

2.47-1 

2.92-1 

3.00-1 

2.47-1 

3.19-1 

3.07-1 

2.5 

2.32-1 

2.72-1 

2.60-1 

2.32-1 

3.02-1 

2.87-1 

3.0 

1.98-1 

2.32-1 

2.18-1 

1.98-1 

2.61-1 

2.45-1 

3.5 

1.61-1 

1.89-1 

1.81-1 

1.61-1 

2.13-1 

2.00-1 

4.0 

1.35-1 

1.56-1 

1.44-1 

1.35-1 

1.80-1 

1.66-1 

4.5 

1.13-1 

1.30-1 

1.19-1 

1.13-1 

1.51-1 

1.39-1 

5.0 

9.23-2 

1.06-1 

9.51-2 

9.23-2 

1.24-1 

1.13-1 

6.0 

5.85-2 

6.16-2 

5.51-2 

5.85-2 

8.08-2 

6.84-2 

7.0 

3.55-2 

3.14-2 

3.08-2 

3.50-2 

5.03-2 

3.71-2 

8.0 

2.10-2 

1.75-2 

1.73-2 

2.10-2 

2.97-2 

2.10-2 

9.0 

1.22-2 

1.01-2 

1.01-2 

1.22-2 

1.77-2 

1.22-2 

10.0 

6.94-3 

5.72-3 

5.70-3 

6.94-3 

1.01-2 

6.94-3 

12.0 

2.12-3 

1.74-3 

1.73-3 

2.12-3 

3.11-3 

2.11-3 

15.0 

3.16-4 

2.34-4 

1.97-4 

3.16-4 

6.23-4 

5.14-4 

17.0 

8.07-5 

5.96-5 

5.01-5 

8.07-5 

1.60-4 

1.33-4 

(Multiplier) 

* 

(9.38*12) 

(5.05*14) 

“ 

(2.18*13) 

(2.16*13) 

1.  *tr  • ACRR  cor*  (DOT) 

2.  *tr  • ACPR  cor*  modified  below  3 x 10"c  MeV  by  bora)  ittcnuatlon 

NOTE:  Run  numbers  correspond  to  Table  2 except  Run  13  Is  SANOIA  Independent  foil  data  for  tC’R 
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Table  4B.  Tabulated  SAND  II  values  of  E*0 


and  E* 0(E).  E is  in  MeV 
is  normalized  to  unit 


Jues  of  E*0fc  • , 

and  0(E)  is  in  neutrons  per  cin  per  MeV.  0?E) 
fluence  above  .01  MeV.  The  multiplier  given  in  the  last  row  may  be  used 
to  obtain  the  absolute  fluence  for  each  corresponding  run. 


RUN  3 

RUN  4 

RUN  9 

RUN  11 

RUN  12 

En(H.V) 

E*t  3 
t *tp 

E*b(E) 

E«>(E) 

E*4(E) 

E +4(E) 

E**tr3 4 

E*P(E) 

1.0-4 

1.29-8 

3.88-4 

3.50-4 

4.53-4 

3.28-4 

3.80-9 

1.12-4 

1.5 

1.32-8 

4.00-4 

3.60-4 

4.73-4 

3.27-4 

4.24-9 

1.16-4 

2.0 

1.40-8 

4.27-4 

3.80-4 

4.98-4 

3.47-4 

4.39-9 

1.23-4 

3.0 

1.48-8 

4.62-4 

4.04-4 

5.33-4 

3.63-4 

4.74-9 

1.31-4 

4.0 

1.56-8 

4.94-4 

4.25-4 

5.57-4 

3.88-4 

5.07-9 

1.39-4 

5.0 

1.63-8 

5.21-4 

4.43-4 

5.68-4 

4.23-4 

5.38-9 

1.50-4 

6.0 

1.68-8 

5.48-4 

4.58-4 

6.01-4 

4.16-4 

5.63-9 

1.52-4 

8.0 

1.88-8 

6.27-4 

5.11-4 

6.75-4 

4.58-4 

6.58-9 

1.70-4 

1.0-3 

2.05-8 

6.94-4 

5.58-4 

7.34-4 

5.04-4 

7.18-9 

1.89-4 

1.5 

2.39-8 

8.27-4 

6.51-4 

8.55-4 

5.91-4 

8.60-9 

2.26-4 

2.0 

2.69-8 

9.47-4 

7.33-4 

9.58-4 

6.72-4 

9.82-9 

2.66-4 

3.0- 

3.30-8 

1.16-3 

9.00-4 

1.15-3 

8.56-4 

1.35-8 

3.72-4 

4.0 

3.85-8 

1.37-3 

1.05-3 

1.35-3 

9.95-4 

1.67-8 

4.45-4 

5.0 

4.34-8 

1.55-3 

1.19-3 

1.52-3 

1.13-3 

2.00-8 

5.22-4 

6.0 

5.01-8 

1.80-3 

1.37-3 

1.75-3 

1.31-3 

2.51-8 

6.23-4 

8.0 

6.35-8 

2.27-3 

1.73-3 

2.20-3 

1.67-3 

3.49-8 

8.49-4 

1.0-2 

7.63-8 

2.71-3 

2.08-3 

2.62-3 

2.03-3 

4.58-8 

1.12-3 

1.5 

1.21-7 

4.31-3 

3.32-3 

4.14-3 

3.27-3 

7.74-8 

1.96-3 

2.0 

1.72-7 

6.10-3 

4.71-3 

5.84-3 

4.66-3 

1.14-7 

3.00-3 

3.0 

2.87-7 

1.01-2 

7.85-3 

9.64-3 

7.84-3 

2.01-7 

5.47-3 

4.0 

4.13-7 

1.44-2 

1.13-2 

1.38-2 

1.14-7 

3.06-7 

8.35-3 

5.0 

5.51-7 

1.91-2 

1.51-2 

1.83-2 

1.52-2 

4.44-7 

1.15-2 

6.0 

7.14-7 

2.47-2 

1.95-2 

2.37-2 

1.98-2 

5.96-7 

1.55-2 

8.0 

1.08-6 

3.71-2 

2.95-2 

3.55-2 

3.00-2 

9.50-7 

2.43-2 

1.0-1 

1.50-6 

5.13-2 

4.11-2 

4.92-2 

4.21-2 

1.37-6 

3.59-2 

1.5 

2.71-6 

9.20-2 

7.42-2 

8.82-2 

7.62-2 

2.55-6 

6.82-2 

2.0 

3.97-6 

1.34-1 

1.09-1 

1.29-1 

1.12-1 

3.65-6 

1.03-1 

2.4 

5.10-6 

1.71-1 

1.40-1 

1.64-1 

1.44-1 

5.10-6 

1.36-1 

3.0 

6.68-6 

2.23-1 

1.83-1 

2.14-1 

1.90-1 

6.68-6 

1.85-1 

3.4 

7.69-6 

2.55-1 

2.10-1 

2.44-1 

2.18-1 

7.69-6 

2.17-1 

4.0 

8.72-6 

2.85-1 

2.38-1 

2.71-1 

2.48-1 

8.72-6 

2.52-1 

4.5-1 

1 .51-6 

3.07-1 

2.58-1 

2.90-1 

2.70-1 

9.51-6 

2.80-1 

5.0 

1.03-5 

3.28-1 

2.80-1 

3.09-1 

2.91-1 

1.03-5 

3.06-1 

6.0 

1.15-5 

3.50-1 

3.10-1 

3.29-1 

3.23-1 

1.15-5 

3.45-1 

7.2 

1.22-5 

3.50-1 

3.31-1 

3.36-1 

3.44-1 

1.22-5 

3.68-1 

8.0 

1.26-5 

3.45-1 

3.41-1 

3.40-1 

3.53-1 

1.26-5 

3.71-1 

9.2 

1.30-5 

3.38-1 

3.52-1 

3.43-1 

3.61-1 

1.30-5 

3.74-1 

1.0 

1.28-5 

3.18-1 

3.49-1 

3.35-1 

3.56-1 

1.28-5 

3.57-1 

1.2 

1.31-5 

3.09-1 

3.60-1 

3.37-1 

3.63-1 

1.31-5 

3.51-1 

1.5 

1.30-5 

3.03-1 

3.62-1 

3.29-1 

3.55-1 

1.30-5 

3.52-1 

1.7 

1.26-5 

2.81-1 

3.53-1 

3.16-1 

3.39-1 

1.26-5 

3.25-1 

2.0 

1.20-5 

2.73-1 

3.37-1 

3.00-1 

3.23-1 

1.20-5 

3.12-1 

2.5 

1.07-5 

2.52-1 

3.04-1 

2.69-1 

2.88-1 

1.07-5 

2.79-1 

3.0 

9.17-6 

2.33-1 

2.64-1 

2.33-1 

2.48-1 

9.17-6 

2.49-1 

3.5 

7.70-6 

1.95-1 

2.21-1 

1.95-1 

2.11-1 

7.70-6 

2.21-1 

4.0 

6.38-6 

1.71-1 

1.85-1 

1.63-1 

1.75-1 

6.38-6 

1.89-1 

4.5 

5.08-6 

1.39-1 

1.47-1 

1.30-1 

1.39-1 

5.08-6 

1.53-1 

5.0 

4.08-6 

1.11-1 

1.17-1 

1.05-1 

1.11-1 

4.08-6 

1.22-1 

6.0 

2.47-6 

6.05-2 

6.62-2 

6.11-2 

6.10-2 

2.47-6 

6.59-2 

7.0 

1.40-6 

3.39-2 

3.63-2 

3.24-2 

3.18-2 

1.40-6 

3.43-2 

8.0 

7.79-7 

1.86-2 

2.01-2 

1.76-2 

1.69-2 

7.79-7 

.1.83-2 

9.0 

4.01-7 

1.02-2 

1.03-2 

9.11-3 

8.69-3 

4.01-7 

9.64-3 

10.0 

2.14-7 

5.53-3 

5.49-3 

4.85-3 

4.61-3 

2.14-7 

5.15-3 

12.0 

5.91-8 

1.54-3 

1.51-3 

1.34-3 

1.26-3 

5.91-8 

1.42-3 

15.0 

8.74-9 

2.25-4 

2.44-4 

1.46-4 

2.08-4 

8.74-9 

2.10-4 

17.0 

2.23-9 

5.60-5 

6.23-5 

3.71-5 

5.32-5 

2.23-9 

5.36-5 

(Multiplier) 

* 

(2.35+13) 

(2.96+13) 

(8.41+13) 

(7.86+13) 

“ 

(2.46+13) 

3.  *tr  • FBR  leakage,  FBR  1-1/2  m above  floor  (10F  calculation) 

4.  atr  • FBR  leakage,  FBR  6 m above  floor  (IDF  calculation) 
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Table  4C.  Tabulated  SAND  II  values  of  E*®trial  and  E*0(E).  E is  in  MeV 
and  ( E ) is  in  neutrons  per  cm^  per  MeV.  pTE)  is  normalized  to  unit 
tluence  aoove  .01  MeV.  The  multiplier  ^iven  in  the  last  row  may  be  used 
to  obtain  the  absolute  tluence  for  each  cor responu ing  run. 


RUN  5 

RUN  6 

RUN  2 

RUN  10 

En<""> 

E**tr5 

E*4(E) 

E-tr6 

E*a(E) 

E.*tr7 

£*♦( E ) 

E+*(E) 

1.0-4 

2.34-1 

2.90-2 

6.71-3 

2.13-4 

2.57-11 

1.62-8 

1.15-8 

1.5 

2.27-1 

2.86-2 

7.14-3 

2.26-4 

5.82-11 

3.74-8 

2.61-8 

7.0 

2.26-1 

3.00-2 

7.54-3 

2.39-4 

1.03-10 

6.57-8 

4.59-8 

3.0 

2.30-1 

3.20-2 

8.95-3 

2.84-4 

2.33-10 

1 . 50-7 

1.04-7 

4.0 

2.39-1 

3.59-2 

1.06-2 

3.37-4 

4.13-10 

2.65-7 

1 .85-7 

S.O 

2.47-1 

4.01-2 

1.22-2 

3.66-4 

6.42-10 

3.98-7 

2.87-7 

6.0 

2.53-1 

3.97-2 

1.36-2 

4.30-4 

9.29-10 

5.95-7 

4.16-7 

8.0 

2.63-1 

4.75-2 

1.61-2 

5.11-4 

1 .92-9 

1.24-6 

8.57-7 

1.0-3 

2.58-1 

4.20-2 

1.96-2 

6.21-4 

3.38-9 

2.17-6 

1.51-6 

1.5 

2.47-1 

4.09-2 

2.82-2 

8.94-4 

9.60-9 

6.16-6 

4.30-6 

2.0 

2.43-1 

4.12-2 

3.65-2 

1.16-3 

1.97-8 

1.26-5 

8.82-6 

3.0 

2.54-1 

4.38-2 

5.09-2 

1.61-3 

5.42-8 

3.37-5 

2.43-5 

4.0 

2.63-1 

4.55-2 

6.44-2 

2.04-3 

1.10-7 

6.82-5 

4.91-5 

5.0 

2.72-1 

4.71-2 

7.75-2 

2.45-3 

1.88-7 

1.17-4 

8.42-5 

6.0 

2.77-1 

4.80-2 

9.08-2 

2.88-3 

2.92-7 

1.81-4 

1 .31-4 

8.0 
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(1.33+13) 
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(6.52*13) 

“ 

(2.19+13) 

(6.53+13) 

5.  $tr  • ANICN  calculation  of  FBR  leakage  neutrons  penetrating  10  cm  polyethylene 

6.  *tr  ■ ANICN  calculation  of  FBR  leakage  neutrons  penetrating  5 cm  Iron 

7.  *tr  ■ Glory  Hole  spectrum  (IDF  calculation) 
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Calculation  (IDF) 
SAKO  11 


Figure  12.  SPR  II  Leakage,  Fast  Burst.  Sand  II  solution  and  trial 
spectrum  (1  DF  calculation) . 
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Leakage  Spectrum,  Steady  State  . 

II  solution  and  trial  spectrum  (1  DF  calculation) . 
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Figure  15.  APRF  Leakage  on  Hardness  Assurance  Table,  Steady  State. 


CeUuutlon  (Modified  1 0F  teak  Age) 


Leakage,  Burst,  APRF  6 m Above  Floor.  SAND  II  solution 
rial  spectrum  (1  DF  calculation) . 


8.  TRANSISTOR  DAMAGE  COEFFICIENT  MEASUREMENTS 


In  order  to  determine  how  well  the  standard  method 
for  determination  of  1 MeV  equivalent  neutron  fluence  works 
for  a i""actical  situation,  sets  of  2N2222A  transistors  were 
exposed  with  each  run  and  resultant  neutron  damage  coeffi- 
cients compared.  Test  samples  were  drawn  from  a lot  of  200 
cotnmerical  devices  purchased  from  one  vendor  and  representing 
two  different  data  codes.  All  200  units  were  individually 


marked  and  characterized  for  D.C.  hpp  and  base  transit  time 
tfe.  hpp  was  measured  on  an  automated  test  set  at  room  tem- 
perature with  VCE  - 10V  and  Ic  - 10uA,  lOOuA,  1mA  and  10mA. 


(Peak  hpE  for  these  devices  occurs  near  10mA) . t^  Was  extra- 
polated from  plots  of  (2nf^)  ^ vs  (Ip)*^.  The  f^,  values  were 
obtained  from  hpE  measurements  made  on  a commercial  test  set 


at  20  MHz. 


An  initial  experiment  was  run  on  the  WSMR  Fast  Burst 
Reactor  to  determine  the  effect  of  a t^  screen  on  the  distri- 
bution of  neutron  damage  coefficients,  K,  where 


1 _ 1 
(hVEH  (bpE^ 


cbK* 


Ten  devices  with  .586  nSec  < t^  < .623  nSec  and  ten 
devices  with  .551  nSec  i t^  s.  .839  nSec  were  exposed  to  the 
same  neutron  fluence.  The  average  ratio  of  the  damage  coef- 
ficient standard  deviation  to  mean  for  the  four  collector  cur- 
rents was  6.37.  for  the  sample  with  the  wider  t^  range,  and 
8.37,  for  the  sample  with  tighter  t^  range.  Therefore,  no 
attempt  was  made  to  preselect  units  on  the  basis  of  a t^ 
screen . 


38 


A minimum  of  seven  2N2222A  units  was  exposed  with 
each  set  of  activation  foils  and  in  four  cases  the  units  were 
both  inside  and  outside  the  Boron-10  ball. 

All  units  were  exposed  passively  and  each  batch  was 
tightly  grouped  to  minimize  spatial  effects.  Post  irradiation 
hpg.  measurements  were  made  about  one  month  after  the  final 
exposure.  Measurements  were  made  at  room  temperature  on  the 
same  test  set  as  the  preirradiation  hp£  measurements.  Com- 
parison of  control  samples  indicated  a repeatability  of  the 
hpp  measurements  to  within  27., 

A significant  exposure  to  total  ionizing  (gamma) 
radiation  resulted  with  some  of  the  neutron  exposures . There- 
fore, a correction  was  made  for  total  (gamma)  dose  effects. 

The  increase  in  base  current  from  both  neutron  effects  and 
total  dose  effects  is  additive,  thus: 

(Al/hpE)ToTAL  “ Cb  K * + ^Al/hFE)y 


The  added  (Al/hpp)  term  can  be  subtracted  out  as 
a correction  factor  to  give  proper  K values.  In  order  to 
determine  the  (Al/hpp)Y  term  for  each  set  of  devices,  expos- 
ures were  made  on  a Co^O  source.  A total  of  36  additional 
2N2222A's  from  the  same  lot  as  used  for  the  neutron  expos- 
urcs  were  passively  exposed  to  3x10  , 10  , 3x10  , and  103 
rad(Si) , 9 devices  at  each  level.  The  (Al/hpp)^  was  deter- 
mined for  each  device  at  V^p  - 10V  and  Ic  - lOuA,  lOOuA, 

1mA,  and  10mA.  A plot  of  the  median  Al/hpp  vs  y was  made  for 
each  collector  current  level  and  the  value  of  (Al/hpp)  cor- 
responding to  the  total  dose  level  of  each  neutron  exposure 


was  subtracted  from  (Al/hpp) <p0TAL  to  8ive  c^e  correcte<*  K 


values.  The  maximum  coefficient  of  variation  (- 


for  the 


(Al/hpp) Y values  ranged  from  ^50"  at  10  A to  ^307,  at  10mA. 


I 
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CaF2  TLDS  were  used  with  each  neutron  exposure  to  determine 
the  total  ionizing  dose.  The  correction  factors  rangeo  from 
less  than  1%  of  1/hpE*TOTAL  tor  most  ot  the  fast-burst 
reactor  exposures,  up  to  10%  for  exposure  #7  which  was  per- 
formed on  the  TRIGA  Marx  F benind  1/4"  Uoral.  The  total  dose 

4 

for  this  exposure  was  9.5x10  rad(Si). 

Table  2,  shown  on  page  fa  of  Section  3,  presents 
the  results  of  the  K*  measurements  for  the  2N2222A  K's  mea- 
sured at  Ic  ■ lOi.A,  lOOuA,  1mA  ana  10mA.  The  standara  devia- 
tion for  each  set  of  7 to  10  units  is  also  given.  Here, 

K*  * (M/hF£)  / (tb) 

The  values  ot  Kt  were  divided  by  where  the  *•>  were 

obtained  from  the  spectrum  measurements  of  0(E),  toldeo  in  with 
D(E)  as  per  Equation  1.  Tne  values  of  K ^ through  k4  are  given  in 
Table  2,  along  with  and  ^eq/ ' lor  each  spectrum.  Here,  Kj  K 

(lOuA),  K2  * K (lOOuA),  K j - K(lmA),  ano  K4  (10mA). 

Note  that  for  K2  through  K4  the  standard  deviation 
on  K is  4-5%  (Kj,  tor  10uA,  is  perhaps  of  little  interest)  inui- 
cating  that  D( E ) represents  the  damage  function  for  2N2222A's 
over  a wide  range  of  both  unperturbed  and  "filtered"  reactor 
spectra.  This  should  make  it  possible  quantitatively  to  inter- 
relate the  radiation  damage  work  with,  say,  a TRIGA  reactor  and 
an  FUR  leakage  spectrum.  The  value  of  £ must  easily  be  measur- 
able for  each  irradiation,  as  is  oiscusseo  in  Section  9 below. 

Since  the  observed  damage  coefficients  are  seen  to 
be  nearly  constant  over  a range  ot  spectral  shapes  where  o /c, 
the  damage  per  neutron  above  0.01  MeV,  was  made  to  range  from 
0.99  to  1.42,  these  data  help  establish,  although  not  con- 
clusively, the  usefulness  ot  D(E)  ano  0 in  quantitatively 


40 


relating  the  damage  induced  by  a basic  (or  modified)  reactor 
spectrum  to  the  damage  expected  of  a threat  spectrum.  One 
would  expect  a good  correlation  at  least  in  the  case  of  larger 
volume  transistors,  where  volume  recombination  effects  predomin- 
ate in  the  damage  mechanism  involved. 


9.  PRACTICAL  METHOD  FOR  MEASUREMENT  OF  4> 

eq 

Table  2 lists  values  of  # /S  for  all  the  measure- 

eq 

ments  of  basic  and  modified  spectra  carried  out  at  the  SANDIA 

SPRII  and  ACPR,  the  APRF  FBR  and  the  Northrop  TRIGA  facilities. 

In  order  to  obtain  # for  any  run  (and  thereby  obtain  a true 

measure  of  damagability) , one  measures  the  absolute  sulfur 

activation  S and  multiplies  by  $ /S  to  get  $ . Here  S is 

ec\  32 

the  absolute  disintegration  rate  (NX)  of  P,  corrected  back 
to  the  time  at  the  end  of  the  irradiation. 

Two  caveats  must  be  interjected  at  this  point.  The 

first  is  that  an  absolute  calibration  (rigorous,  yet  possibly 

" 32 

very  simple)  in  units  of  disintegration  rate  per  S nucleus 

is  needed.  The  second  is  that  the  $ /S  deduced  from  the 

eq 

earlier  work  presented  in  Appendix  A should  not  be  used  for 
the  WSMR  FBR  and  the  GA  TRIGA.  Since  those  data  were  obtained 
without  the  benefits  of  the  recent  improvements  in  foil  counting 
(the  cascade  correction),  and  since  all  reactor  spectra  of  a given 
type  (e.g.,  all  FBR’s  of  roughly  the  sane  dimensions)  are  apprecia- 
bly closer  together  than  the  error  arising  from  the  cascade  correc- 
tion. it  is  by  far  more  accurate  to  use  the  SPRII  leakage  spectrum 

: /S  (Table  1)  for  the  WS'TR  FBR  than  to  obtain  it  directlv 

eq 

from  the  earlier  work  on  the  WSMR  FBR  of  Appendix  A. 

32  32 

The  absolute  P activity  per  S nucleus  can  be 
readily  determined  simply  by  the  method  alluded  to  in  the  text 
for  getting  the  sulfur  activation,  or  by  the  method  suggested 
to  us  by  H.  Wright  of  WSMR.  In  the  latter  method,  the  sulfur 
foils  and  some  aluminum  foils  are  irradiated  in  the  same  14  MeV 

->4 

neutron  flux.  The  aluminum  foils  are  counted  for  the  2.75  MeV  “ Na 
line  with  either  a Ge(Li)  or  Nal  detector  calibrated  against 
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some  NBS  standard  foils  (thorium  and  yttrium),  and  the  absol- 
ute 14  MeV  neutron  flux  is  calculated  from  this  measurement 
with  the  use  of  the  Al(n,a)z  Na  cross  section  at  14  MeV. 
With  a knowledge  of  the  14  MeV  flux  and  the  ^S(n,p)^P 

reaction  cross  section  at  14  MeV,  the  beta  detector  is  cali- 

27  12 

brated  absolutely.  Both  the  Al(n.cx)  and  S(n,p)  cross 
sections  are  known  with  reasonably  good  accuracy  at  14  MeV. 

With  the  method  developed  in  this  program  and 
alluded  to  in  the  text,  one  simply  activates  a nickel  foil 
and  a sulfur  foil  together  in  the  same  reactor  neutron  field 
and  counts  the  absolute  nickel  activation  rate  with  a cali- 
brated Nal  or  Ge(Li)  detector  (see  Appendix  B) . The  sulfur 
32  32 

( P(  activation  rate  per  S-target  atom  is  then  taken  as 
2.92  x the  ^®Ni  activation  rate  per  ^®Ni  foil  atom.  This 
approach  was  found  to  be  good  to  an  estimated  accuracy  of 
i27.  in  the  refolding  calculations  carried  out  with  the  SANDII 
code . 
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The  accuracy  of  unfolding  neutron  spectra  with  threshold  activation  foils  was  investigated 
by  comparing  the  unfolded  spectra  with  "known"  spectra  obtained  from  calculations  verified 
with  more  accurate  spectrometry  techniques.  Because  of  the  lack  of  spectral  data  below  1 MeV 
(only  one  useful  activation  threshold  exists  between  0.01  and  1 MeV),  reliable  results  are 
generally  obtained  only  with  the  a priori  constraint  of  a physically  complete  trial  spectrum.  A 
family  of  trial  spectra  useful  for  reactor  spectra  land  any  nuclear  "evaporation  ' spectra,  such  as 
from  high-energy  (y.n)  reactions/  has  been  formulated.  The  sensitivity  of  the  solution  <dE)  to 
the  specific  unfolding  code  used,  to  the  threshold-foil  cross  sections  and  related  nuclear  data, 
and  to  foil-activation  measurement  accuracy  was  also  investigated. 


I.  INTRODUCTION 


This  paper  summarizes  the  results  of  a pro- 
gram in  which  a standard  neutron  spectrometry 


method  was  selected  and  evaluated  as  to  accuracy. 
The  program  was  carried  out  in  support  of  a set 
of  draft  American  Society  for  Testing  and  Mate- 
rials (ASTM)  Standards  for  characterizing  neutron 
fields  in  terms  of  radiation  damage  effectiveness 
to  silicon  semiconductor  electronic  components. 
Such  a characterization  is  of  interest  in  planning 
irradiation  schedules  for  parts  screening  or  sam- 
ple specification  tests  and  in  comparing  the  re- 
sults of  radiation  damage  studies  carried  out  in 
different  neutron  fields- 

The  method  provides  the  neutron  spectral  in- 
tensity e(£)  for  general  use,  such  as  calculating 
radiation  damage  to  materials  other  than  silicon 
where  the  damage  versus  neutron  energy,  D(E),  is 
known  or  is  calculable  from  known  neutron  cross 
sections. 

O.  SELECTION  OF  THE  THRESHOLD-FOIL  METHOD 

Several  candidate  methods  were  considered  for 
characterizing  a neutron  field  in  terms  of  bulk 

‘Present  address:  Science  Applications,  Inc..  PO.  Box  2351, 
La  Jolla.  California  92038. 


silicon  radiation  damage  effectiveness.  These 
included  time-of-flight  (TOF)  spectrometry,1’2  the 
proton-recoil  scintillator,1'1  the  proton-recoil  • 
proportional  counter, *'*  sHe  and  *Li  coincidence- 
measurement  techniques  (sandwich  counters),9 
the  proton-recoil  telescope,  nuclear  emulsions, 
threshold-foil-activation  spectrometry,  and  the 
use  of  silicon  itself  as  a dosimeter.  The  latter, 
although  potentially  the  most  straightforward,  was 
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eliminated  because  of  difficulties  in  measuring 
the  number  of  dislocation  centers.  The  silicon 
would  have  to  be  kept  at  cryogenic  temperatures 
during  and  after  irradiation,  and  damage  annealing 
from  all  sources  would  have  to  be  correctly 
accounted  (or  The  annealing  can  vary  with  the 
type  of  device  used,  thi  cutt*?*  drain,  tempera* 
ture.  gamma -ray  dose  rate,  and  even  the  neutron 
dose  rate- 

The  threshold -foil  spectrometer  waa  chosen 
for  its  versatility  in  measuring  both  pulsed  and 
steady -state  reactor  spectra,  its  wide  energy 
range  (thermal  through  14  MeV)  and  fluence  rang* 
(3  * 10"  to  3 x 10‘*  Hem1),  freedom  from  gamma - 
ray  interference,  good  reproducibility,  ease  of 
access  to  remote  regions  in  a reactor,  and  rela- 
tively good  economy  It  will  not  work  at  distances 
greater  than  -30  cm  from  a reactor10  in  a low-Z 
shield  such  as  LiH  or  H*0,  nor  can  it  be  used  with 
a Linac  source  (in  the  forward  direction)  due  to 
the  significant  number  of  (y, fission)  reactions. 
With  a Linac -driven  source,  TOT  methods  are 
much  more  accurate."  However,  for  Lmac-pulsed 
reactor  assemblies,  a long  flight  puth,  complex 
dteaway-time  corrections,  and  a low  J (-0.9  or 
less)  are  required.’"*’ 10 


TO.  INTERACTIVE  OPERATION  OF  3AND-IJ 

The  method  of  threshold -foil  spectrometry  can 
be  described  as  follows.  A set  of  foils  is  selected 
with  threshold  energies,  £,.  that  vary  from  0 to 
14  MeV  (see  Table  1).  The  foils  are  weighed, 
covered  with  cadmium  or  l0B  as  necessary,  and 
exposed  to  a neutron  field  having  a spectral 
intensity  ««(£).  The  resultant  radioactivity  is 
measured  with  a high -efficiency  Ge(Li)  or  in- 
trinsic germanium  gamma- ray  detector  of  good 
resolution  to  determine  the  number  ,V  of  radioac- 
tive atoms  produced  during  the  irradiation  of  ,V0 
target  atoms  of  the  specific  isotope.  For  each  foil 
species  x, 

S,  . ,V0  /“  <7,(£)eo(£)<f£  (1) 


**V  V VERBINSKI,  T FUSE.  1.  0 KINGTON,  ind  K.  M. 
HENRY.  "Transport  of  Ftst  Neutrons  Through  Lithium  Hy- 
dride Part  II  Angular  Distribution*,’*  ORNLJ360,  Neutron 
Physics  Division  Annusi  Progress  Report  for  Penod  Ending 
September  I,  1441,  Oak  Ridge  Niuonai  Laboratory  (1961). 

"V  V VERBINSKI.  M S.  BOKHAR1,  J C.  COURTNEY, 
and  G E.  WHITESIDES.  .Vucl  Set.  Eng..  *7.  283  (1967). 

UJ.  M.  YOUNC.  I M NEILL  P.  d’OULTREMONT.  E.  L 
SLAGC1E.  and  C.  a.  PRESKITT.  \ucL  So.  £nf..  48.  1*1 
(1971). 

'»?  d’OU’.TREMONT,  D H.  HOUSTON,  and  J.  C YOUNG. 
"CAGE-BIRO-SPEC,"  Gulf.RT-IOI9J,  Cutf  Radiation  Tech- 
noiogy  (1970). 


where  j,i£)  is  the  cross  section  for  the  reaction. 
The  set  of  measured  specific  activities  given  by 
the  rstio  KmUc)  » l.V  .V„),  is  used  as  input  to  the 
SAND- Q unfolding  code,  along  with  a trial  spec- 
trum a„l£)  that  contains  all  the  physical  features 
(see  Sec.  IV). 

Using  Eq.  (1)  with  ou(£)  in  place  of  the  un- 
known, «*(£),  the  SAND-0  code  calculates  the 
various  ratios  R«(x)  that  are  consistent  with  o„(£) 
by  meins  of  the  SAND-0  Input  library"  of  cross 
sections,  o,(£). 

The  frictional  differences. 


Ao(x)  • 


Rmix)  * Rft (x) 
RJx) 


(2) 


for  all  the  foils  are  calculated  by  the  SAND-0 
code,  and  the  standard  deviation,  S«,  is  generated 
from  the  set  of  Aoix).  If  So  is  <0.05  (or  any  other 
value  the  user  may  input),  then  o „(E)  is  the  solu- 
tion. If  not,  the  code  raises  (lowers)  the  spectrum 
0vi£)  In  the  regions  just  above  the  threshold 
energy  £,(x)  for  which  the  corresponding  Ao(x)  is 
large  and  positive  (negative).  The  process  is 
repeated  to  generate  the  sets  ®i(x),  . . . , ®„lx), 

At(x).  ....  A^x),  and  $ SL  where  S«  s 

0.03.  If  * * 10  or  less,  the  probability  that  o„  has 
little  spurious  structure  is  quite  good. 

If  is  is  large  and/or  if  considerable  spurious 
structure  Is  observed  [as  determined  most  sensi- 
tively by  comparing  E o(E)  plots  of  o„(£)  and 
o„l£)],  then  the  user  must  examine  the  Ao(x)  set 
for  a value  (or  values)  that  is  appreciably  differ- 
ent from  the  remainder  of  the  set,  eliminate  the 
questionable  value  as  not  being  consistent  with  the 
remainder  of  the  set,  and  rerun  SAND-Q  without 
the  corresponding  R„(x).  Alternatively,  the  user 
can  also  restructure  o„(£)  to  yield  a shape  that 
results  in  a better  set  of  Ao(x)  and  that  therefore 
agrees  with  the  measured  activation  set,  A*(x). 


rv.  A PRIORI  CONSTRAINTS 

The  limited  information  provided  by  the  dozen 
or  ao  measured  foil-activity  values,  E«(x),  means 
that  the  threshold -foil  spectrometry  method  re- 
quires the  a priori  information  of  a meaningful 
trial  spectrum,  ««(£).  A simple  one-dimensional 
mockup  of  the  reactor  and  surrounding  scattering 
material  Is  usually  adequate  to  provide  a calcu- 
lated spectrum,  Oc.  that  contains  all  the  physics  of 
the  problem.  Alternatively,  the  parameterized 
family  of  curves  presented  in  Sec.  VU1  can  be 
used  for  o„i£)  measurements  of  an  unshielded 
rsactor  core. 


“For  an  sarly  summary.  s«c  R.  L.  SIMONS  and  W N 
McELROY,  "EvaJuattd  Reference  Cross-Section  Lbrary.” 
BNWL-1312.  Battalia-Northwest  Laboratories  ( 1970). 
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TABLE  l 

Activation  Foil*  (1.27 -cm  dlio) 


in.  , im  . 
Aula,  y)  Au 

’•co<«.Yrco 


"Np(i../)l*U 

"Mpta./l'Zr 


‘“UKs.s')1*^ 

1.0 

a*U  (s,/)1—La 

1.46 

“Th(«./)l,*Bs 

1.73 

“Th(«./),,Zr 

1.75 

MFs(s,p)MMn 
**Nl  (*,  p)‘*Co 
"S(»,P)”P 

2.20 

“Mgle.PI^Ns 

6.3 

"Fe<»,  P)**Mn 

7.5 

WA1  (».  o)MNs 

8.7 

"Ke.ae)1*! 

11.0 

aZr(e.2ii)“Zr 

14 

Gamma/  Reaction 
(Fission  Yield.  %) 

Ti* 

Foil  Mass 
(O 

0.96 

2.696  days 

0.066 

1.00 

3.288  yr 

0.057 

1.00 

0.99 

2.58  h 

0.030 

0.29 

0.96  (6.29) 

40.23  h 

0.281 

0.96  (3.24) 

40.23  h 

0.150 

0.96  (3.69) 

40.23  h 

0.380 

0.92  (6.01) 

16.8  h 

0.380 

0.30 

4.50  h 

0.235 

0.96  (6.02) 

40.23  h 

0.495 

0.256-  (7.91) 

12.8  days 

1.066 

0.92  (4.12) 

16.8  h 

1.066 

1.00 

303  days 

0.142 

0.99 

71.3  days 

0.282 

1.000 

14.3  days 

4.10 

1.00 

13.0  h 

0.030 

0.99 

2.38  h 

0.142 

0.29 

1.00 

15.0  h 

0.257 

0.34 

13.06  days 

0.687 

0.33 

0.99 

78.4  h 

0.108 

Isotopic  j 
Abundance  i 


Cadmium  covar  0.06  to  0.10  cm  thick. 

bVse  **Co  instead  o(  *Au  and  **Mn  for  vary  loaf  irradiations. 

eSt  * 0.01  MeV  with  mB  sphere  [Important  for  aoft  (TRIG A)  spactra,  whara  P(S)  < 0.01  will  otherwise  domiaatel. 
Whan  “U  or  °*Pu  foil  la  03-cavered,  also  covar  a(J  and  ^Np  foils  so  that  accurata  corractlons  can  bo  mads  (or 
*U  and  ”*Pu  lmpurttlas  la  thasa  high  E,  foils. 

°40.23-h  daughter  of  12.80-day  l4%a.  Walt  flvo  days  for  maximum  decay  rata  (ASTM  E 393). 

'Fission  yields  are  for  bombardment  with  fission-spectrum  neutrons.  For  thermal  and  14  MeV,  saa  Ref.  S. 

'Use  "Zr  for  low  duaooa  (3  * 10U  to  3 x 10un/om’).  Use  peak  shape  analysis  or  measure  twice,  seven  days  apart, 
to  strip  oS  740-kaV  aMo  gamma  ray  (Tl/t  - 67  h). 

*alTh  radioactivity  interferes  with  l*La  line. 

‘Requires  separata  detector,  and  calibration  technique  Is  complex. 

‘Maximum  manganese  impurity  - 0.001%,  cadmium -covered.  Omit  **Fe  la.pl  **Mn  for  long  irradiations. 


Not*  from  Table  I that  a "fabricated”  thresh- 
old appears  at  ~0.0l  MeV  as  produced  by  a thick 


B cower  for  the 


Pu  fission  foil.  The 


only  other  useful  threshold  up  to  l MeV  Is  that  for 
the  mNp<  n,/)  reaction,  where  the  threshold  energy 
E,  * 0.3  MeV.  (Data  below  0.01  MeV  are  usually 
of  little  use  la  radiation  damage  work.)  Conse- 


quently, the  output  of  SAND- Q or  any  other  code 
exhibits  an  extremely  large  "confidence  interval" 
or  "region  of  likely  solutions"  below  l MeV. 
unless  It  la  constrained  to  act  as  a small -pertur- 
bation calculation  on  a good  trial  spectrum  to 
provide  a final  spectrum  consistent  with  the  foil 
activations. 
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The  second  constraint  is  that  which  is  built  Into 
SAND-U.  This  is  an  improved  version  of  the 
original  iterative  unfolding  code1'  which  was  re- 
programmed to  retard  the  formation  of  spurious 
structure1'  while  preserving  the  physically  mean- 
ingful structure,  i.e.,  a basic  fission  spectrum 
with  varying  strength  of  1/E  slowing  down  compo- 
nent attached  at  the  low-energy  side,  a resonance 
''dip,"  a monoenergetic  peak  at  ~14  MeV,  or  a 
broad  thermal -neutron  Maxwellian  peak. 

The  third  constraint  is  a user-interaction  con- 
straint of  keeping  the  number  of  iterations,  n, 
down  to  a low  value  (see  Sec.  III).  This  is  accom- 
plished by  eliminating  spurious  (oil  activations 
and/or  modifying  PB(E)  in  regions  above  a (oil 

threshold  energy  £,*  for  the  foil  yielding  the 

largest  disparity  with  the  first  a „(E)  used.  These 

disparities  [boUt)  in  Sec.  Ill)  are  printed  out 

in  SAND-U  for  subsequent  interactions  as  well 
[b,Cr)j.  Also,  the  energy  "band"  where  5 and 
95%  of  the  activation  has  occurred  for  each 
foil  x is  listed  for  each  spectrum  perturbation 
(bo(£)  * $„(£),  <p»(£),  ....  «.(£))  as  an  additional 
guide. 

The  final  a prion  constraint  is  that  of  accurate 
nuclear  data,  i.e.,  an  evaluated  consistent  set  of 
threshold-foil  cross  sections1*  input  to  the  code, 
and  of  accurate  half-lives  and  fluorescent  yields 
(conversion  coefficients,  branching  ratios,  and 
fission  yields)  available  for  processing  the  foil- 
activity  data. 

V.  CALCULATION  OF  DIE)  FOR  SILICON 

The  damage  function  for  silicon.17  D(E),  used  in 
the  present  work  is  shown  in  Fig.  I (dotted  line) 
and  in  Table  II  tlO"'  to  2 * 10’  eV).  It  was  ob- 
tained with  neutron  cross  sections  for  silicon 
from  the  ENDF  B-IV  MATU94  evaluation.1'  De- 
tails of  the  calculational  method  are  described  in 
Ref.  17.  The  results  of  a previous  calculation,  l* 
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Fig.  I.  Silicon  kemi  ior  displacement  and  lifetime  degra- 
dation calculated  for  neutron  energies  from  0.01  to  10  MeV 


also  shown  in  Fig.  I,  were  based  on  neutron  cross 
sections  available  in  1965  and  1966. 

Although  the  uncertainty  of  the  calculation  of 
total  energy  [ionisation,  HE),  and  displacement, 
D(E) j is  mainly  governed  by  the  neutron  cross 
sections,  the  partition  of  total  energy  deposition 
between  /(£)  and  D(E)  is  also  limited  by  the 
accuracy  of  the  Lindhard  factor  for  determining 
the  fraction  of  recoil-atom  and  charged -particle 
kinetic  energy  that  is  available  to  cause  displace- 
ments, and  by  the  Robinson14  functional  rep- 
resentation of  the  Lindhard  factor.  The  Robinson 
representation  needs  to  be  improved  by  fitting  to 
the  very  few  data  points  presently  available  for 
energy  partition  and  by  future  careful  measure- 
ments. It  is  hoped  that  present  D(E)  data  will  be 
updated  when  the  work  on  the  energy  partition  Is 
accomplished. 


'•W  N McELROY  et  al.,  “A  Computer- Automated  Iterative 
Method  for  Neutron  Flu*  Spectra  Determination  by  Foil  Acti- 
vation.'' .aFWL-TR-oT-M.  Vol.  I.  .Air  Force  Weapon*  Laboratory 
liven. 

'*5.  BERC  and  W N McELROY.  "A  Computer- Automated 
Iterative  Method  for  Neutron  Flu*  Spectra  Determination  by 
Foil  Activation.''  AFWL  TRolAl,  Vol.  11,  .Aar  Force  Weapons 
Laboratory  1 1 967) 

”V  c.  ROGERS.  L HARRIS.  lr..  D K.  STEINMAN,  and 
D E.  BRYAN.  IEEE  rnns.  ,Vu cl  Set..  NS- 23.  5,  1316  (19751; 
trrtntm.  NS-BS.  /.  S75  11976). 

“D  C LARSON  and  F G.  PEREY.  "ENDF  B Evaluation  of 
Neutron  and  Gamma- Ray  Production  Croia  Sections  for  Sili- 
con.'' published  ut  ENDF/B  Summary  Documentation,  BNL- 
17541  t ENDF- :on.  National  Neutron  Cross  Section  Center. 
Brookluven  National  Laboratory  1 1973.  plus  updates) 

‘*R-  R.  HOLMES,  unpublished  data  11970). 


VI.  THE  l- MeV  EQUIVALENT  FLUX 
CHARACTERIZATION  OF  A 
NEUTRON  HELD 

A neutron  field  can  be  characterized  by 
(1  MeV),  which  represents  the  fluence  of  1-MeV 
neutrons  required  to  produce  the  same  radiation 
damage  in  bulk  silicon.  Since  nearly  all  the 
neutron  damage  is  produced  by  the  neutrons  above 
0.01  MeV.  even  for  a soft  TRIG  A spectrum  (see 
below),  we  define  as 


*M.  T ROBINSON,  Proc.  Sucltar  Fusion  Rtactors  Conf.. 
p.  364.  British  Nuclear  Energy  Society.  Culham  Laboratory 
(1969);  see  also.  0.  G.  DORAN.  Vu cL  Sci  Enf..  49.  130  (1971). 
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TABLE  II 


Calculated  Silicon  Displacement  and  Ionization  Ktrni 


Energy 

<eV) 

Silicon  Kerma 
(MeV-rnb) 

Energy 

(eV) 

Silicon  Kerma 
IMeV-mb) 

Energy 

(eV) 

Silicon  Kerma 
(MeV -mb) 

Displacement 

lomxanon 

Displacement 

lonixatioQ 

Displacement 

Ionization 

1.000-3* 

8.991-0 

2.370-0 

1.242-3 

1.328-1 

2.728-2 

3.162*5 

3.190*1 

2.831*1 

1.S57-S 

6.604-0 

1.886-0 

1.885*3 

2.288-1 

4.378-2 

3.548*5 

3.132*1 

2.396*1 

3.449-3 

4.830-0 

1.386-0 

2.861*3 

3.426-1 

7.038-2 

3.981*5 

3.260-1 

3.073*1 

V 405-5 

3.362-0 

1.018-0 

4.342*3 

3.131-1 

1.132-1 

4.467*5 

5.441*1 

3.314*1 

1.159-4 

2.613-0 

7.473-1 

6.389-3 

7.604-1 

1.806-1 

5.012*5 

5.715*1 

3.369*1 

2.209-4 

1.917-0 

3.480-1 

1.000-4 

1.041*0 

2.663-1 

3.623*5 

1.332*2 

1.240*2 

4.102-4 

1.406-0 

4.018-1 

1.112*4 

1.122*0 

2.932-1 

6.310*5 

5.409*1 

3.444*1 

7.617-4 

1.031-0 

2.946-1 

1.239-4 

1.203-0 

3.208-1 

7.080*5 

6.056*1 

4.290*1 

1.429-3 

7.363-1 

2.162-1 

1.413*4 

1.279*0 

3.466-1 

7.943*5 

9.354*1 

8.767*1 

2.627-3 

3.334-1 

1.388-1 

1.385*4 

1.368*0 

3.308-1 

8.913*5 

7.309*1 

5.273*1 

4.879-3 

4.080-1 

1.168-1 

1.778*4 

1.478*0 

4.202-1 

1.000*6 

1.017*2 

8.989*1 

9.060-3 

2.996-1 

3.363-2 

1.993*4 

1.387-0 

4.812-1 

1.162*6 

4.438*1 

3.182*1 

1.663-2 

2.192-1 

6.267-2 

2.239*4 

1.727*0 

3.131-1 

1.349*6 

8.994*1 

7.420*1 

3.123-2 

1.614-1 

4.613-2 

2.312*4 

1.874*0 

3.691-1 

1.567*6 

9.936*1 

9.353*1 

3.802-2 

1.178-1 

3.368-2 

2.818-4 

2.037*0 

6.388-1 

1.821*6 

9.654*1 

9.334*1  1 

1.078-1 

8.678-2 

2.480-2 

3.162*4 

2.261*0 

7.183-1 

2.113*6 

1.348*2 

1.309*2 

2.001-1 

8.374-2 

1.821-2 

3.348*4 

2.301*0 

8.128-1 

2.437*6 

1.593*2 

1.879*2 

3.718-1 

4.670-2 

1.334-2 

3.981*4 

2.323*0 

9.391-1 

2.833*6 

2.003*2 

2.283*2 

6.901-1 

3.444-2 

9.824-3 

4.467  .4 

3.309*0 

1.127*0 

3.315*8 

1.187*2 

1.457.2 

1.282-0 

2.329-2 

7.191-3 

3.012*4 

3.897*0 

1.360-0 

3.830*8 

1.212*2 

1.653*2 

2.380-0 

1.382-2 

5.308-3 

3.623*4 

4.738*0 

1.694*0 

4.472*6 

1.795*2 

2.615*2  1 

3.612-0 

1.343-2 

4.303-3 

6.310*4 

6.304*0 

2.384*0 

5.193*6 

1.932*2 

3.127*2  | 

5.482-0 

1.283-2 

3.517-3 

7.080*4 

7.247*0 

2.723*0 

6.034*6 

1.965*2 

5.236*2 

6.319-0 

1.100-2 

2.911-3 

7.943*4 

6.481*0 

2.301*0 

7.009*6 

2.014-2 

1.061*3 

1.263-1 

9.727-3 

2.439-3 

8.913*4 

3.735*0 

2.272*0 

3.142*6 

1.734-2 

1.680*3 

1.916-1 

9.036-3 

2.088-3 

1.000*5 

3,287.0 

2.152*0 

9.437*6 

2.006*2 

2.272*3 

2.9C8-1 

9.127-3 

1.878-3 

1.122*5 

3.024*0 

2.104-0 

1.099*7 

1.981*2 

2.386*3 

4.413-1 

1.006-3 

1.812-3 

1.239*5 

4.008*0 

1.728*0 

1.276-7 

2.134*2 

2.345*3 

6.898-1 

1.220-2 

1.947-3 

1.413*5 

2.643*0 

1.173-0 

1.482*7 

2.134-2 

3.086*3 

1.017.2 

1.393-2 

2.343-3 

1.585*5 

5.231*0 

2.398*0 

1.722*7 

2.213*2 

3.853*3 

1.543-2 

2.205-2 

3.131-3 

1.778*3 

3.220*1 

2.461*1 

2.000*7 

2.297*2 

3.073*3 

2.341-2 

3.162-2 

4.532-3 

1.995*3 

1.086*2 

5.253*1 

3.333*2 

4.628-2 

6.883-3 

2.239*3 

7.566*1 

3.748*1 

5.393-2 

6.833-2 

1.078-2 

2.312*5 

6.041*1 

3.089-1 

8.164*2 

1.022-1 

1.707-2 

2.818*5 

3.423*1 

2.863-1 

'Rud  as  1.000  * 10*'. 


ia  m«v 

/ *{E)D(E)dE 

“ "“dU  M.V)  <3> 

Because  of  the  rapid  variation  of  D(E)  near 
1 MeV,  we  define  D (1  MeV)  as  the  average  value 
between  0.83  and  1.15  MeV 

It  is  also  useful  to  define  a neutron- spectrum- 
shape  parameter  o^o.  which  is  the  number  of 
1-MeV  neutrons  per  unit  fluence  of  the  neutron 
field  being  specified  required  to  product  the  same 
damage,  t.e.. 


li  MtV 

/ o(  £)£>(£)<*£ 


0(1  MeV)  / o (£)d£ 

Ml 

By  averaging  £X£)  over  the  energy  interval  of  0.35 
to  1.15  MeV,  o is  slightly  larger  than  unity  for 
a GODIVA  spectrum.  This  is  intuitively  prefer- 
able to  a value  less  than  unity  [which  would  result 
from  using  D (1  MeV)  instead  of  D (0.35  to  1.15 
MeV)),  considering  that  the  average  energy  ,s 
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>1  MeV  and  that  D(E)  increases  with  energy.  By 
excluding  the  neutrons  below  0.01  MeV  (which  con- 
tribute only  ~1'$  to  0,,/e  even  (or  the  soft  TRIGA 
spectrum),  the  value  of  o.,,  o becomes  independent 
of  unimportant  changes  in  neutron  flux  at  sub- 
cadmium and  sub-‘°B  energies. 

Vn.  MEASUREMENTS  OF  THREE 
STANDARD  SPECTRA 

VII  ji.  Description  of  Measurements  and 
Method  of  Accuracy  Assessment 

Measurements  were  made  of  three  different 
reactor  spectra,  the  glory-hole  and  leakage  (50 
cm  from  core)  spectra  of  a bare  GODIVA-type 
fast  burst  reactor  (FBR)  at  White  Sands  Missile 
Range  (WSMR)  and  the  leakage  spectrum  from  the 
General  Atomic  TRIGA  reactor.  The  first  two 
were  measured  by  WSMR  personnel1'  m 1974  and 
by  IRT  personnel  in  1975  to  obtain  a measure  of 
the  precision  of  the  threshold-foil  spectrometry 
method.  The  two  groups  used  different  foil  sets, 
Ge(Li)  gamma-ray  spectrometers,  and  National 
Bureau  of  Standards  (NBS)  calibrated  gamma-ray 
reference  standards. 


Each  spectrum  was  unfolded  by  using  the 
corresponding  “standard"  spectrum  (see  the  Ap- 
pendix) as  a trial  spectrum,  o „(£),  in  the  SAND-U 
unfolding  code.  Agreement  between  the  SAND-Q 
result,  o .(E),  and  Oi,  m this  case  gives  a mea- 
sure of  the  accuracy  of  the  foil  activation  mea- 
surements, the  published  branching  ratios  or 
fluorescent  yields  given  in  Table  I,  and  the 
threshold-foil  crosa-section  library  provided  with 
the  SAND-il  code.”  The  sensitivity  of  the  thresh- 
old-foil method  to  ou,  to  the  type  of  unfolding  code 
used,  and  to  the  individual  foil  cross  sections  (or 
foil  activations)  is  discussed  in  Secs.  VH1.  DC.  and 
X. 

va  B.  Results 

Table  m lists  the  foil  activations  for  the  five 
measurements  described  above.  The  FBR  foil 
activations  measured  by  IRT  were  measured  by 
WSMR  personnel  as  well,  using  the  same  foil  set 
but  different  detectors  and  NBS  gamma- ray  cali- 
bration standards.  The  two  agreed  with  an  aver- 
age standard  deviation  of  ~4.5x>,  although  the  wait 
(decay)  times  were  appreciably  different  in  most 


"J  MEASON.  H.  WRIGHT.  I HOCAN.  and  J.  HARVEY. 
IEEE  Trans.  \ucL  Set.  NS- S3.  4.  1330  (1975) 


USAND-I1  u available  tt  'he  Radiation  Shielding  Intorma- 
uon  Center,  Oak  Ridie  National  Laboratory,  Oak  Ridge,  Ten- 
nessee. 


TABLE  IB 


Foil  Activations  for  Five  Measurements:  Fission  Foils  • Fission/  Nucleus:  Others  • 
Disintegration/ (s-nucleus)  as  Input  to  SAND-U  In  Time -Integrated 
Mode  (l.e. . corrected  to  Instantaneous  Irradiation) 


FBR  (IRT) 

FBR  (WSMR) 

TRIGA  (CRT) 

60  cm 

Glory  Holo 

50  cm 

Glory  Hole 

./-Tube 

'".Aui-i.y)* 

1.79-17* 

8.05-17* 

2.06-17* 

4.55-17* 

‘"Au(e.y) 

2.83-17 

8.00-17 

3.78-17 

1.13-18 

“u  (»./■) 

9.94-12* 

2.34-10 

1.24-11* 

1.45-09* 

5.91-12* 

a,Pu(e./) 

1.71-11* 

1.82-09* 

"><-.,/) 

6.42-12* 

1.49-10 

9.56-12 

1.03-09* 

3.37-13* 

‘“ln(e,»') 

3.92-17 

8.71-16 

6.08-17 

2.01-18 

*"u  <»./) 

1.29-12* 

3.34-11 

1.97-12* 

1.88-10* 

8.67-14* 

a*Th(n./) 

3.81-13 

8.05-12 

4.85-13 

3.93-11 

1.83-14* 

**Fe  (n.P) 

1.08-20 

2.58-19 

1.54-20 

1.37-18 

**Nl  (*,p> 

5.77-20 

1.49-18 

9.18-20 

8.07-18 

3.01-21 

"SO..*) 

t. 72-19 

2.57-19 

8.97-21' 

Mg(e.p) 

9.42-20 

2.23-18 

1.33-19 

1.25-17 

7.31-21 

**Fe(e,  p) 

4.90-19 

9.24-18 

5.95-19 

5.38-17 

^Ai  (Ha  a) 

4.57-20 

1.06-18 

6.47-20 

3.38-21 

"Ke.lel 

2.58-21 

4.79-21 

4.60-19 

2.28-22 

*Zr(n,2n) 

1.73-19 

‘Cadmium -covered  (~d.050  cm);  read  as  1.79  * 10'”. 
*‘°B-covered  (-1.45  g/cmV 
c Deduced.  SAND-U  run. 
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c ases.  This  makes  the  companion  a "worst 
case"  situation  in  that  sense. 

Figure  2 shows  the  results  of  the  measured 
leakage  spectrum  at  SO  cm  from  the  center  of  the 
WSMR  FBR  (solid  circles!  using  the  197S  1RT  foil 
activation  measurements  as  input  to  the  SAND-Q 
unfolding  code.  The  spectrum  calculated  with  the 
IDF  code  (see  the  Appendix)  was  used  as  the 
trial  spectrum  (solid-line  histogram  Fig.  2).  The 
agreement  in  shape  between  the  SAND -a  result  and 
the  IDF  calculation  is  quite  good,  indicating  that 
the  gamma -ray  activation  measurements,  the  gam- 
ma-ray and  fission -fragment  yields  of  Table  I, 
and  the  SAND-Q  foil -activation  cross-section  li- 
brary are  reasonably  accurate.  The  SAND-Q 
spectrum  was  arbitrarily  normalized  to  the  60-cm 
calculation  at  0.3  to  1.3  MeV  for  shape  com- 
parison. 

Table  IV  presents  the  values  of  g^'o  for  the 
WSMR  measurements  (spectrum  A)  as  well  as  the 
IRT  results  (spectrum  B).  The  two  values  of 
o^/o  (1.080  and  1.060)  bracket  the  IDF  value  of 
1.077,  which  shows  good  agreement  in  spectral 
shape  for  all  three  cases.  The  flux  above  0.01 
MeV  per  sulfur  activation  varied  from  the 
mean  of  3.59  x 10*"  "P  nuclei  per  "S  nucleus  per 
watt-second.  This  gives  a measure  of  the  ac- 
curacy of  specifying  the  silicon- semiconductor 
radiation  damage  once  the  spectral  shape  is 
accurately  known.  If  it  is  not,  the  sulfur  pellet  is 
a poor  choice  for  monitoring  radiation  damage 
(see  Sec.  X1I1.C). 

The  IDF  calculations!  results  at  14  and  400  cm 
from  the  FBR  center  are  also  shown  in  Fig.  2 


NEUTRON  ENERGY  (MeV) 

Fla  2.  Companion  of  threshold-foil  spectrum  t SAND-11 
tolid  circles)  with  IDF  calculation  (solid-line  histogram)  at 
60  cm  from  tha  WSMR  FBR  canter.  The  IDF  multi  are  alio 
shown  at  14  and  400  cm  to  illuatme  the  importance  of  \,E 
moderation  component  (flat  line)  due  to  wall  scattering. 


(both  arbitrarily  normalized  to  the  60-cm  spec- 
trum at  0.5  to  1.3  MeV  for  shape  comparison)  to 
illustrate  the  increasing  importance  of  the  mod- 
erated component  due  to  wall  and  floor  scattering 
as  the  distance  from  the  reactor  increases. 

The  measurement  of  the  WSMR  FBR  glory-hole 
spectrum  is  shown  in  Fig.  3 along  with  the  IDF 
calculation  used  as  os  in  the  SAND-Q  unfolding 
code.  The  1974  (WSMR)  and  1975  (IRT)  results 
were  practically  identical  (Table  IV.  spectra  C 
and  D)  yielding  o^/o  • 1.037  and  1.023,  respec- 
tively. These  values  again  bracket  the  IDF  result 
of  o,q/0  * 1.030.  The  flux  above  0.01  MeV  per 
sulfur  activation  agreed  within  2*$  for  the  two 
measurements  in  this  case. 

A measurement  for  the  TRIGA  swimming-pool 
reactor  was  made  in  the  remote-access  ■/-tube11 
that  abuts  the  reactor  core.  A shield  of  7.6-cm 


Fig.  3.  The  SAND-ll  result!  for  the  WSMR  FBR  ^lory-hole 
spectrum  with  »„t£)  * IDF  calculation  Table  IV  C and  Cl 


“G.  B.  WEST.  "Calculated  Fluxes  and  Cross  Sections  for 
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TABLE  IV 

Data  for  Individual  Spectra 


Unfolded  with  SPECTRA  code  (lladlallon  Shielding  Information  Center.  Oak  nidge). 
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lead  and  0.6-cm  Boral  (thermal-neutron  trans- 
mission »5.6  x 10'4)  separates  the  core  from  the 
7- lube  cavity  used  (or  radiation  damage  studies. 
The  measured  spectrum  is  shown  m Fig.  4 along 
with  the  best  available  calculated  spectrum  that 
was  used  as  oB(£)  m the  SAND-D  unfolding  code. 
In  this  case,  •«(£)  is  a composite  of  the  NIOBE 
code  calculation,  which  had  a low-energy  cut-off 
of  0.8  MeV,  and  the  GAZE  calculation”  (see 
Appendix).  This  measurement  (spectrum  £ of 
Table  IV)  yielded  a value  of  e^/  e that  was  within 
'3.4%  of  the  calculations!  value.  The  SAND-n 
result  using  the  GAZE  code  as  o»(£)  was  S% 
higher  in  kyo  than  the  corresponding  GAZE 
value  (spectrum  F of  Table  IV).  Thus,  the 
measurements  are  In  best  agreement  with  the 
composite  calculated  spectrum,  which  was  experi- 
mentally verified  above  0.8  MeV  with  an  indepen- 
dent spectrometry  method. 

VIII.  PARAMETERIZING  «„<£> 

FOR  REACTOR  SPECTRA 

A GO  DIVA  spectrum  (spectrum  5 of  the  SAND- 
Q library  of  trial  spectra)  is  proposed  as  the 
basic  reactor  spectrum,  but  with  a varying 


strength  of  l.  £ moderating  component  that  de- 
pends on  the  amount  of  moderator  present.  This 
GO  DIVA  spectrum  without  any  moderating  compo- 
nent (£iau  • 0)  is  adequate  as  o„l£)  for  the  glory- 
hole  spectrum  (Table  IV  H).  A comparison  with 
the  SAND-11  result  with  o»(£)  ■ IDF  calculation 
shows  agreement  within  2%  in  0**/*  and  1%  in 
0 > 0.01  MeV. 

At  0.5  m from  the  reactor,  neutron  moderation 
from  the  concrete  floor  '2  m below  yielded  a 1.  £ 
component  that  can  be  parameterized  by  flMa  ■ 
0.01  MeV.  Figure  3 shows  the  results  of  three 
SAND-n  runs  with  o„l£)  • GO  DIVA  * 1 £ modera- 
tion component  Joined  at  £ » 0.003,  0.01,  and  0 02 
MeV  (Table  IV  I.  J,  and  K).  The  best  fit  of  the 
SAND-n  output  to  o„(£)  occurs  for  £,««  ■ 0.01 
MeV.  where  o,^  o » 1.080  for  the  SAND-n  output 
versus  1.079  for  o*l£)  with  £,«.  * 0.01  MeV 
(Table  IV  J).  This  compares  to  1.060  for  SAND-U 
with  o*(£)  « IDF,  and  a^/ o « 1.077  for  the  IDF 
spectrum.  In  addition,  values  of  o above  10  keV 
agree  within  2%. 

For  the  TRIGA  7-tube  spectrum,  where  the 
NIOBE  calculation  above  0.8  MeV  agrees  well  with 
the  GODIVA  spectrum,  we  find  £|OI,  » 0.2  MeV 
(Table  IV  G).  Here,  agreement  with  the  baseline 


Fig.  4.  The  TRIGA  /-lube  spectrum  SAND-11  results  with  e«Af>  “ NIOBE  (>0.8  MeV>  * CODIVa  (0.3  (o  0.8  MeV)  * |,£; 
three  gold  foils  (cadmium-covered).  Corrections  for  foil  thickness  (0.0025  cm)  ue  shown  tTable  (V  G! 
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Fig.  5.  Three  unfolded  result!  of  the  WSMR  FBR  leakage  diu  with  Pn(£1  * "GODIVa"  (Stan  number  J of  SAND- 1 1 library 
of  tnal  spectral  and  I/£  tail  (flat  line  in  above  plot)  joining  at  £,oul  ■ 0.005,  0.01.  and  0.02  MeV  Note  good  agreement  for  £,„.«• 
0.0 1 MeV 


case  (Fig.  4 and  Table  IV. E)  is  1$  in  and 
in  o. 

In  summary,  £ » 0,  0.01,  and  0.2  MeV  for  the 
FBR  glory  hole,  FBR  0.5-m  and  TRIGA  ./-tube 
spectra,  respectively,  where  the  degree  of  moder- 
ation becomes  progressively  greater. 

IX.  PROBLEMS  WITH  LOW-ENERGY  NEUTRONS 

The  three  low-energy  curves  of  Fig.  4 
represent  SAND-Q  runs  with  the  three  differ- 
ent gold-foil-resonance  self -absorption  correc- 
tions shown  (1.0,  2.0,  and  2.5).  Note  that  the 
gold-foil  self-absorption  correction  of  2.0  for  the 
4.95-eV  resonance  gives  best  agreement  with  a 
1 /£  spectrum  (Fig.  4).  This  is  for  nearly  for- 
ward-directed  neutrons  m the  ./-tube  and  is  in 
qualitative  agreement  with  the  factor  of  2.5  of 
Sanford  and  Seckwger1*  for  isotropic  flux  0.0025- 
cm-thick  foil.  The  isotropic  flux  correction  is,  of 
course,  always  the  larger  of  the  two. 


!*G-  S.  SANFORD  ind  J.  H.  SECXINGER.  "Thickne!! 
Correction!  for  Neutron-Activated  Gold  Foili."  AN L'545, 
Argonne  NanonaJ  laboratory  1 19691. 
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The  ratio  of  bare- to-cadmium -covered  gold 
foils  (Table  111)  shows  that  the  TRIGA  ./-tube 
thermal  flux  is  appreciable.  Although  the  TRIGA 
J- tube  is  Boral-Uned,  the  Boral  leakage  plus  the 
scattering  materials  inside  the  ./-tube  add  a 
noticeable  amount  of  thermal  component.  This 
stresses  the  importance  of  covering  all  thermal- 
activation  foils  with  cadmium  if  a 10B  covering  is 
not  used.  This  holds  especially  true  for  foils 
where  activations  of  foil  impurities  by  thermal 
neutrons  compete  with  finite-threshold  activations. 
i.e.,  ,JNa  in  "Al(>«.a)1,Na,  **Mg  in  “Fetn1p)“Mn. 
,J*U  in  ’"Udt./).  J*Pu  in  ,5TNp(*/),  etc.  For  the 
soft  TRIGA  spectrum,  this  effect  is  so  serious  as 
to  rule  out  the  use  of  most  of  these  reactions 
unless  the  impurities  are  in  the  few  parts  per 
million  range,  and  impurity  foils  are  exposed  with 
the  same  cover  to  provide  reliable  impurity 
correction  data. 

X.  SENSITIVITY  TO  a^El 

Two  SAND-n  unfolding  results  were  obtained 
that  illustrate  that  the  solution  does  not.  m gener- 
al, converge  to  a physically  meaningful  spectral 
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shape  unless  o»(£)  brings  the  a pnon  constraint 
of  a physically  meaningful  spectrum  into  the 
calculation. 

Figure  6 shows  the  results  of  a SAND- II  un- 
folding operation  for  the  WSMR  FBR  leakage 
spectrum  at  SO  cm  from  the  reactor.  Using 
o«(£)  • GODtVA  (solid  line,  Fig.  8)  and  omitting 
the  low-energy  gold  foil  data  (cadmium-covered), 
the  results  obtained  (represented  by  solid  circles) 
agree  with  ®,(£)  ■ GODIVA  (Table  IV. L)  because 
the  calculated  *"U  and  M*Pu  activations  (5  to  95% 
~nnts)  are  produced  in  the  0.18-  to  4.3-MeV  and 
6.2-  to  4.4-MeV  bands,  respectively.  If  the  gold 
foil  data  are  added,  the  spectrum  designated  by 
the  symbol  "X”  results  (Table  IV .M).  This  is  a 
clearly  unreal  spectrum.  Here,  o^  o has  dropped 
from  1.090  (Table  IV  A)  to  0.818,  the  total  flux  is 
8%  high,  and  e^,  has  dropped  to  30%  of  the  proper 
value. 

The  next  example  illustrates  that  the  unfolding 
process  is  not  necessarily  a converging  process 
in  that  it  may  get  progressively  worse  without  the 
a pnon  constraint  of  a physically  real  spectrum 
for  «,(£).  A comparison  of  the  1975  TRIGA 
J-tube  measurement,  as  unfolded  with  SAND-11 
using  «,(£)  ■ NIOBE  * GAZE  (Table  IV. E),  was 
made  with  some  1973  results  [six-foil  data  un- 
folded with  the  SPECTRA  code  using  o»(£)  » 
constant].  The  value  of  «,y 0 tor  the  1973  spec- 
trum is  0.343.  When  this  spectrum  was  used  as  a 
trial  spectrum  for  a SAND -II  unfolding  code  run 


Fig.  6.  Results  of  unfolding  50-crn  WSMR  FBR  Issksge 
spectrum  with  l/£  moderation  component  missing  from 
#*<£);  i.e..  o»i£)  * "CODIVA”  (solid  line!  Without  gold-foil 
data.  SAND-11  output  igrees  with  glory-hole  spectrum.  When 
low-energy  foil  Jsu  are  added  (gold  foil)  and  with  no  l/£ 
component  in  the  SAND-II  results  ate  a catastrophic 

failure  and  show  the  need  of  good  u pnon  data  (#,(£)! . 


with  the  1975  foil  data,  the  resulting  spectrum 
(Table  IV  N)  yielded  o*,  o » 0 359,  which  shows  a 
further  deterioration.  This  is  compared  to  a,*  a * 
1.019  (Table  IV  E).  Also,  the  total  flux  above 
10  keV  is  65%  high  and  o,,  is  9%  high. 

XI.  SENSITTVmr  TO  UNFOLDING  CODE 

The  WSMR  FBR  glory-hole  data  were  unfolded 
with  the  SPECTRA  code  with  «„(£)  • IDF  calcula- 
tion. The  results,  presented  as  Fig.  7 and  Table 
IV  O,  are  nearly  identical  to  those  of  Fig.  3 and 
Table  IV. C.  The  values  of  o agree,  and  both 
duxes  0 and  e,^  are  the  same  within  1%. 

A comparison  of  Figs.  3 and  7 shows,  however, 
that  the  feature  of  the  SAND-II  code  designed  to 
retard  the  formation  of  spurious  structure  is 
clearly  evident.  The  SPECTRA  result  oscillates 
more  about  the  IDF  calculation  above  0.5  MeV. 


Fig.  ’ Comparison  of  WSMR  FBR  glory -Hot*  spectrum 
unfolded  with  SPECTRA  code  with  e,i£)  ■ I DF  calculation. 
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XII.  VARIATIONAL  STUDY 

A number  of  variational  studies  of  the  effects 
of  cross -section  uncertainties  on  unfolding  have 
been  made  over  the  past  several  years.”  In  some 
of  these  studies,  the  reaction  cross  section  for  a 
given  foil  was  varied  at  random  over  the  neutron 
energy  region  above  threshold.  These  showed  that 
a random  variation  over  many  energy  groups  has 
almost  no  net  effect  on  the  output  spectrum.  In 
the  work  described  here,  the  entire  cross-section 
curve,  a,l£),  was  shifted  bodily.  This  change  can 
be  representative  of  a cross-section  normaliza- 
tion error,  an  incorrect  gamma-ray  branching 
ratio,  or  any  uncertainty  in  the  fission-fragment 
yield  used  in  calculating  back  from  count  rate 
observed  to  activations  actually  produced.  This 
cross -section  change  was  simulated  by  changing 
the  input  activation  number  used  in  SAND- II.  In 
fact,  this  variational  study  simultaneously  evalu- 
ates the  effects  of  counting,  foil  weighing,  and 
self-absorption-correction  errors,  as  well  as 
cross-section  uncertainties- 

Figure  S shows  the  result  of  a series  of 
SAND-II  unfolding  operations  with  nine  foils,  with 
the  cross  sections  varied  by  15  and  25%,  indi- 
vidually for  most  cases,  and  in  pairs  for  the  ”*U 
and  J,Mg  foils.  The  perturbed  spectra  in  Fig.  8 
are  shown  compared  to  the  unperturbed  spectrum 
(solid  line).  The  threshold  energy  of  the  foil 
cross  section  being  varied  is  depicted  by  a verti- 
cal bar  and  horizontal  arrow  located  below  the 
corresponding  curve.  Note  that  the  effect  of  vary- 
ing one  of  the  activations  (or  cross  sections)  is 
very  local  for  those  foils  that  have  thresholus  £, 
above  1 MeV,  where  the  E,  values  are  relatively 
close  together. 

The  variational  data  (Tables  IV.  P through  Z) 
indicate  that  for  foils  with  thresholds  below  2 
MeV,  where  few  thresholds  exist,  a 25%  increase 
in  cross  section  or  foil  activation  produces  ~2% 
decrease  in  Oeq/®l  for  foils  with  thresholds  just 
above  2 MeV,  where  the  neutron  flux  is  still  high 
and  the  thresholds  are  close  together,  a 25%  in- 
crease in  cross  section  of  one  foil  produces  ar. 
— 1%  or  less  increase  in  o«q/«;  for  foils  with 
thresholds  well  above  2 MeV,  a 25%  increase  in 
cross  section  results  in  a 0 to  *0.4%  ®,q/o  change. 

The  flux  above  0.01  MeV  only  varies  by  0 to  3% 
for  25%  variations  in  o except  for  plutonium, 
which  alone  covers  the  0.01-  to  0.5-MeV  region. 
Here,  a 25%  variation  in  cross  section  yielded  a 
10%  variation  in  o.  The  value  of  varies  by  0 to 
4%  for  all  cases  except  plutonium  (8%). 

It  is  clear  that  a flux- measuring  method  em- 


aW  N.  McELROY.  Hanford  Engineering  Development 
Laboratory.  Private  Communication  (197$). 


ploying  a large  number  of  foils  will  provide 
greater  accuracy  than  one  employing  one  or  two 
since  error  averaging  and  even  error  cancellation 
is  involved  with  large  numbers. 

The  threshold-foil  method  of  flux  (and  radia- 
tion-damage) measurement  is  most  sensitive  to 
cross-section  and  counting  errors  in  foils  with 
thresholds  below  2 MeV  because  few  threshold 
values  exist  in  this  energy  region.  However,  this 
energy  region  is  covered  by  fission  foils.  Here, 
cross-section  data  are  now  approaching  5%  in 
accuracy,*  so  that  the  contribution  of  these  foils 
to  the  overall  uncertainty  in  ®,  o,,,  and  o«,/® 
should  be  small  compared  to  the  importance  of  the 
a priori  constraint  of  a good  calculated  value  for 
the  trial  spectrum,  ® „(£),  used  in  the  unfolding. 
Above  a few  MeV,  where  the  cross-section  errors 
are  probably  twice  as  large,  there  are  several 
thresholds  available,  and  appreciable  cancellation 
of  errors  is  expected.  Also,  the  spectrum  is 
rapidly  falling  there,  so  that  the  contribution  to 
®rq  is  small. 

Xm.  SUMMARY  AND  CONCLUSIONS 

XI1I-A.  Capabilities  and  Limitations  of 

Foil  Motivation  Spectrometry 

Threshold-foil  neutron  spectrometry  is  a use- 
ful tool  if  enough  information  is  input  to  the 
unfolding  process.  This  consists  of  the  a priori 
constraint  of  a good  trial  spectrum,  ®w(£),  es- 
pecially below  1 MeV  (but  above  0.01  MeV),  where 
at  most  two  “windows”  (intervals  between  suc- 
cessive £,  values)  exist.  In  addition,  a limited 
amount  of  operator  interaction  may  be  required  to 
detect  and  eliminate  any  spurious  activation  value. 

Thus,  the  threshold -foil  activation  method 
works  well  only  when  it  carries  out  a small  per- 
turbation (correction)  on  a good  spectral  estimate 
to  obtain  ®(£)  that  is  consistent  with  all  threshold- 
foil  data.  The  SAND-11  code  appears  to  be  unique- 
ly qualified  to  produce  Ihe  minimum  necessary 
perturbation,  because  this  improved  version  of 
SAND  retards  the  formation  of  spurious  structure 
while  preserving  the  known  input  structure  (14- 
MeV  peak,  resonance  dip,  etc.).  However,  the 
version  of  SAND-II  used  for  this  work  does  not 
include  the  option  of  assigning  a varying  statisti- 
cal weight  to  each  foil  activation,  and  this  is  a 
serious  shortcoming  that  should  be  corrected. 

XUIB ■ Accuracy  and  Precision 

The  characterization  of  a neutron  field  in 
terms  of  its  silicon  radiation  damage  effective- 


HProc  Consultants  Mtg.  on  Nuclear  Data  for  Reactor  Neu- 
tron Dosimetry.  INDC  (NDSV56/U,  September  10-12,  1973, 
International  Atomic  Energy  Agency,  Vienna  (1973). 
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9 FOILS  

STANOARO 


’"Pu  X 1.15 
»*Pu  x I 25 


n,U  < 1 15 
■>•0  x 1.25 


”S  < 1.15 
”S  x 1 25 


;4Mg  a 1.25 


“PeXI.25 


a*U  x 1.25  ANO 
’•Mg  a 1.25  — 


n»U  x 1.25  ANO 
’•Mg  X 0.75 


NEUTRON  ENERGY  (VeV) 

Fi*.  t.  Effect  of  virying  foil  activation  or.  equivalently,  the  cross  section  on  SAND-11  spectra.  Solid  curve  ts  for  unperturbed 
data  m each  case  (the  “standard”  run). 


ness  is  accomplished  with  aa  extensive  param- 
eter, 0^  and  an  intensive  parameter,  o^/'o.  The 
former  quantifies  the  total  radiation  damage  m 
terms  of  the  equivalent  fluence  of  1-MeV  (*15%) 
neutrons,  and  the  latter  quantifies  the  damage- 
producing  capability  of  a unit  of  fluence,  again  in 
terms  of  the  equivalent  1-MeV  (*15%)  fluence. 


The  definition  of  0«s  and  0*^0  in  Eqs.  (3)  and 
(4)  involves  only  o(£)  and  DiE).  Since  DiE ) ap- 
pears in  the  numerator  and  D (1  MeV  * 15%)  in  the 
denominator,  only  the  shape  of  D(£)  is  important 
for  specifying  the  accuracy  limitations  of  the 
present  means  of  characterizing  the  neutron  field. 
Thus,  assuming  the  shape  of  D(E)  is  correct  [the 
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determination  of  the  accuracy  of  D(£)  is  beyond 
the  scope  of  this  paper],  this  leaves  the  measure- 
ment of  the  spectral  flux  o(£)  as  the  controlling 
factor  in  characterizing  the  radiation  damage 
effectiveness  of  the  neutron  field.  The  accuracy 
of  determining  the  spectral  flux  depends  mostly  on 
the  accuracy  of  the  foil  activation  measurements 
and  data  processing,  on  the  foil  cross  sections  and 
gamma- ray  branching  ratios,  and  on  the  trial 
spectrum  o „(£),  used  in  the  unfolding  codes. 

The  cross-section  set  is  probably  known  to  an 
accuracy  of  10  to  15%  for  nonfission  foils  and 
approaches  5%  for  the  fission  foils’9  that  provide 
nearly  all  the  spectral  data  from  0.01  to  2 MeV. 
The  foil-activation  measurements  and  branching 
ratios  are  given  a suggested  uncertainty  of  15%  or 
smaller.  If  a large  number  of  foils  are  used,  the 
variational  study  therefore  indicates  that  can 
be  obtained  to  an  accuracy  of  -10%  and  o«q/o  to 
-5%.  This  was  independently  checked  by  unfolding 
foil  measurements  of  several  ‘‘known’’  spectra 
with  os  * oCmi coined-  The  calculations  (oc)  were 
known  to  better  than  5%  in  terms  of  o«,/o. 

The  precision  of  the  threshold-foil  method 
depends  mostly  on  foil-measurement  accuracy.  It 
Is  only  slightly  better  than  the  accuracy  because 
the  precision  or  repeatability  is  probably  the 
limiting  factor  in  determining  thp  accuracy  of 
specifying  o(£),  and  o«q/<a,  once  the  constraint 
of  a good  starting  spectrum,  o„(£),  is  accepted 
and  applied. 

X1U.C.  Spectral  Index  and  Other  Methods  of 
Characterizing  Silicon- Damage 
Effectiveness  of  a Neutron  Field 

Once  the  spectral  shape  o(£)  of  a neutron  field 
has  been  determined,  the  absolute  fluence  can  be 
measured  accurately  with  any  single  threshold- 
activation  foil  whose  threshold  value,  £,,  is  low 
enough  to  sample  the  energies  of  importance  yet 
high  enough  to  receive  negligible  activation  from 
the  unimportant  thermal  and  epithermal  flux.  For 
the  three  reactor  spectra  presented  here,  the 
fraction  of  the  total  damage  sampled  by  a boron- 
covered  2*U  or  2"Pu  foil  is  100%  (£,  =»  0.01  MeV), 
"Np  m 85%  (£,  » 0.5  MeV),  2MU  “ 50%  (£,  * 
1.45  MeV),  and  "S  * 15%  (£,  * 2.9  MeV).  While  a 
boron-covered  2MU  or  23®Pu  foil  can  give  a good 
measure  of  and  0,,/®,  because  of  a low  thresh- 
old energy  and  flat  cross  section  up  to  several 
MeV,  the  common  practice  of  using  °S  foils 
(because  of  handling  and  licensing  convenience  and 
the  economy  of  the  required  beta  counter)  to 
estimate  the  silicon  radiation  damage  can  lead  to 
appreciable  error  If  the  neutron  spectrum  is  not 
accurately  known.  Calculated  spectra  have  often 
been  used  to  obtain  a (£)  which,  in  turn,  was  used 
to  obtain  the  plutonium-to-sulfur  ratio  (spectral 


index  of  Table  IV).  However,  most  early  calcula- 
tions have  underestimated  o(£)  at  energies  above 
2.9  MeV.  Thus,  the  spectral-index  characteriza- 
tion can  lead  to  large  errors  in  estimating  the 
total  neutron  damage,  o,q,  unless  the  spectral 
shape  is  accurately  known.  Note  that  in  the 
variational  study  (spectra  P through  2,  Table  IV), 
the  SI  varied  by  ~2f%,  while  the  range  of  variation 
of  ®«,  and  ®e,/®  was  -8%. 


APPENDIX 

THE  STANDARD  SPECTRA 

The  WSMR  FBR  spectrum  (Figs.  2 and  3)  was 
calculated  in  the  glory  hole  and  at  a number  of 
points  outside  the  FBR  with  the  IDF  code”  (a 
modification  of  the  DTF-IV  code).2*  The  ENDF/ 
B-IV  cross  sections  were  used,  and  the  slant- 
back-scattering  from  the  floor  (-2  m away)  was 
simulated  by  a 4ir  concrete  wall,  18  cm  thick 
and  12  m away.  The  20-cm-diam  x 19-cm-high 
core,  made  of  10%  Mo— 93.2%  23*V-enriched  urani- 
um fuel,  was  homogenized  with  stainless -steel 
structural  material  into  a spherical  geometry. 
The  flux -weighted  cross  sections  were  calculated 
with  the  GGC-5  code.2*  The  IDF  calculations 
utilized  S,  quadrature,  P,  scattering,  30  energy 
groups,  and  6 material  zones  [glory  hole;  air; 
liner  (stainless  steel);  core  (93%-23*U-enriched  ♦ 
10%  Mo  -i-  homogenized  stainless -steel  bolts); 
shroud  (thin  ‘°B  A1  * resin);  air;  and  4ir  concrete 
wall]. 

A recent  IDF  calculation  for  a split-bed  sub- 
critical  reactor  assembly  utilizing  ENDF/B-IV 
cross  sections  is  shown  in  Fig.  A.l.  It  is  in  very 
good  agreement  with  the  accurate  TOF  (Ref.  7) 
and  the  proton-recoil  proportional  counter  mea- 
surements,’,  above  0.01  MeV  where  nearly  all  of 
the  radiation  damage  occurs  (see  Fig.  1). 

The  TRIGA  reactor  "standard  spectrum”  was 
taken  as  a composite  spectrum  of  a fine-group- 
structure  NIOBE  code  calculation11  above  0.8  MeV 
(verified  experimentally  with  a *Li-7Li  diode- 
sandwich  spectrometer  utilizing  shielded  diodes11) 
overlapped  with  the  GAZE -code  reactor  criticality 


nThe  IDF  calculations  and  the  GGC-5  calculation  of  the 
flux-weighted  cross  sections  were  earned  out  by  C Rind  flench 
and  D.  Houston  at  IRT  Corporation,  San  Diefo,  California. 

aK.  D.  LATHROP.  "DTF-IV  Code,  A FORTRAN-IV  Pro- 
gram for  Solving  the  Multi-Group  Transport  Equation  with 
Anisotropic  Scattering,''  LA-3373,  Los  Alamos  Scientific  Lab- 
oratory (1964).  (IDF  is  an  IRT  modification  of  the  DTF-IV 
code.) 

”D.  R.  MATHEWS  et  al.,  “GGC-5,  A Computer  Program  for 
Calculating  Neutron  Spectra  and  Group  Constants,"  GA-8871, 
General  Atomic  Company  (1971). 
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NEUTRON  ENERGY  ItV) 

Fig.  A.I.  rime-of-fll*hr  and  proton-recod  tpectra  lor  Linac-pulstd  lubcnticii  assembly  compared  to  I OF  calculations  iENDF' 
B-IV  cron  sections). 


calculation11  In  the  0.3-  to  3-M«V  region.  The 
GAZX-code  reaults  were  used  below  0.3  MeV. 
The  GAZE  spectrum  was  checked  Independently 
with  ‘"Pu  (l#B-covered),  '"Np,  *"U,  and  "S  acti- 
vation, with  respective  neutron  threshold  energies 
of  -0.01.  0.9,  1.45,  and  2.9  MeV  (Ref.  30).  A 
comparison  of  messured  neutron  Integral  fluxes 
above  each  of  the  four  respective  thresholds  with 
those  given  by  the  GAZE  calculation  yielded  a 
constant  ratio  (within  j ■ 1.4 %). 


*1.  A.  SAYEG.  Compiler,  "Neutron  uid  Gemma  Dosim- 
a tty  Measurements  at  the  AFRRl-DASA  TRICA  Reactor," 
AFFRI  CR6J-6.  Air  Force  Radiobiology  Reaearch  Iniutute, 
Defense  Atomic  Support  Agency  1 1965). 
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DESIGNATION:  E XXI 


Standard  Method  for 

IRRADIATING  A STANDARD  SET  OF 
NEUTRON  THRESHOLD  ACTIVATION  FOILS 
FOR  RADIATION  HARDNESf  TESTING 


1 . SCOPE 

1.1  This  method  describes  a procedure  for  irradiating 
a standard  set  of  threshold  activation  foils  to  be  used  for 
neutron  spectrum  unfolding  covered  in  ASTM  Method  E XX3, 
"Standard  Method  for  Unfolding  Neutron  Spectra  for  Radiation 
Hardness  Testing."  It  is  intended  to  be  used  in  conjunction 
with  ASTM  Method  E XX2,  "Measuring  Foil  Activities  for  Radia- 
tion Hardness  Testing,"  for  providing  the  foil-activation  data 
used  as  input  to  the  SAND  II  neutron  spectrum  unfolding  code. 

1.2  This  method  presents  a standard  set  of  foils  that 
have  been  used  at  many  facilities,  and  describes  the  flux- 
uniformity,  the  neutron  self-shielding,  and  flux-depression 
corrections  that  need  to  be  considered  in  choosing  the  foil 
thickness,  the  covers,  and  the  locations  of  the  foils. 

1.3  In  this  method,  considerations  that  apply  rather 
generally  to  neutron-activation  detectors  are  discussed  in 
ASTM  Method  E 261,  "Measuring  Neutron  Flux  by  Radioactivation 
Techniques . " 

1.4  Background  information  on  the  detailed  methods  for 
individual  threshold  foil  detectors  of  widest  use  is  given  in 
the  following  ASTM  Methods. 
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E 262  Measuring  Thermal  Neutron  Flux  by 
Radioactivation  Techniques 

E 263  Measuring  Fast-Neutron  Flux  by 
Radioactivation  of  Iron 

E 264  Measuring  Fast-Neutron  Flux  by 
Radioactivation  of  Nickel 

E 265  Measuring  Fast-Neutron  Flux  by 
Radioactivation  of  Sulfur 

E 266  Measuring  Fast-Neutron  Flux  by 
Radioactivation  of  Aluminum 

E 343  Fast  Neutron  Flux  by  Analysis  of 
Molybdenum-99  Activity  from 
Uranium-238  Fission 

E 393  Measuring  Fast-Neutron  Flux  for 
Analysis  for  Barium-140  Produced 
by  Uranium- 2 38 

E 419  Selection  of  Neutron  Activation 
Detector  Materials. 

2 . APPARATUS 

2.1  Precision  balance  (0.1  mg  accuracy). 

2.2  Cadmium  covers  (0.05  to  1 mm  thick). 

2.3  cover  (''•1-1. 8 g/cm^  ^®B) 

2.4  A 1.27  cm  (0,5  inch)  i.d.  aluminum  "cup"  with 
cover  for  ‘■*■'1  (as  ICH^CONHj)  . 

3.  THRESHOLD  ACTIVATION  FOILS 

3 . 1 Redundancy 

3.1.1  The  set  of  foils  listed  in  Table  I is  part  of 
the  self-consistent  set  that  resides  in  the  SAND  II  cross- 
section  library.  These  foils  were  selected  in  an  experimental 
evaluation  program  (1)  in  which  three  different  "known"  spectra 
were  compared  to  corresponding  threshold-foil  spectrometry 
measurements  unfolded  with  the  SAND  II  code.  This  set 


contains  two  redundant  1/v  foils  {all  three  should  be 

exposed,  counted,  and  input  to  SAND  II  for  a first-pass 

(ACTIVATION)  run  to  select  the  average  one) , and  a redundant 

E..  ■ 0 fission  foil.  233U  and  23^Pu  have  very  similar 
c 235 

cross-section  shapes.  The  U foil  presents  much  less  of  a 

239 

safety  hazard  than  Pu,  and  is  cheaper.  It  is  very  useful, 
when  measuring  soft  (TRIGA)  spectra,  in  determining  the  cor- 
rection for  the  235U  impurity  in  the  238U  foil  (readily  avail- 
able with  ^400  ppm  *'33U  impurity).  Although  the  32S  foil  is 
widely  used  as  a monitor  foil,  it  is  not  part  of  the  set  be- 
cause it  requires  an  entirely  different  counter  (for  betas) 

and  an  involved  calibration  technique.  It  has  about  the  same 
58  58 

threshold  as  the  Ni(n,p)  Co  foil  react  and  should  only 
be  used  when  the  AQ  (the  SAND  II  activation  error"  as  cal- 
culated for  #tr(E)  as  input:  see  ASTM  Method  E XX3)  of  the 

two  are  within  a few  percent (^5  percent)  of  each  other;  other- 
wise, SAND  II  may  be  driven  to  many  iterations  in  trying  to 
find  a solution  $(E)  that  is  compatible  with  both  activations. 
This  usually  leads  to  an  unacceptable  solution. 

3 . 2 Foil  Impurities 

3.2.1  Foil  impurities  are  especially  serious  for  a 

moderated  source  (TRIGA  reactor)  when  the  impurity  leads  to 

the  same  reaction  product  by  way  of  thermal -neutron  capture. 

Some  of  these  foils  (impurities)  are  238U(233U) , 27A1(  3Na) , 

^Fe(3^Mn) . For  a soft  spectrum,  such  as  the  TRIGA  J-tube 

235 

spectrum  (boral  shielded) , the  U fissions  (Cd  covered) 

238  238 

were  »100  times  greater  than  the  U;  therefore,  the  U 

235 

must  have  no  more  than  ^200  ppm  of  U to  reduce  the  error  to 

2 percent.  Higher  impurities  can  be  tolerated  for  GODIVA- 

type  reactors  where  the  low-energy  flux  is  much  less  Intense, 

235 

or  with  TRIGA-type  reactors  if  the  U foil  data  are  used 

238 

for  correcting  the  U activity.  In  this  case,  the  percent 
4“3<*U  in  238U  must  be  accurately  known.  For  the  38Fe(n,p)^Mn 
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reaction  the  33Mn  impurity  must  be  no  more  than  10  ppm 
(with  Cd  cover)  for  use  with  a TRIGA  spectrum,  and  100  ppm 
at  50  cm  from  a GODIVA-type  reactor  (^2  meters  off  the  con- 
crete floor)  for  42  percent  background  from  the  33Mn(n , y)  ^Fe 
reaction.  Similarly,  a manganese  foil  (Cd  covered)  can  be 
used  to  correct  the  ^Fe  data  if  the  impurity  correction  is 
small  (i20  percent  of  total  n,p  activation)  and  the  percent 
Mn  in  Fe  is  accurately  known. 

3.3  The  Influence  of  Nuclear  Data  on  the 


3.3.1  Since  the  total  number  of  interactions  must  be 
deduced  from  the  absolute  ganssa-'-a/  count  with  good  accuracy 
(say  5 percent  per  foil) , the  foils  selected  must  have  gamma- 
ray  yields  known  to  a similar  accuracy  or  better.  These  cor- 
rections include  conversion-electron  production,  branching 
ratio  to  a given  energy  level,  and  fission  yield. 

1 40 

3.3.2  The  1593-keV  gamma-ray  line  from  La  produced  in 

232  232 

Th  fission  is  not  useful  because  of  interference  from  Th 

radioactivity.  This  has  often  led  to  the  use  of  the  537-keV 
140  140 

line  from  the  Ba  precursor  of  La,  which  is  listed  in  the 

Table  of  Isotopes  (2)  as  having  an  intensity  of  0.34  gamma  per 

^®Ba  decay.  A recent  evaluation  (3)  has  shown  this  to  be 

25.7  percent,  which  is  in  much  closer  agreement  with  some  work 

235 

of  Ref.  1,  in  which  U fission  yields  were  compared  by  wav 

140 

of  four  fission-product  gamma  rays  (1593  keV  for  La,  537  keV 
for  U0Ba,  743  keV  for  97Zr,  and  668  keV  for  132I)  . 

3.3.3  The  choice  of  gamma-ray  line  thus  directly  influences 

the  accuracy  of  determining  the  specific  activity  produced 

during  the  neutron  irradiation.  It  also  influences  the  final 

choice  of  foil  thickness,  in  that  the  selection  of  a low- 

energy  gamma-ray  line  may  lead  to  a large  gamma-ray  self- 

232 

absorption  correction  in  counting.  For  example,  the  Th 
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foil  of  Table  I has  a maximum  attenuation  of  22  percent, 
or  an  average  correction  of  ^11  percent,  for  the  537  keV 
line.  This  represents  an  upper  limit  for  foil  thickness. 

Thus,  the  gamma  scattering,  as  well  as  the  neutron  self- 
shielding discussed  below,  will  influence  the  foil  selection. 

3.4  Foil  Encapsulation 

3.4.1  Fission  foils  may  be  encapsulated  in  a hermetic- 
ally sealed  container  to  avoid  oxidation  and  loss  of  materials , 

and  for  health-hazard  requirements.  The  23^Pu  foil,  if  used 

23S  - 

instead  of  the  much  safer  U foil,  will  require  special 
encapsulation  and  periodic  wipes  to  check  for  leakage  of  the 
material.  Copper  encapsulation  has  been  found  satisfactory 
for  235U,  238U,  237Np,  and  232Th  foils.  It  is  made  -v-O.l  to 
0.25  mm  thick  at  the  flat  surfaces  and  is  soldered  at  the 
periphery . 

3.5  Foil  Diameter 

3.5.1  All  foils  are  to  be  1.27  cm  (0.5  inch)  or  smaller 

in  diameter  to  simplify  foil-size  effects  in  calibrating  agains 
point  gamma-source  standards  (see  ASTM  Method  E XX2) . A cor- 
rection (small)  can  be  made  for  this  geometric  effect. 

4.  IRRADIATION  PROCEDURES 

4.1  Foil  Covers 

4.1.1.  Cadmium  covers  of  0.5  ran  thickness  are  prescribed 
for  all  fission  foils  and  1/v  detectors,  and  for  detectors  such 
as  238U,  38Fe,  8®Ni,  and  27A1,  where  traces  of  impurities 
(235U.  33Mn,  8®Co,  and  23Na)  that  yield  the  same  reaction 
product  via  thermal-neutron  capture  can  lead  to  unmanageably 
large  corrections.  When  these  corrections  are  greater  than 
5 percent,  the  irradiation  should  include  cadmium-covered  foils 
made  of  these  impurities  (233U,  33Mn,  6®Co,  23Na,  etc.).  The 
corrections  can  then  be  made  with  good  accuracy  if  the  percent 


impurity  in  the  foil  is  accurately  known.  If  not,  a 
thermal-neutron  activation-analysis  type  of  irradiation  will 
be  required.  Cadmium  foils  are  not  needed  in  the  glory  hole 
of  a fast-burst  reactor  with  little  or  no  moderator  material 
(<0.5  g,  say)  near  or  inside  the  cavity. 

4.1.2.  Covers  of  ^B  are  useful  on  fission  foils  for 

the  case  of  measuring  a soft  TRIGA  spectrum,  especially  if 

a boral  shield  is  used  to  surround  the  irradiation  cavity. 

If  the  boral  shield  is  a good  4tt  shield,  and  if  a negligible 

amount  of  moderator  is  contained  within  the  shield,  the  boral 

shielding  can  be  properly  accounted  for  bv  the  SAND  II  unfold 

ing  code.  In  this  case,  the  cover  is  not  as  important. 

The  boral  thickness  is  entered  for  each  and  every  foil,  and 

a 1/E  spectrum  is  entered  down  to  3 x 10"'  MeV , the  cadmium 

cutoff  energy;  or  lower.  If  no  cover  is  used,  and  if  the 

cavity  is  only  partially  shielded,  it  will  be  difficult  to 

predict  the  neutron  spectrum  from  10'^  MeV  down  to  ^3  x 10” ' 

MeV,  where  the  cadmium  covers  become  effective.  In  this  case 

it  is  important  to  place  all  the  1/v  foils,  the  foils  with 

important  1/v  impurities  (see  above) , and  all  fission  foils 

in  a boral  "box”  or  a ^B  cover.  For  best  results,  a ^B 

2 10 

cover  of  1 to  1.8  g/cm  of  (93  percent)  B is  used.  This 

cover  thickness  is  then  input  to  SAND  II.  In  this  wav,  the 

-7  -2 

fraction  of  activations  arising  from  3 x 10  MeV  to  10  “ MeV 
neutrons  will  be  both  greatly  reduced  and  more  accurately 
calculated  for  the  SAND  II  unfolding  process.  The  known 
spectrum  outside  the  boral  (a  1/E  spectrum)  is  used  down  to 
3 x 10‘7  MeV  or  lower. 

4.1.3  The  ^B  covers1  may  be  replaced  by  cadmium  covers 
for  up  to  one  meter  from  a GODIVA-tvpe  reactor  that  is  a 
few  meters  above  the  concrete  floor,  or  for  the  glcrv  hole 
where  no  low-energy  neutrons  are  found,  and  where  the  ^B 


covers  generally  cannot  be  placed  in  any  case.  If  the 

covers  are  used  with  a directional  source,  such  as  outside  . 

the  GODIVA  reactor,  the  fission-foil  activation  will  require 

a correction  for  scattering  by  ^°B.  The  correction  can  either 

be  done  experimentally,  with  pure  finite-threshold  fission 
237  232 

foils  ( Np  or  Th)  that  contain  negligible  zero- threshold 

impurities  which  yield  the  same  gamma-rav  lines,  or  with  a 

sophisticated  calculation  that  weights  each  inscattering  by 

the  track  length  through  the  material.  These  corrections 

2 10 

are  the  order  of  10  percent  for  a 1.65  g/cm  B cover  and 
a thin  1/2-inch-diameter  fission  foil  (4) . 

4. 2 Foil  Interferences 

O35 

4.2.1  A strong  resonance  absorber  such  as  a thick  “ U 
foil  cannot  be  placed  in  front  of  a 1/v  absorber,  and  thick 
foils  with  their  covers  should  not  be  stacked  so  as  to  result 
in  a large  and  unmanageable  scattering  correction.  For  iso- 
tropic flux,  the  interfering  foils  should  not  scatter  more 
than  about  10  percent  of  the  flux,  as  given  by  the  simple 
hand  calculation: 


(1) 


where  ItT  is  the  sum  of  the  products  of  the  macroscopic  total 
cross  sections  and  thicknesses  for  those  foils  that  are  stacked 
in  front  of  the  deepest  foil,  i.e., 


£tT  - St(l)T1  + Et(2)T2  + • • • ■ (2) 

For  a beam  geometry,  che  corresponding  "interference"  must  be 
held  to  no  more  than  ^4  percent  because  of  the  more  complex 
scattering  correction. 
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4.3  Foil  Self  Shielding 

4. 3. 1.1  The  correction  for  self  shielding  is  appreciable 
only  for  the  0 . 0025-cm-thick  gold  foil  (with  its  highly  absorb- 
ing resonance  at  "^5  eV)  , being  about  a factor  of  two  for  epi- 
cadmium  neutrons  (5,6).  The  effect  on  the  unfolded  spectrum, 
♦(E),  of  varying  the  gold-foil  activation  by  a factor  of  two  was 
mostly  local,  appearing  as  a dip  or  bump  at  the  5-eV  resonance, 
and  changing  ♦ by  only  1 percent. 

4 . 4 Flux  Uniformity 

4.4.1  If  the  foils  cannot  all  be  located  in  the  same 
region,  or  in  a region  of  uniform  flux  (as  determined  from 
symmetry  considerations) , they  can  be  spread  out  over  a larger 
volume  of  varying  flux  but  of  constant  neutron-spectrum  shape. 

If  the  flux  varies  by  more  than  3 percent  from  point  to  point, 
flux  monitors  should  be  used  at  the  various  cases.  Around 

a GODIVA-tvpe  reactor,  thick  sulfur  foils  can  be  located  near 
individual  foils.  Where  space  is  more  limited,  thin  nickel 
{^Ni(n,p)^Co}  , iron  { ^Fe(n , p)  ■’Sin)  , or  even  aluminum 
{^'Al(n,a)  - wa}  , monitors  can  be  used;  considerations  dis- 
cussed in  ASTM  Method  E 261  apply. 

4 . 5 Flux  Depression 

4.5.1  At  low  energies,  flux  depression  is  important 

for  bare  thermal-neutron  detectors  near  cadmium-covered  discs, 
if  both  are  embedded  in  a moderator.  At  high  energies,  it  is 
important  for  the  same  situation  if  the  moderator  contains 
reactor  fuel.  However,  in  sizeable  cavities  used  for  radiation 
damage  studies,  the  cadmium  covers,  as  well  as  the  foils, 
generally  subtend  a negligibly  small  solid  angle  at  the  point 
of  any  surrounding  moderator  of  fuel.  For  a GODIVA  reactor 
glory  hole,  the  foil-interference  considerations  discussed 
in  Sec.  4.2  will  completely  dominate. 
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FOIL  WEIGHING 


5.1  All  foils  should  be  weighed  with  a 0.1-mg  precision 
balance  that  is  checked  against  known  weights  over  the  range 
of  the  foil  masses.  The  encapsulated  foils  are,  of  course, 
carefully  weighed  before  encapsulation. 

6.  PURITY  CERTIFICATION 

6.1  The  foil  purity  analysis  results  must  be  kept  on 

permanent  record  for  use  in  making  foil  impurity  corrections. 

These  impurities  must  be  known  to  an  accuracy  dictated  by  the 

magnitude  of  the  correction,  which,  in  turn,  will  depend  on 

the  kind  of  neutron  spectrum  being  measured  (see  Sec.  3.2). 

235 

If,  for  example,  the  percentage  impurity  (say  400  ppm  U 
^38 

in  “ U)  is  known  to  an  accuracy  of  10  percent  in  a foil,  and 

235 

the  separate  impurity  foil  ( U)  is  irradiated  the  same  way 

as  the  other,  then  the  impurity  effect  can  be  reduced  to 

235  '>38 

10  percent  of  its  seated  value  (40  ppm  U in  ~ U,  for  this 

example)  bv  calculating  the  correction.  In  this  case,  up  to 
23S 

2000  ppm  of  'U  impurity  could  be  tolerated  for  a TRIGA 
spectrum  (see  Sec.  3.2). 

7.  PRECISION 

7.1  Foil  weighing,  neutron  self  shielding,  and  flux 
uniformity  influence  the  precision  of  those  aspects  of  the 
flux  measurement  to  which  this  method  is  addressed.  An 
estimated  value  of  3 percent  is  suggested  as  the  effects  of 
these  corrections  on  the  measurement  of  specific  activities, 
assuming  the  scattering  of  the  stack  is  kept  to  less  than 

10  percent  for  an  isotropic  flux,  and  to  less  than  4 percent 

for  a directed  flux.  Assuming  no  thick  covers  are  used 

in  the  directed-f lux  case,  and  secondary  monitors  such  as 
58  58 

thin  Ni(n,p)  Co  reaction  foils  are  used  when  the  foils 
are  widely  distributed  in  space,  these  factors  will  have 
only  about  a 1 percent  effect  on  $eq- 
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8 . ACCURACY 


8.1  The  factors  that  affect  the  precision  (see  Sec.  7) 
will  directly  impact  the  accuracy.  In  addition,  the  selection 
of  the  foils  and  the  gamma-ray  line  to  be  measured  for  each 
foil  impact  the  accuracy  by  way  of  the  uncertainties  in  gamma- 
ray  i".censity  (gamma/reaction)  and  fission  yield.  These  may 
be  as  arge  as  10  percent  in  individual  cases,  but  should,  on 
the  average,  contribute  no  more  than  5 percent  to  the  uncer- 
tainty in  the  specific  activations.  Assuming  a reasonable 
amount  of  foil-to-foil  randomness  in  this  error,  the  contri- 
bution to  the  uncertainty  of  specifying  <>e(^  is  estimated  to 
be  no  more  than  2 to  3 percent. 
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TABLE  I 


Activation  Foils  (1.27-cm  diam) 


Reaction 

E, 

IMeV) 

Gamma/ Reaction 
(Fission  Yield.  %) 

Foil  Maaa 

(8) 

Isotopic 

Abundance 

R) 

Notes 

‘"Aui*.  vi^Au 

0 

412 

0.95 

2.696  days 

0.056 

100 

A 

"CoU.yf’Co 

0 

1173 

1.00 

5.258  yr 

0.057 

100 

1333 

1.00 

"Ma(«,y)"MB 

0 

847 

0.99 

2.58  h 

0.030 

100 

a.b 

0.29 

3*l){*./)‘'°La 

0<0.01)c 

0.96  (6.29) 

40.23  h 

0.281 

100 

A,d.* 

a*Pu(e,/)1<0La 

0(0.01)c 

0.96  ( 5.24) 

40.23  h 

0.150 

100 

Aed.* 

"Npin./l^La 

0.5 

1596 

0.96  (5.69) 

40.23  h 

0.580 

100 

A,d,« 

3TNp(i.,/)"Zr 

0.5 

743 

0.92  (6.01) 

16.8  h 

0.530 

100 

A,f,« 

‘“ln(«.a')lis"ln 

1.1) 

335 

0.50 

4.50  h 

0.255 

95.7 

a*U(«t./)‘“La 

1.45 

1593 

0.96  (6.02) 

40.23  h 

0.495 

100 

A.d,« 

“Thfx./J^Ba 

1.75 

537 

0 256  (7.91) 

12.8  days 

1.066 

100 

*.|.e 

“Th  <*./i"Zr 

1.75 

743 

0.92  (4.12) 

16.8  h 

1.066 

too 

a /.* 

MFe(i«.*>14Mn 

2.20 

835 

1.00 

303  days 

0.142 

5.82 

**Nl  (».;>)5,Co 

2.9 

ato 

0.99 

71.3  days 

0.282 

67.8 

nS  <*,/>)”  P 

2.9 

3etas 

1.000 

14.3  days 

4.10 

95.0 

a 

**Mg(it,^)MNa 

6.3 

1369 

1.00 

15.0  h 

0.030 

79 

A 

**Fe  (»./>)MMn 

7.5 

347 

0.99 

2.58  h 

0.142 

91.7 

Kt 

1811 

0.29 

”Al  (a.  a)MNa 

8.7 

1369 

1.00 

15.0  h 

0.257 

100 

A 

2i»)l*I 

11.0 

386 

0.34 

13.05  days 

0.657 

100 

667 

0.33 

"Zrl*,  2n)*Zr 

910 

• 0.99 

78.4  h 



0.108 

SI. 5 

'Cadmium  cover  0.05  to  0.10  eta  thick. 

‘*Co  instead  o(  ‘"au  and  uMn  for  very  long  irradiations. 

c E,  * 0.01  MeV  with  “B  sphere  [important  for  soft  (TRICA)  spectra,  where  <#( E)  < 0.01  will  otherwise  dominatel. 
•Vhen  “U  or  °*?u  foil  is  108-covered.  also  cover  aU  and  ^Np  foils  so  that  accurate  corrections  can  Oe  made  for 
"U  and  3*Pu  impurities  in  these  high  E,  fo.ls. 

d40.23-b  daughter  of  12. 80-day  ‘“Ba.  iVatt  five  days  for  maximum  decay  rate  (ASTM  E 303). 

'Fission  yields  are  for  bombardment  with  fission-spectrum  neutrons.  For  thermal  and  It  MeV,  see  Ref.  5. 
fUse  "Zrfor  low  fluence  (3  * 10u  to  3 x tO^n/cra1).  Use  peak  shape  analysis  or  measure  twice,  seven  days  apart, 
to  strip  off  740-keV  *Mo  gamma  ray  (rl/#  » 67  h). 

*3,'I7i  radioactivity  interferes  with  140 La  line. 

h Requires  separata  detector,  and  calibration  technique  is  complex. 

‘Maximum  manganese  impurity  » 0.001%.  cadmium -cove  red.  Omit  MFe  (s.pi^Mn  for  long  irradiations. 
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DESIGNATION:  E XX2 


Standard  Method  for 

MEASURING  FOIL  ACTIVITIES  FOR 
RADIATION  HARDNESS  TESTING 


I.  SCOPE 

1.1  This  method  describes  a standard  procedure  for 
measuring  the  absolute  gamma-ray  emission  rate  from  a stan- 
dard set  of  neutron  threshold-activation  foils  used  in  spec- 
trum measurements,  and  for  calculating  the  specific  activity 
of  the  foil  from  the  gamma-rav  data. 

1.1.1  The  measuring  procedure  takes  into  account  cor- 
rections for  finite  foil  size  and  thickness  in  detector  cali- 
bration, count  rate  and  pulse-pileup  losses,  and  background 
measurements  for  complex  decay  schemes  such  as  exist  in  some 
fission-foil  lines. 

1.1.2  The  data  processing  takes  into  account  complex 
background  subtraction  for  fission-foil  lines,  corrections 
for  irradiation,  wait  and  count  times,  and  corrections  for 
gamma-ray  branching  ratios,  conversion  electrons,  fission 
yields,  and  gamma-ray  self  absorption  in  the  foils. 

1.2  This  method  is  intended  to  be  used  in  conjunction 
with  ASTM  Method  E XXI,  "Irradiating  a Standard  Set  of  Neut- 
ron Threshold  Activation  Foils  for  Radiation  Hardness  Test- 
ing," for  providing  the  foil-activation  input  data  that  are 
required  in  unfolding  neutron  spectra,  as  described  in  ASTM 
Method  E XX3,  "Unfolding  Neutron  Spectra  for  Radiation  Hard- 
ness Testing."  This  method  is  part  of  the  set  of  draft 
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standards  which  include  ASTM  Method  E XX4,  "Characterizing 
Neutron  Spectra  in  Terms  of  1-MeV  Equivalent  Fluence  for 
Radiation  Damage  in  Silicon,"  and  ASTM  Method  E XX5 , "Measur- 
ing the  Relative  1-MeV  Silicon-Equivalent  Fluence  with  Fast 
Neutron  Monitors." 

1.3  In  this  method,  considerations  that  apply  rather 
generally  to  neutron-activation  detectors  are  discussed  in 
ASTM  Method  E 261,  "Measuring  Neutron  Flux  by  Radioactivation 
Techniques,"  while  background  information  on  the  detailed 
methods  for  the  individual  threshold-foil  detectors  of  widest 
use  is  given  in  the  following  ASTM  Methods. 

E 262  Measuring  Thermal  Neutron  Flux  by 
Radioactivation  Techniques 

E 263  Measuring  Fast-Neutron  Flux  by 
Radioactivation  of  Iron 

E 264  Measuring  Fast-Neutron  Flux  by 
Radioactivation  of  Nickel 

E 265  Measuring  Fast-Neutron  Flux  by 
Radioactivation  of  Sulfur 

E 266  Measuring  Fast-Neutron  Flux  by 
Radioactivation  of  Aluminum 

E 343  Fast-Neutron  Flux  by  Analysis  of 
Molybdenum-99  Activity  from  238u 

E 419  Selection  of  Neutron  Activation 
Detector  Materials 

2 . APPARATUS 

2 . 1 Ge(Ll)  or  Intrinsic  Germanium 
Camma-Ray  Detector 

2.1.1  Ten  percent  or  greater  photopeak  efficiency  for 
6^Co , referred  to  a 7.62  cm  by  7.62  cm  (3  inch)  diameter 
Nal(Tl)  detector,  both  at  25  cm  from  source.  Multichannel 
analyzer,  about  4000  channels  total,  100  mHz  or  faster, 

2 to  2.5  keV  resolution.  Precision  tail-pulser,  ‘'-60  cps. 
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with  inputs  to  gamma-ray  detector  preamplifier  and  to  scaler. 

2 . 2 Foil  and  Source  Holders 

2.2.1  For  accurately  positioning  (and  measuring)  the 
center  of  each  foil,  and  the  center  of  the  gamma-rav  standard. 
Required  precision  is  ^0 . 2 mm  or  better  in  distance  from  face 
of  detector,  and  0.5  mm  or  better  in  lateral  alignment.  The 
apparatus  must  have  provision  for  inserting  a 1 . 27-cm-thick 
lead  shield  (see  Sec.  2.3). 

2.3  Lead  Slab 

2.3.1  A lead  slab  1.27  cm  thick,  positionable  against 
Ge(Li)  or  intrinsic  germanium  cryostat  (flat  face)  . 

2 . 4 NBS  11-Line  Gamma-Ray  Standard 

2.4.1  An  NBS  11-line  gamma-rav  standard,  less  than  one 
year  old. 

3 . DETECTOR  CALIBRATION 

3.1  The  detector  must  be  calibrated  to  give  the  correct 
efficiency  for  gamma  rays  emitted  from  a 1 . 27-cm-diameter 
source,  using  a smaller  diameter  source  as  a standard.  This 
is  accomplished  by  measuring  the  count  rate  under  each  peak 

of  the  standard,  first  near  the  face  of  the  detector  cryostat, 
and  then  10  cm  farther  back. 

3.2  This  measurement  is  made  at  both  places  with  the 
11-line  NBS  gamma-ray  source.  This  gives  the  efficiencv 
e(x)  at  some  unknown  distance  x from  the  effective  center  of 
the  detector,  and  e(10  + x)  at  10  cm  + x.  (At  10  cm,  the 

1 . 27-cm-diameter  foil  size  leads  to  a correction  of  <2  per- 
cent, as  compared  to  the  source  which  is  ^0.5  cm  diameter.) 
Plot  the  efficiency  at  each  position  on  log-log  plotting 
paper  and  read  off  the  efficiency  at  each  position  for  an 
energy  E^ , corresponding  to  one  of  the  hottest  foils,  sav 
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Au(n, y) °Au  (412  keV) , and  calculate  the  ratio 
Rcal  * e(x)/c(10+x).  Repeat  the  measurement  at  both  posi- 
tions with  the  gold  foil  and  get  Multiply  the  e(x) 

by  Rgo^^/Rcai  to  get  the  correct  efficiency  for  a 1.27-cm- 
diameter  foil  at  x.  Repeat  the  procedure  for  one  or  two 
higher  gamma-rav  energies  with  another  hot  foil  (say  ^Mn 
with  gamma-ray  lines  at  847  and  1881  keV) . Now  slide  the 
e(x)  curve,  as  measured  with  the  11-line  standard,  to  pass 
through  these  three  corrected  points.  This  represents  the 
proper  efficiency  curve  for  a 1 . 27-cm-diameter  foil  at  the 
distance  x from  detector  "center." 

3.3  From  count-rate  measurements  at  x and  10+x,  it  was 
found  that  the  effective  center  of  a particular  13  percent 
Ge(Li)  gamma-ray  detector  (1)  is  2.1  cm  from  the  front  face 
of  the  cryostat.  Thus,  the  center  of  the  standard  source 
and  of  all  the  foils  (the  center  being  at  half  the  foil 
thickness)  must  be  located  to  within  0.2  mm  to  keep  the 
positioning  uncertainty  down  to  27„. 

3.4  The  above  same  calibration  must  be  repeated  for 
■737 

the  “ Np  source  behind  the  1.27-cm  lead  shield. 

3.5  A simpler  approach  can  be  used  if  the  true  axis  of 
the  detector  is  located.  A careful  measurement  and  x-y  plot 
of  count  rate  versus  position  at  the  cryostat  face  will 
accomplish  this.  The  point  source  is  then  simply  moved  off 
axis,  0.7  of  the  way  to  the  periphery  of  the  1 . 27-cm-diameter 
foil.  The  efficiency  measured  here  is  approximately  the  area- 
weighted  efficiency  over  the  1 . 27 -cm-diameter  foil. 

4.  COUNTING  PROCEDURE  AND  AREA  ANALYSIS 

4.1  Set  the  amplifier  gain  so  that  2 MeV  corresponds  to 
4000  channels  on  the  analyzer,  and  the  1-MeV  point  at  about 
2000  channels  so  that  the  zero-channel  offset  is  small.  Set 
the  precision  pulser  so  that  it  falls  at  %1.9  to  1.95  MeV, 
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well  above  the  ^Mn  1811-keV  line.  Place  the  foil  in  the 
counting  position.  Run  the  analyzer  on  clock  time  with  the 
pulser  running.  The  ratio  of  the  pulser  events  appearing  in 
the  peak  at  vl.9  to  1.95  MeV  to  the  pulses  generated  gives  the 
correction  for  the  combination  of  analyzer  deadtime  and  pulse- 
pileup  losses  from  the  peak.  The  analyzer  can  be  run  on  live 
time  instead,  but  in  this  case  the  true  clock  time  must  be 
measured.  Alternatively,  the  pulser  events  can  be  counted 
with  a scaler,  care  being  taken  that  the  analyzer  and  scaler 
are  turned  on  and  off  together.  , 

4.2  The  peak  analysis  routine  must  be  the  same  for 
pulses  from  the  pulser  as  for  those  from  both  calibration 
and  foil  events.  In  one  method,  the  counts  are  plotted  for 
the  peak  and  the  nearby  region  (v5  keV  or  more  on  each  side) . 

By  fitting  a straight  line  through  the  baseline,  the  base- 
line area  can  be  subtracted  from  the  peak  area.  In  counting 
fission  foils,  the  peak  shape  must  be  examined  carefully  for 
the  presence  of  a very  close  neighboring  peak,  in  which  case 
peak-shape  analysis  must  be  used.  The  counting  statistics 
must  be  very  good  for  the  peak -shape  analysis  to  be  accurate. 

The  analysis  can  either  be  done  by  hand,  or  with  the  SAMPO 
code  (2).  The  analyzer  should  be  run,  whenever  possible, 
until  the  peak  area  is  10,000  counts  or  greater. 

5 . BACKGROUNDS 

5 . 1 Room  Background 

5.1.1  This  must  be  kept  to  a minimum  by  selecting  a 

low-background  counting  area,  stacking  a lead  shield  of  at 

least  5-cm  thickness  around  the  detector,  and  moving  away 

all  sources  not  being  counted.  A room-background  run  is 

taken  and  is  definitely  required  for  long  (overnight)  counts, 

90  89 

such  as  for  the  Zr(n,2n)  Zr  reaction,  which  has  a threshold 
energy  Et  at  about  14  MeV,  and  consequently,  a low  specific 
activity,  or  for  ^Ni(n,p)^Co,  where  the  half-life  is  very 
long. 
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5 . 2 Fission-Foil  Background 

5.2.1  Fission-foil  backgrounds  must  be  carefully 
measured  because  the  foils  are  radioactive  with  a very  com- 
plex gamma-ray  emission  spectrum,  and  because  they  are  re- 
used due  to  their  high  replacement  cost.  If  they  had  been 
irradiated  within  several  half-lives  before  the  next  planned 
irradiation  (of  at  least  the  same  or  higher  fluence) , the 
background  should  be  measured  more  than  once  to  separate  out 
the  normal  line  {537,  743,  or  1593  keV  (see  Table  I)},  with 
relatively  short  half-life,  from  any  long-lived  radioactivity 
contribution  due  to  natural  radioactivity  (or  to  other  fission 
fragments) . The  background  peak  area  that  will  still  be  pres- 
ent during  the  counting  time  after  the  next  irradiation  is 
then  calculated  by  adding  the  time-dependent  area  of  the 
gamma-ray  line  with  known  half-life  to  that  of  the  steady- 
state  component. 

6.  DATA  ANALYSIS 

6.1  Correct  the  peak  areas  for  analyzer  deadtime  losses 
and  pulse-pileup  losses  by  multiplying  the  area  by  the  ratio 
of  the  number  of  precision  pulses  generated  during  the  count- 
ing period  to  the  number  found  in  the  pulser  peak  at  *vl.9  MeV. 
Divide  by  detector  efficiency,  gamma-ray  intensity  (gamma/ 
reaction,  Table  I) , and  fission  yield  (for  fission  foils) . 
Correct  for  gamma-ray  self  absorption  by  using  the  approxi- 
mate expression 

+Et( M/2) 

0O  - <J>e  (1) 

where  $ is  the  measured  activation  decay  rate  at  the  time  of 

counting,  the  corrected  activation,  Z^  the  macroscopic 

2 ^ 

gamma-ray  total  cross  section  (cm  /g) , and  M the  thickness 
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(g/cm  ) of  the  foil  and  encapsulation  material.  Equation  (1) 
is  accurate  to  ~ 1 percent  in  0Q  for  0o/0  <1.2.  The 

number  N of  atoms  made  radioactive  (or  in  the  case  of  fission 
foils,  the  number  of  fissions  induced)  by  the  neutron  irradia- 
tion is  then  given  by: 


«A  >t(i)  elt(wl  [l-e-  ‘tlclj-1 


(2) 


Vf 


where  is  the  net  peak  area  corrected  for  detector  efficiency, 
t(i)  the  neutron  irradiation  time,  t(w)  the  wait  time  between  the 
end  of  the  irradiation  interval  and  the  beginning  of  the  counting 
interval,  and  t(c)  is  the  duration  of  the  counting  interval,  Y^ 
is  the  X/reaction,  and  Y£  is  the  fission  yieli. 

6.2  For  the  1593-keV  line  of  140La  (T.  - 40.23  hours), 

1 40  * / *• 

which  is  produced  Dy  the  Ba  decays  (Tjy2  m 12.8  days),  Eq.  2 
becomes 


r - X Bt(i)-j_1  r -XBt(e)  -i  _1 

»LXBt(i)^L-XB>  [X-e  J LX“e  j 


XLYf  * 0 


r -xBt(w)  -xLt(w)i 

.96  [c  - e j 


(3) 


where  NB,  X0,  and  are  the  number  of  fissions  produced,  the  decay 
constant  for  X40Ba,  and  the  fission  yield  for  X4oBa  (see  Table  I), 
respectively,  N.  is  the  peak  area  divided  by  detector  efficiency  for 
the  1593-KeV  14®La  gamma  rays,  XL  the  decay  constant  for  140La, 
and  0.96  is  the  (r/reaction)  for  the  1593-keV  gamma  ray  from  1(0^. 
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6.3  For  wait  times  of  10  days  or  more,  Eq.  3 reduces  to: 


V<w) 


NL\0t(i)  * (0.869)  C 


L J (4) 


Yf  * 0.96 


with  ~1  percent  error  at  10  days. 


6.4  The  value  of  N or  NQ  (Eqs.  2,  3,  or  4)  is  corrected 
for  any  significant  neutron  self  shielding  or  flux  depression  during 
irradiation,  and  R ■ N/NQ  (radioactive  atoms  per  target  atom)  is 
computed  from  the  measured  target  mass  and  isotopic  abundance  for 
the  reaction  (Table  1). 


7.  PRECISION 


7.1  The  precision  is  limited  by  the  counting  statistics, 
by  the  reproducibility  of  the  NBS  calibration  source  and  of  the 
location  of  the  foils  and  the  standard  source  with  respect  to  the 
detector,  and  by  the  reproducibility  of  the  peak-area-analysis 
routine  used  for  the  foil-countinq  and  background  measurement. 

7.2  A calibration  of  one  NBS  11-i.ine  source  mix  (88-,  122-, 
165- , 279-,  392-,  514-,  662-,  898-,  1173-,  1333-,  and  183b-keV  lines) 
against  another  resulted  in  a standard  deviation  of  1 percent  tor 
the  reproducibility  of  the  calibration  points.  The  source  location 
(including  source-si2e  correction)  leads  to  a precision  uncertainty 
of  ~ 3 percent  (if  a good  alignment  scheme  is  employed  for  both 
distance  and  lateral  alignment);  bcause  of  the  nearness  of  the 
foils  to  the  detector,  the  weaker  foils  must  be  counted  next  to 
the  face  of  the  aetector  cryostat.  Peak-area  analysis  ana  background 
subtraction  are  somewhat  inseparable.  If  ~1  percent  or  better 
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counting  statistics  can  be  achieved,  the  associated  pre- 
cision is  ^2  percent  for  non-fission  foils  and  2 to  5 per- 
cent for  fission  foils,  depending  on  the  line  chosen,  the 
irradiation  fluence,  and  the  magnitude  of  the  background. 

7.3  An  overall  precision  of  5 to  6 percent  is  sug- 
gested as  the  best  achievable  precision,  combining  all  the 

above  factors.  If  the  16.3-hour  *^Zr  line  of  743  keV  (4)  is 

23'1  '>37 

used  for  analysis  of  “Th  and  Np  fission,  the  error  may 
be  much  greater,  unless  shape  fitting  is  utilize!  and  measure- 
ments are  made  at  different  times  to  subtract  out  the  effect 
99 

of  the  Ho  gamma-ray  line  (67  hour) , which  has  nearly  the 
same  gamma-rav  energy  (740  keV) . 

8 . ACCURACY 

8.1  The  accuracy  of  counting  and  data  analysis  reflects 
directly  on  the  precision  (see  Sec.  7),  and  the  accuracy  of 
the  nuclear  data.  At  this  time,  it  is  felt  that  the  uncer- 
tainty in  the  branching  ratios  of  the  gamma-rav  decay  schemes, 
including  the  conversion  electron  coefficients,  is  ^5  per- 
cent for  most  cases,  and  may  be  as  large  as  10  percent  for 
some  individual  cases.  The  fission  yields  are  probably  known 
to  5 percent.  The  half-lives  are  known  to  1 percent  for  most 
cases,  and  will  not  be  important  except  for  waiting  times  of 
more  than  2 to  3 half-lives.  An  overall  accuracy  of  ^10  per- 
cent is  suggested  as  the  value  probably  achievable  for  measur- 
ing specific  activations,  after  some  experience  has  been 
achieved  with  the  source-location  and  area-analysis  schemes 
and  reliable  methods  have  been  developed;  a value  of  15  per- 
cent is  more  likely,  and  will  have  about  a 5 percent  effect 

on  the  accuracv  of  measuring  * (see  ASTM  Methods  E XX3  and 

eq 

E XX4) . 


.81 


REFERENCES 


(1)  Verbinski,  V.  V.,  Lurie,  N.  A.,  and  Rogers,  V.C., 
"Threshold-Foil  Measurements  of  Reactor  Spectra  for 
Radiation  Damage  Applications,"  Nuclear  Science  and 
Engineering  65,  316  (1978) . 

(2)  Routti,  J.  T.,  "SAMPO,  A Fortran  IV  Program  for  Com- 
puter Analysis  of  Gamma  Spectra  from  Ge(Li)  Detectors, 
and  Other  Spectra  with  Peaks,"  UCRL-19452  (1969). 

(3)  Wright,  H.  L. , Applied  Sciences  Division,  White  Sands 
Missile  Range,  New  Mexico  88002  (private  communication)'. 

(4)  The  value  of  743  keV,  given  in  the  "Chart  of  Nuclides" 
(ERDA,  Division  .of  Isotopes  Development,  1970)  agrees 
with  measurements  of  Ref.  1.  The  value  of  747  keV, 
listed  in  the  "Table  of  Isotopes"  does  not,  (Lederer, 

C.  M. , Hollander,  J.  M. , and  Perlman,  I.,  John  Wiley  & 
Sons,  Inc.,  New  York,  1967). 

(5)  Meek,  M.  E.,  and  Rider,  B.  F.,  "Compilation  of  Fission 
Product  Yields,"  NEDO-12154-1,  Vallecitos  Nuclear 
Center,  Pleasanton,  California  94566 


TABLE  I 


Activation  Foils  (1.27-cm  diam) 


Reaction 

E, 

(MeV) 

Ey 

(keV) 

Gamma/  Reaction 
(Fission  Yield.  %> 

7\/i 

Foil  Mass 

«8) 

Isotopic 

Abunoanc 

(%» 

'"AuOi.yt^Au 

0 

412 

0.95 

2.696  days 

0.056 

100 

**Co  (e.  y)*°Co 

0 

83£1 

1.00 

5.258  yr 

0.057 

100 

O 

1.00 

**Mn  (n.  y)**Mn 

0 

847 

0.99 

2.58  h 

0.030 

100 

1811 

0.29 

mtl  (a./)“*La 

0(0.01)4 

1596 

0.96  (6.29) 

40.23  h 

0.281 

100 

"PU  (*./)“°La 

0(0.01)c 

1596 

0.96  (5.24) 

40.23  b 

0.150 

100 

"Npie./j'^La 

0.5' 

1596 

0.96  (5.69) 

40.23  h 

0.580 

100 

mXp(*.nmZr 

0.5 

743 

0.92  (6.01) 

16.8  h 

0.580 

100 

"*Ini/i,s')lll“ln 

t.O 

335 

0.50 

4.50  h 

0.255 

95.7 

**U<n./>,-La 

1.45 

1593 

0.96  (6.02) 

40.23  h 

0.495 

luo 

a*Th(n./)‘*’Ba 

1.75 

537 

0.256  (7.91) 

12.8  days 

1.066 

100 

“Th|i«./),TZr 

1.75 

743 

0.92  (4.12) 

16.8  h 

1.066 

100 

,‘Fe(e.0)MMn 

2.20 

835 

1.00 

303  days 

0.142 

5.82 

**Rl  (et0)**Co 

810 

0.99 

71.3  days 

0.282 

67.3 

"S(>i.P)3,P 

Betas 

1.000 

14.3  days 

4.10 

95.0 

“Vg^.pi^Na 

6.3 

1369 

1.00 

15.0  h 

0.030 

79 

**F*  (e,p)‘*Mn 

7.5 

847 

0.99 

2.58  h 

0.142 

91.7 

1811 

0.29 

^Al  (e,  a)MNa 

8.7 

1369 

< 1.00 

15.0  h 

0.257 

100 

l*I (",  2ji)13*I 

11.0 

386 

0.34 

13.05  days 

0.657 

100 

667 

0.33 

“Zrt*.  Isi^Zr 

14 

910 

• 0.99 

78.4  h 

1 

0.108 

51.5 

Note* 


*.(.« 


».8.c 


‘Cadmium  cover  0.05  to  0.X0  cm  Stick. 

bUse  ’*Co  instead  of  ir,Au  and  s>Ma  for  very  long  irradiations. 

*£,_»  O.Ot  MeV  with  “8  sphere  (impoctnnt  for  soft  (TRICAR  spectra.  where  S(£)  < 0.01  will  otherwise  dominaiel. 

When 


U or  a(Pu  foil  Is  1#8-covered,  also  cover  a U and  " Np  foils  so  that  accurate  corrections  can  be  made  for 
U and  a,Pu  impurities  m these  taigta  E,  (oils. 


Jt0.23-b  daughter  of  12. 80-day  l“3a.  Wait  five  days  for  maximum  decay  rate  iASTM  E 393). 

*F(ssioa  yields  are  for  bombardment  with  fission-spectrum  neutrons.  For  thermal  and  14  MeV,  see  Ref.  5. 

(Use  "Zrfor  low  nuance  (3  x I0uto  3 x lO^n/cm1).  Use  peak  shape  analysis  or  measure  twice,  seven  days  apart, 
to  strip  off  740-keV  “Mo  gamma  ray  ir,„  • 67  h). 

^’Th  radioactivity  interferes  with  ‘“La  lice. 

''Requires  separate  detector,  and  calibration  tecnnique  is  complex. 

' Maximum  manganese  Impurity  • 0.001%,  cadmium -cove  red.  Omit  ’,Fe(*.ptJ*Mn  for  long  irradiations. 
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DESIGNATION:  E XX 3 


Standard  Method  for 

UNFOLDING  NEUTRON  SPECTRA 
FOR  RADIATION  HARDNESS  TESTING 


1.  SCOPE 

1.1  This  method  describes  a set  of  standard  procedures 
for  unfolding  a neutron  spectrum  +>(E)  with  the  SAND  II  code, 
using  as  input  data  a set  of  threshold- foil  activations,  with 
threshold  energies  Et  effectively  ranging  from  cadmium  cutoff 
(••-  3 x 10" 7 MeV)  to  vl4  MeV . 

1.2  The  selection  and  exposure  of  a set  of  activation 
foils  is  covered  in  ASTM  Method  E XXI,  "Method  of  Irradiating 

a Standard  Set  of  Neutron  Threshold  Activation  Foils  for  Radia- 
tion Hardness  Testing."  In  addition,  a minimum  set  of  reliable 
foils  is  presented  in  that  method. 

1.3  The  measurement  of  the  specific  activities  of  the 
threshold-activation  foils  is  covered  in  ASTM  E XX2 , "Measur- 
ing Foil  Activities  for  Radiation  Hardness  lesting." 

1.4  The  neutron  spectrum  $ (E)  as  measured  by  this 
method  may  have  a wide  variety  of  uses.  Nevertheless,  this 
method  is  specifically  addressed  to  characterizing  $(E)  for 
the  purpose  of  comparing  spectra  used  in  radiation  hardness 
testing.  It  is  also  possible  to  express  $(E)  in  terms  of 

a single  parameter  4»ga  that  gives  a measure  of  its  effective- 
ness in  producing  radiation  damage  in  silicon  semiconductor 
devices.  This  is  discussed  in  ASTM  Method  E XX4 , 
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'Characterizing  Neutron  Spectra  in  Terms  of  1-MeV  Equivalent 
Fluence  for  Radiation  Damage  in  Silicon." 

1.5  Periodic  updating  of  nuclear  data  will  improve 
the  accuracy  of  the  threshold-foil  activation  cross  sec- 
tions, branching  ratios,  and  decay  constants,  as  well  as 
calculations  of  the  damage  function  D(E) . These  must  be 
utilized  to  update  the  input  data  for  this  method  from 
time  to  time  to  improve  its  accuracy,  which  will  require 
updating  these  standard  methods. 

1.6  All  E = 0 foils  (both  1/v  and  fission  foils)  are 
10  c 

cadmium  or  B covered,  which  excludes  thermal -neutron  measure- 
ments. This  facilitates  a more  reliable  solution  at  higher 
energies,  where  practically  all  of  the  radiation  damage  occurs 
in  silicon. 

1.7  Although  codes  other  than  SAND  II,  such  as 
SPECTRA,  CRYSTAL  BALL,  and  OPTIMO  (also  available  at  Radia- 
tion Shielding  Information  Center  (RSIC)  (1)},  can  give 
results,  the  SAND  II  code  was  chosen  because  of  the  built- 
in  feature  of  retarding  the  formation  of  spurious  structure 
(see  Sec.  11.1).  Also,  it  provides  in  convenient  format 
the  Aq's  needed  to  evaluate  the  agreement  between  the  foil 
data  and  the  trial  spectrum.  This  plays  an  important  part 
in  rejecting  erratic  foil  data  and  in  choosing  a more 
realistic  trial  spectrum;  a user  interaction  that  is 
important  to  unfolding  accuracy. 

2 SIGNIFICANCE 

2.1  This  method  provides  the  neutron  fluence  $(E) 
for  general  use,  such  as  for  calculating  radiation  damage 
to  materials  other  than  silicon,  where  the  damage  versus 
neutron  energy,  D(E)  is  known  or  is  calculable  from  known 
neutron  cross  sections. 
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2.2  This  method  and  the  supporting  methods  (Sec.  1) 
constitute  a set  of  draft  standards  for  characterizing 
diverse  neutron  fields.  This  characterization  is  of  interest 
in  planning  irradiation  schedules  for  parts  screening  or  for 
sample-specification  tests,  and  in  comparing  radiation  damage 
studies  carried  out  in  different  neutron  fields. 

3 . DEFINITIONS 

3.1  The  description  of  terms  relating  to  dosimetry 
are  found  in  ASTM  Definitions  E 170,  "Terms  Relating  to 
Dosimetry . " 

3.2  As  is  the  saturated  activity  of  the  sample  (in 

units  of  reciprocal  time),  as  defined  in  ASTM  Method  E 261-70, 

"Standard  Method  for  Measuring  Neutron  Flux  by  Radioactivation 

Techniques."  R is  the  measured  number  of  activations  per 

nucleus  produced  by  irradiating  Nq  target  nuclei  (of  the 

isotope  that  can  produce  the  reaction)  for  an  irradiation  time 

t . . For  constant  reactor  power,  R is  defined  as  the  measured 

ratio  Rm  * A t./N  . 
m s l o 

4.  SUMMARY  OF  THE  METHOD 

4.1  A set  of  foils  is  selected  with  thresholds  Et  that 
vary  from  0 to  14  MeV,  the  foils  are  weighed,  covered  with 
cadmium  or  as  necessary,  and  exposed  to  a neutron  field 

having  a known  spectral  distribution  cp (E)  . The  resultant 
radioactivity  is  measured  to  determine  the  number  N of  radio- 
active atoms  produced  during  the  irradiation  of  Nq  target  atoms 
of  the  chosen  isotope.  For  each  foil  species  x, 


N 


N 


o 


ax(E)<f>(E)  dE, 


(1) 


where  <t>x(E)  is  the  cross  section  for  the  reaction  being 


86 


measured. 


by  the  ratio  Rffl(x) 


The  set  of  measured  specific  activities,  given 
Rm(x)  * (N/N0)x.  *-s  used  as  input  to  the  SAND  II 


unfolding  code,  along  with  a trial  spectrum  4>tr(E)  that  con- 
tains all  the  physical  features.  These  features  may  include 
a "dip"  due  to  a resonance,  a peak,  or  a fission  spectrum 
with  1/E  moderated  component  (tail) . 

4.1.1  Using  Eq . 1 with  4>tr(E)  in  place  of  <t> (E>  » the 
SAND  II  code  calculates  the  various  Rtr(x),  utilizing  the 
SAND  II  input  library  of  cross  sections  °X(E) • In  calculat- 
ing the  various  Rtr(x)  . the  input  spectrum  4>Cr(E)  is  normal- 
ized so  as  to  minimize  the  standard  deviation  S associated 

o 

with  the  "zeroth  iteration"  differences  A (x)  - (R  -R_..)/R„ 

o m tr  m 

between  the  measured  and  calculated  specific  activities. 

4.1.2  For  foils  where  Aq(x)  is  largest,  the  SAND  II 
code  adjusts  the  spectrum  in  the  energy  interval  above  E£(x) , 
so  as  to  produce  a first-iteration  solution,  4>^(E)  , with 
reduced  A^(x)  and,  consequently,  with  reduced  S^.  These 
iterations  are  continued  until  the  n-th  iteration,  where 

Sn  <5  percent. 

4.2  SAND  II  Limitations 

4.2.1  The  threshold-foil  spectrometry  method  is  not 
truly  a spectrometry  technique,  but  rather  a mild  perturba- 
tion technique.  As  such,  it  provides  a reliable  output  only 

with  considerable  constraints  applied,  especially  for  the 

237 

region  below  1 MeV,  where  only  one  true  threshold  ( Np, 

0.5  MeV)  exists.  This  constraint  is  that  for  reactor 
spectra,  the  epicadmium  solution  resembles  a GODIVA  spec- 
trum with  1/E  tail  fitted  (Sec.  6.1.2).  Since  this  provides 
only  one  adjustable  parameter,  that  of  varying  the  1/E  fit- 
ting point  until  the  solution  resembles  the  trial  spectrum, 
it  is  possible  for  this  limiting  case  of  reactor  spectra  to 
obtain  a reliable  solution  with  one  real  threshold,  an 
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21S  2T<J 

artificial  one  at  10  keV  (*  3U  or  Pu  in  boron) , and  a 
1/v  detector  providing  information  below  1 MeV.  This  is 
carried  out  with  some  interaction  with  the  code. 

5.  INTERACTION  WITH  SAND  II  CODE 

5.1  Examination  of  Rm(x)  and  $tr(E) 


5.1.1  If  the  solution  4>n(E)  Is  physically  meaningful, 

it  will  have  the  same  characteristics  as  $cr>  This  happens 
when  no  erroneous  values  of  Rm  are  input  to  SAND  II,  and  when 
the  starting  spectrum  $tr(E)  contains  all  the  physical  char- 
acteristics of  the  true  spectrum.  It  is  the  task  of  the  user 
to  examine  the  set  of  AQ(x)  to  either  reject  single  spurious 
values  of  the  corresponding  Rm(x)  . or  to  adjust  $tr(E)  In 
regions  above  the  threshold  energy  Ec(x)  where  two  or  more 
corresponding  values  of  are  large  and  agree  in  sign 

and  magnitude. 

5.2  Choosing  Cr for  Unfolding  Reactor  Spectra 

5.2.1  Most  reactor  spectra  used  for  radiation  damage 
studies  can  be  represented  by  a fission  spectrum  with  properly 
normalized  l/E  slowing-down  component.  The  trial  spectrum 
can  be  represented  adequately  in  this  case  by  a fission  spec- 
trum with  a fitted  l/E  tail.  A trial  fitting  is  first  used, 

197 

and  if  the  Aq(x)  for  the  cadmium-covered  1/v  detector  ( Au, 
59Co,  or  55Mn)  and  the  flat  detector  (239Pu  or  i35U  boron- 
covered  fission  detectors)  are  large  and  positive  (or  negative) , 
then  the  l/E  tail  is  normalized  at  a higher  (or  lower)  energy 
to  the  fission  spectru,  and  SAND  II  rerun.  Prescribed  values 
are  given  in  Sec.  6.1  for  the  l/E  normalizing  point  for  three 
different  fission  spectra  commonly  used  in  radiation  damage 
work . 
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5.2.2  If  0tr(E)  is  too  low  in  the  energy  region  above 
a few  MeV,  for  a reactor  spectrum,  the  Aq(x)  will  be  positive 
and  may  increase  with  Et  of  the  foils.  For  this  case,  there 
are  enough  foils  used  in  this  prescription  for  SAND  II  to 
achieve  a good  solution  in  only  a few  iterations,  so  that 
<p cr (E)  need  not  be  adjusted  here.  A more  common  problem  here 
is  that  one  of  the  AQ(x)  is  erroneous,  requiring  many  itera- 
tions to  achieve  a solution.  This  is  seen  as  a physically 
unreal  oscillation  in  the  solution  <|>n(E),  in  which  case  the 
spurious  activation  Rm(x)  must  be  removed  and  SAND  II  rerun 
without  this  activation. 


5 . 3 Recognition  of  a Well-Behaved  Solution 

5.3.1  When  the  Rm(x)  and  4>tr(E)  are  self  consistent, 
the  SAND  II  perturbation  code  will  converge  in  a few  iterations 
(n  = 1 to  10),  and  4>n(E)  will  be  much  like  <J>tr(E)  . 

5 3.2  The  prescriptions  set  forth  in  Section  6 are 
designed  to  achieve  such  a solution.  It  is  most  easily 
recognized  by  comparing  plots  of  <f>tr(E)  and  $n(E)  on  log- 
log  paper,  but  with  E*<J>(E)  plotted  instead  of  <t>(E).  In  this 
way,  the  1/E  slowing-down  region  appears  as  a flat  line,  and 
the  slope  of  4>(E)  above  a few  MeV  is  not  nearly  as  steep. 
Deviations  of  d>n (E)  from  <J>tr(E)  are  therefore  much  more 
apparent  in  the  E*$(E)  plots. 

6.  SAND  II  OPERATING  PROCEDURE 


SAND  II  Inputs 


6.1.1  The  SAND  II  code  is  operated  in  the  TIME 
INTEGRATED  (time  independent)  mode,  which  standardizes  the 


activation-foil  input  format  for  both  fast-burst  and  steady- 
state  irradiations.  The  inputs  are  the  trial  spectrum 
the  specific  activities  Rm(x)  of  the  foils,  and  the  foil- 
cover  data. 


6.1.2  Trial  Spectrum  $tr 

6. 1.2.1  For  a GODIVA-type  reactor,  select  spectrum 
No.  5 of  the  SAND  II  library  of  trial  spectra  for  a glory- 
hole  spectrum. 

6. 1.2. 2 For  a spectrum  at  50  cm  from  a GODIVA-type 

reactor  ^1.5  meters  above  a concrete  floor,  use  trial  spec- 

„ o 

trum  No.  5 with  a 1/E  component  fitted  at  10  MeV  {this 

will  be  compatible  with  the  1/v  foil  (^Co,  ^Mn,  or  ^^Au) 

235  239 

and  the  boron-covered  U or  Pu  low-energy  fission  foil, 
and  will  avoid  distortions  in  the  solution  $(E)  for  E >0.01 
MeV  that  result  from  omitting  the  1/E  component) . 


6. 1.2. 3 
at  0.15  MeV. 


For  a TRIGA  spectrum,  fit  the  1/E  component 


6. 1.2. 4 To  obtain  the  trial  spectrum,  first  run 

SAND  II  with  the  trial  spectrum  No.  5 option  called  for.  The 
SAND  II  output  will  include  a printout  of  start  spectrum  No.  5 
(620-point  spectrum  from  10"^  to  18  MeV).  Normalize  the  1/E 
spectrum  at  the  proper  point,  as  prescribed  above  and  input 
the  resulting  hybrid  spectrum  at  AE/E  - 10  percent  energy 
intervals  above  the  normalizing  point  energy  (i.e. .,  at  0.15, 
0.015,  0.0015  MeV,  etc.,  down  to  ■vlO"10  MeV).  The  hybrid 
spectrum  is  input  as  SPECTRUM  TABULAR,  as  described  in  the 
input  instructions. 

6. 1.2.5  For  cases  other  than  the  three  mentioned  above, 
choose  a normalizing  point  at  0.03  MeV  for  the  1/E  tail,  and 
proceed  with  the  fitting  operation  as  described  in  Sec.  5.2. 

6.1.3  Threshold-Foil  Data 

6. 1.3.1  The  threshold-foil  isotope,  the  type  of  reac- 
tion, the  specific  activity  R , and  the  foil-cover  data  are 
required  inputs  to  the  SAND  II  code.  For  example, 

U235F  BORON  0.101 


and 


AU197G  CADMIUM  0.00480 

235  235 

correspond  :o  the  U(n,f)  reaction  with  the  U foil 

covered  with  1.68  g/cm2  of  10B,  and  the  197Au(n , y) 198 

197 

reaction  with  the  Au  foil  covered  with  a 0 . 1-ctn- thick 
cadmium  foil.  A specific  activation  of  5.9-12  corresponds 
to  Rm(  JU)  • 5.9  x 10  fissions  produced  per  *"J3U  atom 
present  in  the  target.  According  to  the  structure  of  the 
code,  whereas  is  input  for  fission  foils,  for  non-fission 
foils  Rm  is  multiplied  by  X «■  0.693/T^  of  the  radioactive 
species  for  input  to  SAND  II.  Three  or  more  1/v  detector 
foils  such  as  97Au,  ^9Co,  and  "^Mn  can  be  used  with  a SAND  II 
trial  run,  but  the  one  with  the  value  of  AQ  most  nearly  the 
average  of  the  three  should  be  chosen  for  the  final  run  of 
SAND  II,  the  others  being  rejected  to  keep  the  number  of  itera- 
tions of  SAND  II  (and  therefore  the  spurious  structure)  to  a 
minimum . 

6.2  Acceptance  Criteria  for  a Good  Solution  $n(E) 

6.2.1  If  an  E*$(E)  plot  of  $n(E)  shows  the  same  general 
shape  as  a similar  type  plot  of  $tr(E),  with  4>Cr(E)  selected 
as  in  Sec.  5.1.1,  then  4>n(E)  is  likely  to  represent  a good 
solution.  This  is  usually  accompanied  by  n being  small,  as 
mentioned  in  Sec.  5.4.  If  <j>n (E)  exhibits  a shape  very  much 
unlike  $tr,  the  user  must  examine  the  Aq(x)  (given  by  the 
SAND  II  printout)  for  spurious  values  of  Aq(x) , and,  there- 
fore, Rm(x) . Any  spurious  value  is  rejected  and  the  SAND  II 
code  is  rerun  without  this  foil.  On  the  other  hand,  when  more 
than  one  foil’  in  a given  region  of  threshold  energies  E_(x) 
shows  a large  positive  value  of  .\o(x)  and  these  values  agree 
in  sign,  the  $cr(E)  must  be  increased  in  that  region  for 
positive  A (x)  . and  vice  versa.  In  some  cases,  such  as  at 
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high  threshold  energies  where  $tr(E)  obtained  from  a cal- 
culation has  often  been  found  to  be  too  low,  the  SAND  II 
code  may  adjust  the  spectrum  in  very  few  iterations.  This 
happens  when  the  Aq(x)  are  progressively  larger  in  mag- 
nitude with  increasing  Ec(x)  of  the  foil,  and  the  activations 
Rm(x)  are  consistent  with  a smooth,  nonoscillating  solution 
^(E)  in  this  region.  In  this  case,  the  SAND  II  solution  is 
valid  {the  code  need  not  be  rerun  with  $tr(E)  adjusted}. 

7.  MINIMUM  FOIL  SET 

7.1  Table  I lists  foils  that  have  been  successfully 
used.  A 1/v  foil  is  required  for  the  low-energy  region. 

IQ  7 

Au,  cadmium  covered,  is  suggested  with  appropriate  epi- 

235 

cadmium  self-shielding  correction.  Use  U,  or  possibly 
239  10 

*Pu,  in  UB  as  the  "Et  ■ 10  keV"  fission  foil,  along  with 
237Np  (Et  ■ 0.5  MeV)  , 238u  (corrected  for  238U  impurity  with 
the  235U-foil  data),  115In,  54Fe(n.p)54Mn,  60Ni(n,p)60Co, 
24Mg(n,p)24Na,  27Al(n,a)24Na,  and  90Zn(n, 2n)892r . 

8.  AVAILABILITY  OF  THE  SAND  II  CODE  AND 
CROSS-SECTION  LIBRARY 

8.1  The  SAND  II  program  tape,  including  a large  catalog 
of  trial  spectra,  is  available  from  RSIC.  Documentation  of  the 
code  is  also  available  at  RSIC  (1) . 

9.  THE  SAND  II  CROSS-SECTION  LIBRARY 

9.1  The  present  cross-section  library  tape  contains  an 
evaluated  self-consistent  set  of  cross-section  data.  This  is 
referred  to  as  the  1974  evaluated  cross-section  library  tape, 
and  constitutes  part  of  this  standard.  When  this  tape,  or 
any  other  foil-activation  data,  becomes  updated,  it  is  under- 
stood that  these  ASTM  standard  methods  will  similarly  require 
updating. 


92 


10.  PRECISION 


10.1  Using  Che  prescriptions  outlined  above  and  in 
the  two  accompanying  ASTM  Methods  (E  XXI  and  E XX2 , 
mentioned  in  Sec.  1),  both  a glory-hole  spectrum  and  a spec- 
trum at  SO  cm  from  a GODIVA-cype  reactor  were  measured  (2,3) 
by  two  different  experimental  groups  in  an  effort  at  evaluat- 
ing the  reproducibility  of  the  method.  Using  different  foil 

I 

packets,  but  the  same  types  of  theshold  foils,  both  spectra 
agreed  within  less  than  2 percent  in  terms  of  the  hardness 
parameter,  A _/$  (see  ASTM  Method  E XX4)  . Different  Ge(Li) 
gamma-ray  detectors  were  used,  each  being  cross  calibrated 
with  NBS  standard  gamma-ray  sources  having  a quoted  absolute 
accuracy  of  2 percent.  A standard  consisting  of  an  11-line 
source  mixture  was  used  at  both  laboratories,  indicating  the 
precision  of  NBS  standard-source  intensities  is  1 percent  for 
the  11  lines.  Assuming  the  prescriptions  in  these  ASTM 
standard  methods  are  carefully  followed,  a precision  of  10 
percent  should  be  easily  achievable  in  *>eq/$,  and  5 percent 
highly  probably  with  additional  work  (such  as  checking 
activation  ratios  against  those  of  similar  spectra  found  in 
the  literature  (2)  to  help  diagnose  erratic  foil  activations). 

11.  ACCURACY 

11.1  The  accuracy  is  specified  here  in  terms  of  the 

1-MeV  equivalent  flux,  $ , and  the  hardness  parameter, 

/«$,  both  of  which  are  defined  in  ASTM  Method  E XX4.  It 
was  determined  in  a study  wherein  three  different  types  of 
reactor  spectra  were  Measured  (2) . These  spectra  were  found 
to  belong  to  a parametric  set.  wherein  they  can  be  approxi- 
mated by  a GODIVA  spectrum  with  1/E  tail  fitted  according  to 
the  degree  of  moderation  present  see  Sec.  6.1.1,  selecting 
$Cr(E)  . The  accuracy  evaluation  presented  here  is  limited 
to  spectra  belonging  to  this  set,  which  probably  covers  all 


reactor  spectra  used  for  silicon  radiation  damage  studies. 

11.2  Comparison  with  Known  Spectra 

11.2.1  The  SAND  II  results  for  $ /#  as  obtained  from 

eq 

♦ tr(E)  generated  according  to  Sec.  6.1.1,  agreed  with  $eq/$ 
as  obtained  from  calculated  values  of  $(E)  within  better  than 
3 percent.  The  calculations  were  verified  by  time-of-flight 
measurements  (2).  In  addition,  the  SAND  II  unfolding  was  car- 
ried out  using  the  calculated  <j>(E)  as  the  trial  spectrum,  $tr- 
In  this  case,  the  SAND  II  result  agreed  with  the  calculated 
spectrum  to  within  vl  percent  in  e(q / , indicating  that  the 
threshold-foil  data  were  consistent  with  the  calculation.  This 
and  other  data  show  chat  the  accuracy  is  very  sensitive  to  the 
choice  of  4>._(E),  especially  for  E <1  MeV,  where  only  one  use- 

cr  t -i  7 

ful  threshold  exists  Np(n.f))  above  10  keV. 

11.2.2  It  is  not  possible  to  check  the  accuracy  of  i>e(^ 
against  calculations  without  doing  more  exact  calculations, 
geometrically,  and  more  absolute  flux  monitoring. 

11.3  Sensitivity  to  Unfolding  Method 

11.3.1  A reactor  spectrum  was  unfolded  with  both  the 
SAND  II  and  the  SPECTRA  {also  available  at  RSIC  (1)}  codes. 

The  results  agreed  within  2 percent  in  eq / 4> . and  3 percent  in 

indicating  that  the  choice  of  codes  is  not  a major  factor. 
Most  of  the  difference  came  from  some  spurious  structure  in 
the  SPECTRA  solution  that  is  relatively  suppressed  in  the 
SAND  II  method. 

11.4  Variational  Studies  on  o (E)  and  R (x) 

x m 

11.4.1  The  SAND  II  unfolding  of  a GOD IVA- type  spectrum 
was  carried  out  by  varying  the  input  activations  of  individual 
foils,  R (x) , and  of  pairs  of  foils  (2).  This  simulates  a 
bodily  shift  (or  renormalization)  in  the  cross-section  curve. 
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ox(E) . The  results  showed  that  for  the  spectrum  at  50  cm 
from  a GODIVA-type  reactor,  which  was  165  cm  above  a concrete 
floor,  a variation  of  25  percent  in  the  activation  of  any  one 
foil  (out  of  nine)  with  threshold  below  2 MeV,  where  few 
thresholds  exist,  resulted  in  ''•2  percent  increase  of  ♦gq/#- 
For  foils  with  thresholds  between  2 and  6.5  MeV,  where  the 
neutron  flux  is  high  and  many  thresholds  exist,  a 25  percent 
change  in  ox(E)  of  any  single  foil  produced  a 1 percent  change 
in  ♦gq/^-  Above  6.5  MeV,  where  the  neutron  flux  is  very  low, 
the  corresponding  sensitivity  was  0.5  percent  for  one  foil. 

With  two  values  of  a (E)  varied  by  25  percent,  / <p  changed 

x eci  238 

by  as  much  as  4 percent  in  one  case  {+25  percent  for  U(n,f) 
and  -25  percent  for  2Sig(n,p)},  and  2 percent  in  another  (+25 
percent  for  both  “38U  and  ‘‘Slg)  • A total  of  nine  foils  was 
used  {58Ni(n.p),  24Mg(n,p).  127I(n,2n),  5eFe(n.p),  239Pu(n,f), 
238U(n , f) , 235U (n , f) , 197Au(n, y) , in  Cd.  and  237Np(n,f)}. 

11.4.2  The  variations  in  $ were  appreciably  greater 
than  the  corresponding  variations  in  $ /$.  They  were 

23Q  ect 

* +12  percent  for  the  Pu  foil,  2.7  percent  for  *°U, 
and  <0.5  percent  for  foils  with  thresholds  above  2 MeV.  Tne 
238U  and  24Mg  foils  produced  a 2 percent  increase  in  when 

both  were  raised  by  25  percent,  but  „ changed  only  0.3  per- 

238  e" 

cent  when  the  U activation  was  increased  25  percent  and  the 

24Mg  foil  decreased  by  the  same  amount.  The  large  sensitivity 

239  235 

to  changes  in  Pu  (or,  similarly,  U)  activation  (or  cross 
section)  arises  from  this  being  the  foil  with  accuracy  of 
2 percent.  A standard  consisting  of  an  11-line  source  mix- 
ture was  used  at  both  laboratories,  indicating  the  precision 
of  NBS  standard-source  intensities  is  1 percent  for  the  11  lines. 
Assuming  the  prescriptions  in  these  ASTM  standard  methods  are 
carefully  followed,  a precision  of  10  percent  should  be  easily 
achievable  in  $eq/$.  and  5 percent  highly  probably  with  addi- 
tional work  {3uch  as  checking  activation  ratios  against  those 


of  similar  spectra  found  in  the  literature  (2)  to  help 
diagnose  erratic  foil  activations). 

11.  ACCURACY 

11.1  The  accuracy  is  specified  here  in  terms  of  the 

1-MeV  equivalent  flux,  $ , and  the  hardness  parameter, 

both  of  which  are  defined  in  ASTM  Method  E XX4.  It 

eq 

was  determined  in  a study  wherein  three  different  types  of 
reactor  spectra  were  measured  (2) . These  spectra  were  found 
to  belong  to  a parametric  set,  wherein  they  can  be  approxi- 
mated  by  a GODIVA  spectrum  with  1/E  tail  fitted  according  to 
the  degree  of  moderation  present  {see  Sec.  6.1.1,  selecting 
The  accuracy  evaluation  presented  here  is  limited 
to  spectra  belonging  to  this  set,  which  probably  covers  all 
reactor  spectra  used  for  silicon  radiation  damage  studies. 

11.2  Comparison  with  Known  Spectra 

11.2.1  The  SAND  II  results  for  $ / $ as  obtained  from 

eq 


$tr(E)  generated  according  to  Sec.  6.1.1,  agreed  with  $eq/$ 
as  obtained  from  calculated  values  of  $(E)  within  better  than 
3 percent.  The  calculations  were  verified  by  time-of-flight 
measurements  (2).  In  addition,  the  SAND  II  unfolding  was 
carried  out  using  the  calculated  $(E)  as  the  trial  spectrum, 
$Cr.  In  this  case,  the  SAND  II  result  agreed  with  the  cal- 
culated spectrum  to  within  vl  percent  in  '*’eq/4>.  indicating 
that  the  threshold-foil  data  were  consistent  with  the  cal- 
culation. This  and  other  data  show  that  the  accuracy  is  very 

sensitive  to  the  choice  of  $ (E) , especially  for  E <1  MeV, 

cr  237 

where  only  one  useful  threshold  exists  { Np(n.f)}  above 
10  keV. 

11.2.2  It  is  not  possible  to  check  the  accuracy  of 
♦ against  calculations  without  doing  more  exact  calculations, 
geometrically,  and  more  absolute  flux  monitoring. 
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11.3  Sensitivity  to  Unfolding  Method 

11.3.1  A reactor  spectrum  was  unfolded  with  both  the 
SAND  II  and  the  SPECTRA  {also  available  at  RSIC  (1)>  codes. 

The  results  agreed  within  2 percent  in  $eq/$.  and  3 percent 
in  ♦gq.  indicating  that  the  choice  of  codes  is  not  a major 
factor.  Most  of  the  difference  came  from  some  spurious 
structure  in  the  SPECTRA  solution  that  is  relatively  sup- 
pressed in  the  SAND  II  method. 

11.4  Variational  Studies  on  a (E)  and  R (x) 

x m 

11.4.1  The  SAND  II  unfolding  of  a GODIVA-tvpe  spectrum 
was  carried  out  by  varying  the  input  activations  of  individual 
foils.  Rffl(x)  and  of  pairs  of  foils  (2).  This  simulates  a 
bodily  shift  (or  renormalization)  in  the  cross-section  curve, 
ox(E) . The  results  showed  that  for  the  spectrum  at  50  cm 
from  a GODIVA-type  reactor,  which  was  165  cm  above  a concrete 
floor,  a variation  of  25  percent  in  the  activation  of  any  one 
foil  (out  of  nine)  with  threshold  below  2 MeV,  where  few  thresh 
olds  exist,  resulted  in  ^2  percent  increase  of  $eq/$-  For  toils 
with  thresholds  between  2 and  6.5  MeV,  where  the  neutron  flux 
is  high  and  many  thresholds  exist,  a 25  percent  change  in 
ox(E)  of  any  single  foil  produced  a 1 percent  change  in  d> eq / . 
Above  6.5  MeV,  where  the  neutron  flux  is  very  low,  the  cor- 
responding sensitivity  was  0.5  percent  for  one  foil.  With  two 

values  of  a (E)  varied  bv  25  percent,  $„/$  changed  bv  as  much 

x 238 

as  4 percent  in  one  case  (+25  percent  for  U and  “ :ig)  . A 

total  of  nine  foils  was  used  (38Ni(n,p),  "Slg(n.p),  I(n,2n), 

56Fe(n,p),  239Pu(n,f),  238U(n,f),  235U(n,f),  197Au(n,  ).  in  Cd. 

and  “37Np(n , f) } . 

11.4.2  The  variations  in  were  appreciably  greater 

than  the  corresponding  variations  in  A /$.  Thev  were 

?->q  ecl  „ '•38 

A$gq  - +12  percent  for  the  " Pu  foil,  2.7  percent  for  “ U, 
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and  <05  percent  for  foils  with  thresholds  above  2 MeV . The 
*^33  24 

U and  - Mg  foils  produced  a 2 percent  increase  in  when 

both  were  raised  by  25  percent,  but  changed  only  0.3  per- 
238 

cent  when  the  U activation  was  increased  25  percent  and  the 

"*Mg  foil  decreased  by  the  sane  amount.  The  large  sensitivity 
•>39  235 

to  changes  in  “ ?u  (or,  similarly,  U)  activation  (or  cross 
section)  arises  from  this  being  the  foil  with  the  greatest 
sensitivity  to  neutrons  between  0.01-MeV  and  the  0.5-MeV 


threshold  of  ‘^Np. 


Fortunately,  the  cross-section  data  for 


039  ->35 

rission  foils  such  as  “ Pu  and  “ U are  accurately  known, 
and  a 25  percent  uncertainty  in  either  the  cross  section,  or, 
equivalently,  the  foil  count,  is  excessively  large.  Also,  the 
specific  activations  are  considerably  greater  than  for  the 
ocher  fission  foils,  and  are  thus  more  accurately  determined. 


11.5  Overall  Accuracy  Assessment 

11.5.1  The  probable  accuracy  in  q and  $ / $ , as 

r ■ eq  eq  ' 

estimated  from  the  three  contributing  factors  discussed 
above,  is  10  percent  tor  reactor  spectra  if  the  prescriptions 
outlined  in  these  methods  are  carefully  followed.  An  accuracy 
of  o to  7 percent  may  be  achievable  with  added  effort,  such 
as  comparing  activation  ratios  to  those  for  similar  spectra 
presented  in  the  literature  (2)  to  help  detect  erroneous 
values . 
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TABLE  I 


Activation  Foils 


(1.27-cm  diam) 


Reactico 

t, 

(MeV) 

Ey 

(keV) 

Gamma/  Reaction 
(Fission  Yield.  %) 

M/I 

. _ 

Foil  Maas 
<0 

Isotopic 

Abundance 

(%> 

Notes 

lflAul»."y)U*Au 

0 

412 

0.96 

2.696  days 

0.056  ’ 

too 

a 

"Co  (n,  >|SOCo 

0 

(173 

1.00 

5.258  yr 

0.057 

100 

».b 

1333 

1.00 

**Mn  («,  y)J*Mn 

0 

847 

0.99 

2.58  h 

0.030 

100 

a. 6 

1811 

0.29 

“Ufe./V’l.a 

0(0.  onc 

1596 

0.96  (6.29) 

40.23  h 

0.281 

100 

A.  4.* 

a,Pu(i»./)w#La 

0(0.01)c 

1596 

0l96  (5.24) 

40.23  h 

0.150 

100 

a,4.c 

wNp(»,/)l'°La 

0.5 

1596 

0.96  (5.69) 

40.23  h 

0.580 

100 

•,4.« 

aTNp(«.A),"Zr 

0.5 

743 

0.92  (6.01) 

16.8  h 

0.580 

100 

l“ln(n,s')‘l,"ln 

1.0 

335 

0.50 

4.50  h 

0.255 

95.7 

B*U  (e,/)H0La 

1.45 

1593 

0.96  (6.02) 

40.23  h 

0.495 

1U0 

Avd.* 

a,Th(s./)‘*Ba 

1.75 

537 

0.250  (7.91) 

12.8  days 

1.066 

100 

*•!•* 

a*Th(s,/)*7Zr 

1.75 

743 

0.92  (4.12) 

16.8  h 

1.066 

100 

A./,* 

MFe0i.A»MMn 

2.20 

835 

1.00 

303  days 

0.142 

5.82 

**NI  (*./>>J,Co 

2.9 

810 

0.99 

71.3  days 

0.282 

67.8 

5,S(>i./»5,P 

2.9 

Betas 

1.000 

14.3  days 

4.10 

95.0 

A 

*,Mg(".P)**Na 

6.3 

1369 

1.00 

15.0  h 

0.030 

79 

* 

**Fe  (*.£)wMn 

7.5 

847 

0.99 

2.58  h 

0.142 

91.7 

A.J 

1811 

0.29 

”Al(s,  a)MNa 

8.7 

1369 

1.00 

15.0  h 

0.257 

100 

A 

,aI(».2»)lBI 

11.0 

386 

0.34 

13.05  days 

0.657 

100 

667 

0.33 

"Zr(e,li|“Zr 

14 

910 

. 0.99 

78.4  h 



0.108 

51.5 

'Cadmium  cover  0.03  to  0.10  cm  thick. 

hUae  ”co  Instead  of  ‘"au  add  ’*Mn  for  very  long  irradiations. 

c E,  * 0.01  MeV  with  UB  sphere  [important  (or  so(t  (TRICA)  spectra,  where  «(£)  < 0.01  will  otherwise  dominatel. 
When  "U  or  B,Pu  foil  Is  l°B-covered.  also  cover  “U  and  “Np  (oils  so  that  accurate  corrections  can  be  made  for 
“u  tod  B,Pu  impurities  in  these  high  E,  (oils. 

4 40. 23 -h  daughter  o(  12.80-day  "%a.  Wait  five  days  (or  maximum  decay  rata  (ASTM  E 393). 

.'Fission  yields  are  (or  bombardment  with  (lssion-spectrum  neutrons.  For  thermal  and  14  MeV,  see  Ref.  S. 

'Use  "Zrfor  low  (luence  (3  * I0llto  3 x lO^n/cm1).  Use  peak  shape  analysis  or  measure  twice,  seven  days  apart, 
to  strip  off  T40-kaV  *Mo  gamma  ray  (7*,,,  • 67  h). 

*®*Th  radioactivity  interferes  with  ‘“La  line. 
h Requires  separata  detector,  and  calibration  technics  is  complex. 

'Maximum  manganese  Impurity  • 0.001%,  cadmium-covered.  Omit  J,Fe  1/i.prMn  for  long  irradiations. 


100 


DESIGNATION:  E XX4 


Standard  Method  for 

CHARACTERIZING  NEUTRON  SPECTRA  IN 
TERMS  OF  1-MeV  EQUIVALENT  FLUENCE 
FOR  RADIATION  DAMAGE  IN  SILICON 


1 . SCOPE 

1.1  This  method  describes  a standard  procedure  for 
characterizing  a neutron  field  with  spectrum  0(E)  in  terms 
of  the  <j>eq,  the  fluence  of  neutrons  at  ^1  MeV  required  to 
produce  the  same  radiation  damage  in  some  material  x where 
the  radiation  damage  function  Dx(E)  is  known. 

1.2  This  method  describes  a standard  procedure  for 
characterizing  the  shape  of  a neutron  spectrum  0(E)  in  terms 

of  the  hardness  parameter.  0e^/0.  It  is  defined  as  the  fluence 
of  neutrons  at  ''-l  MeV  required  to  produce  the  same  radiation 
damage  as  one  unit  of  fluence  of  neutrons  of  spectral  distri- 
bution 0(E)  . 

1.3  Although  these  standard  procedures  are  applicable 
to  characterizing  the  spectrum  in  terms  of  radiation  damage 

to  any  material  x where  Dx(E)  is  known,  the  prescriptions  out- 
lined here  address  themselves  specifically  to  radiation  damage 
in  silicon;  the  damage  function  for  silicon,  Dg^(E),  is 
rapidly  varying  near  1 MeV.  This  method  consequently  pre- 
scribes an  alternative  definition  for  0 (silicon) . 

eq 

1.4  This  method  is  part  of  a set  of  standards  dealing 
with  neutron  spectrometry,  and  is  intended,  to  serve  as  a means 
of  characterizing  different  neutron  spectra,  0(E),  used  in 
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radiation  damage  studies  and  parts  testing  of  electronic 
components.  A standard  method  of  deducing  #(E)  by  unfolding 
neutron  threshold-foil  activation  data  is  covered  in  ASTM 
Method  E XX3,  "Standard  Method  for  Unfolding  Neutron  Spectra 
for  Radiation  Hardness  Testing."  The  procedures  for  select- 
ing, irradiating  and  counting  the  foils,  the  tabulation  of 
some  of  *e  required  nuclear  constants,  and  the  method  of 
calculating  specific  activations  of  the  foils  are  presented 
in  ASTM  Methods  E XXI,  "Irradiation  a Standard  Set  of  Neutron 
Threshold  Activation  Foils  for  Radiation  Hardness  Testing," 
and  E XX2,  "Measuring  Foil  Activities  for  Radiation  Hardness 
Testing."  Once  4>e(j  and  4>eq/$  have  been  measured  and  calculated, 
it  is  important  for  the  user  of  a neutron  irradiation  facility 
to  be  able  to  measure  for  subsequent  irradiations.  This 
is  covered  in  ASTM  Method  E XX5 , "Measuring  the  Relative  1-MeV 
Silicon-Equivalent  Fluence  with  Fast-Neutron  Monitors." 

1.5  A flow  chart  shown  in  Figure  1 lists  the  sequence 
of  operations  that  lead  to  $(E)  . 4>eq,  'J>eq/$*  and  4>eq/monitor , 

2.  SIGNIFICANCE 

2.1  In  neutron  radiation  damage  studies,  where  the  damage 
is  neutron-energy  dependent,  it  is  convenient  to  parameterize  a 
neutron  field  in  terms  of  a single  number  called  the  1-M#V 
equivalent  fluence,  t ^eq’  quantifies  it  by  a 

single  number  and  a standard  deviation  of  that  number.  In 
addition,  it  is  useful  to  assign  a radiation-damage  quality 
to  a given  spectral  distribution,  a parameter  referred  to  here 
as  the  hardness  parameter.  This  hardness  parameter  ia  simply 
<J>  /t,  or  the  number  of  1-MeV  neutrons  that  produce  the  same 

radiation  damage  as  a unit  fluence  of  neutrons  having  Che 
spectral  distribution  $(E). 
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3. 


SUMMARY  OF  METHOD 

3.1  The  parameter  4^  is  defined  as 


L 


0.01  MeV 


♦(E)  D(E)dE 


D(1  MeV) 


(1) 


where  ♦(£)  is  the  spectral  distribution  of  the  neutrons  and 
D(E)  the  radiation  damage  per  unit  fluence.  D(1  MeV)  is  taken 
as  the  average  0(0.85-1.13  MeV)  because  D(E)  is  rapidly  varying 
near  1 MeV.  The  integral  does  not  Include  neutrons  below  0.01 
MeV  because  even  for  a soft  TRIGA  spectrum,  less  than  1 per- 
cent of  the  damage  in  bulk  silicon  is  contributed  by  neutrons 
below  this  energy.  The  0.01-MeV  cutoff  greatly  simplifies 
the  computation  of  the  hardness  parameter,  ♦ /♦: 


♦ /♦  ■ 
yeq  r 


0 

h 


01  MeV 


♦(E)  D(E)dE 


D(1  MeV) 


(2) 


0.01  MeV 


♦(E)dE 


The  "neutron  count" , given  by  the  integral  in  the  denominator 
of  Eq.  2,  includes  only  about  half  of  the  epicadmium  neutron 
flux  for  a TRIGA-type  reactor,  and  inclusion  of  the  entire 
epicadmium  flux  would  make  ♦ /♦  very  dependent  on  neutrons 
below  0.01  MeV,  where  almost  no  damage  occurs.  If  the  integral 
extended  down  very  low  energies,  to  include  thermal  neutrons, 
the  situation  would  be  much  worse.  The  accurate  assessment 
of  the  value  of  8eq/$  would  then  require  an  accurate  measure 
of  neutrons  in  the  thermal  region  - an  unnecessary  complication 
and  a great  source  of  meaningless  variation  of  ♦ /♦,  since  a 

small  amount  of  thermal -neutron  shielding  would  change  the  para- 
meter drastically  without  changing  ♦(£)  above  0.01  MeV,  where 
almost  all  the  damage  occurs. 
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4. 


CALCULATIONS 


4.1  The  output  of  the  SAND  II  unfolding  code  is  given 
as  a 620-point  spectrum,  from  10'^  ^ 0 ig  MeV,  or  260  points 
in  the  region  of  0.01  to  18  MeV.  A recent  calculation  of  the 
silicon  displacement  KERMA,  D(E) , has  been  carried  out  with 
the  latest  neutron  cross  sections  for  silicon  (1,2),  and  is 
given  at  200  energy  points  covering  roughly  the  same  neutron- 
energy  interval.  In  folding  #(E)  and  D(E) , a simple  numerical 
integration  is  carried  out  as  follows. 


18  MeV  620 

J ♦(£)  D(E)dE  - ^ *(Et)  D(Ei)AEi  , (3) 

i-361 


where  i - 361  corresponds  to  E^  - 0.01  MeV  of  the  SAND  II 
output,  and  the  D(E^)  are  the  values  given  in  Table  I for 
the  same  energy  mesh  as  the  SAND  II  output  above  0.01  MeV. 
They  rperesent  values  obtained  from  the  original  table  (1,2) 
by  interpolation  and  some  group  averaging,  where  necessary. 

4.2  The  value  given  by  Eq . 3 is  divided  by  D(1  MeV)* 
6(0.85-1.15  MeV  ■ 84  MeV  -mb)  to  provide  $ as  per  Eq . 1. 
#eq/4>  is  calculated  by  simply  carrying  the  summation 


18  MeV 


*1^1 


(4) 


0.01  MeV 


to  obtain 


r 

Jo. 


<p(E)dE 
01  MeV 
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in  the  denominator  of  Eq.  2.  The  error  of  the  simple  sum- 
mation is  less  than  1 percent  for  the  620-point  spectrum 
that  vanishes  at  the  upper  limit  of  integration  (18  MeV) . 

5 . PROBLEMS 

5.1  D(E)  is  strongly  varying  near  1 MeV.  Its  value, 

as  given  in  Ref.  1,  is  102  MeV  mb  at  1 MeV,  while  the  average 
value  from  0.85  to  1.15  MeV  is  84  MeV 'mb.  Although  this  dif- 
ference does  not  impact  inter- laboratory  comparisons  of  silicon- 
damage  effects  if  a standard  value  of  D ( 1 MeV)  is  selected,  it 
leads  to  an  unrealistically  small  value  for  $ /$;  a fluence 

of  0.89  neutrons  of  1 MeV  becomes  equivalent  to  a unit  fluence 
of  GODIVA-type  neutrons  (50  cm  from  reactor) . The  GODIVA  50-cm 
spectrum  is  the  hardest  of  three  spectra  measured  in  an  evalua- 
tion program  (2)  supporting  these  draft  specifications.  If, 
however,  D (0.85  to  1.15  MeV)  is  used  in  Eq.  2 instead  of 
D(1  MeV),  the  1-MeV  equivalent  fluence  per  unit  fluence  now 
becomes  1.06;  a much  more  intuitively  meaningful  figure  for 
neutrons  whose  average  energy  is  greater  than  1.0  MeV,  weighted 
with  a damage  function  Dg^(E)  that  increases  with  E. 

6.  RECOMMENDED  DEFINITION 

6.1  Use  5(0.85  to  1.15  MeV)  - 84  MeV -mb  as  D(1  MeV) 
in  Eqs.  1 and  2. 

6.2  Use  Table  I for  D(E^) , where  the  E^  mesh  agrees 
with  that  of  SAND  II. 

7.  PRECISION 

7.1  The  estimated  precision  in  calculating  $ and 
$eq/$  with  the  present  method  is  appreciably  better  than 
1 percent,  assuming  4>  (E)  and  D(E)  are  exact,  and  that  the 
SAND  II  620-point  spectrum  is  used  in  obtaining  the  fluence 
above  0.01  MeV  (Eq.  2)  (260  points  above  0.01  MeV). 
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8 . ACCURACY 


8.1  The  accuracy  of  the  numerical  integration  is 
estimated  to  be  better  than  1 percent.  The  accuracy  of 
D(E)  can  be  usefully  discussed  in  three  parts. 

8.2  Shape  of  D(E) 

8.2.1  In  the  calculation  of  a and  a /♦,  only  the 
shape  of  the  D(E)  curve  is  important,  as  related  to  inter- 
laboratory comparisons  of  these  quantities,  because  D(E) 
appears  in  both  the  numerator  and  denominator  of  Eqs . 1 and 

2.  The  shape  of  D(E)  is  probably  known  well  enough  to  contri- 
bute no  more  than  3 to  5 percent  to  the  relative  values  of 
and  $eq/$  obtained  for  different  spectra  41(E)  . 

3 . 3 The  Absolute  Value  of  D(E) 

8.3.1  The  latest  D(E)  data  (1,2)  were  obtained  with 
the  most  recent  evaluated  cross-section  set  for  silicon, 
which  provides  the  total  energy  deposited  in  Si  to  an 
accuracy  that  would,  in  itself,  impact  $ and  $eq/$  by 

no  more  than  a few  percent  (say  2 to  3 percent)  because  of 
its  effect  on  the  shape  of  D(E) . However,  the  calculation 
of  the  partition  of  this  energy  between  ionization  and  dis- 
placement is  much  less  well  known.  The  partition  was  cal- 
culated with  the  Lindhard  theory  (3,4),  which  has  not  been 
modified  to  agree  with  experimental  results  (5) . Con- 
sequently, the  D(E)  values  may  change  by  the  order  of 
10  percent  when  appropriate  adjustments  are  made  to  obtain 
good  agreement  with  past  and  future  measurements. 

8.4  D( 1 MeV) 

8.4.1  The  value  of  a „ and  a will  vary  about 

eq  eq 

20  percent,  depending  upon  the  interval  over  which 

5(1  MeV  t dE)  is  averaged.  There  is  a great  deal  of 

arbitrariness  in  the  absolute  value  of  and  * /$  in 

eq  eq 


10b 


chi*  respect.  in  chat  the  definition  does  not  influence 
intercomparisons  of  these  quentities  for  most  reactor-type 
spectre  used  in  radiation  damage  studies.  The  interval  0.85 
to  1.15  MeV  was  chosen  as  a compromise  between  a very  wide 
energy  interval,  over  which  *(E)  may  encounter  significantly 
large  variations,  and  a delta  function  chat  results  in  an 
intuitively  low  value  for  ♦ /♦  because  of  the  rapidly  vary- 

ing (and  therefore  unrepresentative)  value  of  D(E)  at  or  very 
near  1 MeV  (see  Discussion,  Sec.  5). 
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TABLE  I.  D(E)  in  MeVmb  (E  in  MeV) 


E 

0 

E 

0 

E 

0 

E 

0 

1.000-2 

1 . 054 

9.600-2 

S.356 

9.20-1 

96.1 

5.20*0 

193.3 

l.OSO-2 

1 . 086 

1 . 000- 1 

5.298 

• 9.60-1 

104.8 

5.30*0 

168.1 

1.100-2 

1.112 

1.050-1 

5.269 

1.00*0 

72.5 

5.40*0 

182.5 

1.150-2 

1.144 

1.100-1 

5.182 

1.10*0 

62.6 

5.50*0 

176.4 

1.200-2 

1.164 

1.150-1 

4 . 803 

1.20*0 

44.0 

5.60*0 

170.6 

1.275-2 

1.208 

1.200-1 

4.367 

1.30*0 

69.9 

5.70*0 

175.5 

1 . 350-2 

1.252 

1.275-1 

3.734 

1.40*0 

92.6 

5.80*0 

181.9 

1.42S-2 

1.287 

1.350-1 

3.202 

1.50*0 

98.7 

5.90*0 

188.6 

1.500-2 

1.322 

1.425-1 

2.911 

1.60*0 

101.9 

6.00*0 

196.5 

1.600-2 

1.377 

1.500-1 

2.911 

1.70*0 

120. S 

6. 10*0 

189.2 

1.700-2 

1.429 

1.600-1 

7.569 

1.80*0 

96.6 

6.20*0 

177.6 

1.800-2 

1.479 

1.700-1 

16.30 

1.90*0 

166.8 

6.30*0 

166.5 

1.900-2 

1.563 

1.800-1 

66.95 

2.00*0 

177.0 

6.40*0 

155.7 

2.000-2 

1.572 

1.90-1 

101.9 

2.10*0 

184.8 

6.50*0 

147.3 

2.100-2 

1.630 

2.00-1 

10S.7 

2.20*0 

142.6 

6.60*0 

151.4 

2.200-2 

1.688 

2.10-1 

91.1 

2.30*0 

133.0 

6.70*0 

155.7 

2.500-2 

1.747 

2.20-1 

81.5 

2.40*0 

154.3 

6.60*0 

161.0 

2.400-2 

1.805 

2.30-1 

71.3 

2.50*0 

157.2 

6.90*0 

165.5 

2.550-2 

1.913 

2.40-1 

64.6 

2.60*0 

142.6 

7.00*0 

170.9 

2.700-2 

1.965 

2.55- 1 

59.4 

2.70*0 

157.2 

7.10*0 

175.0 

2.800-2 

2.038 

2.70-1 

SS.6 

2.80*0 

197.9 

7.20*0 

177.3 

3.000-2 

2.140 

2.80-1 

54.4 

2.90*0 

185.4 

7.30*0 

179.3 

3.200-2 

2.212 

3.00-1 

53.0 

3.00*0 

142.6 

7.40*0 

181. 1 

3.400-2 

2.387 

3.20-1 

52.1 

3.10*0 

127. S 

7.50*0 

183. 4 

3.600-2 

2.515 

3.40-1 

51.5 

3.20*0 

123.4 

7.60*0 

183. 1 

3.800-2 

2.678 

3.60-1 

51.2 

3.30*0 

118.8 

7.70*0 

181.4 

4.000-2 

2.887 

3. SO-  1 

52.1 

3.40*0 

122.3 

7.80*0 

180.2 

4.250-2 

3.0S7 

4.00-1 

52.7 

3.50*0 

124.6 

7.90*0 

178.7 

4.500-2 

3.348 

4.25-1 

53.6 

3.60*0 

123.4 

8.00*0 

177.3 

4 . 750-2 

3.581 

4.50-1 

54.7 

3.70*0 

122.6 

S. 10*0 

176.1 

S. 000-2 

3.930 

4.75-1 

55.6 

3.80*0 

121.7 

S . 20*0 

173.0 

5.250-2 

4.221 

5.00-1 

59.7 

3.90*0 

121.7 

8.30*0 

174.7 

5.500-2 

4.512 

5.25-1 

80.1 

4.00*0 

123.1 

S.40*0 

174.4 

5.750-2 

S.  094 

5.50-1 

147.0 

4.10*0 

124.3 

8.50*0 

1 73 . S 

6.000-2 

S.  706 

5.75-1 

113. S 

4.20*0 

129.8 

8.60*0 

173.2 

6.300-2 

6.404 

6.00-1 

61.1 

4.30*0 

148.5 

S.70*0 

172.9 

6.600-2 

6.986 

6. 30- 1 

55. 0 

.4.40*0 

1 68 . 8 

s.so*o 

1'3.2 

6.900-2 

7.248 

6.60-1 

56. S 

4 . S0*0 

180.2 

S.  90*0 

2 

7.200-2 

7.132 

6.90-1 

59.4 

4.60*0 

1SS.4 

9.00*0 

1S1.6 

7.600-2 

6.  783 

7.20-1 

64.6 

4.70*0 

lSS.o 

9.10*0 

1S6.6 

8.000-2 

6.462 

7.60-1 

76.6 

4.80*0 

192.7 

9.20*0 

190.4 

8.400-2 

6. 171 

8 . 00- 1 

92.6 

4.90*0 

195.0 

9.30*0 

193.0 

8.800-2 

5.822 

8.40-1 

96.0 

5.00*0 

194 .5 

9.40*0 

3 9S.2 

9.200-2 

5.560 

8.80-1 

S2.7 

S.10»0 

193.9 

9.50*0 

2CC.3 
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E D E 'D 


9. 

60*0 

199 

.4 

9. 

70*0 

196 

.8 

9. 

60*0 

19S 

.0 

9. 

90*0 

197 

.1 

1. 

00*1 

196 

.5 

1. 

01*1 

195 

.6 

1. 

02*1 

19S 

. 0 

1. 

03*1 

195 

. 3 

1. 

04  *1 

195 

.6 

1. 

05*1 

196 

7 

• * 

1. 

06*1 

196 

.5 

1. 

07*1 

196 

.8 

1. 

08*1 

197 

.1 

1. 

09*1 

197 

.4 

1. 

10*1 

193 

.0 

1. 

11*1 

199 

.1 

1. 

12*1 

200 

.6 

1. 

13*1 

202, 

.0 

1. 

14*1 

203 

.5 

1. 

15*1 

204 

.9 

1. 

16*1 

206 

.4 

1. 

17*1 

207 

.8 

1. 

18*1 

209, 

,6 

1. 

19*1 

210, 

.5 

1. 

20*1 

210, 

.6 

1. 

21*1 

211. 

,0 

1. 

22*1 

211. 

,6 

1. 

23*1 

211. 

9 

1. 

24*1 

209. 

6 

1. 

2S*1 

212. 

5 

1. 

26*1 

212. 

8 

1. 

27*1 

213. 

1 

1. 

28*1 

213. 

4 

1. 

29*1 

213. 

f 

1. 

30*1 

213. 

7 

1. 

31*1 

213. 

7 

1. 

32*1 

214. 

0 

1. 

33*1 

214. 

0 

1. 

34*1 

214. 

*> 

to 

1. 

35*1 

214. 

•> 

1. 

36*1 

214. 

5 

1. 

37*1 

214. 

5 

1. 

38*1 

214. 

8 

1. 

39*1 

214. 

8 

1. 

40*1 

214. 

8 

1.41*]  214.8 

1.42*1  214.8 

1.43*1  214.8 

1.44*1  214.8 

1.45*1  215.1 

1.46*1  215.1 

1.47*1  215.1 

1.48*1  215.4 

1.40*1  215.7 

1.50*1  216.3 

1.51*1  216.6 

1.52*1  217. 

1.53*1  217. 

1.54*1  218. 

1.55*1  2i8. 

1.56*1  218. 

1.57*1  219. 

1.58*1  219. 

1.59*1  220. 

1.60*1  220. 

1.61*1  220. 

1.62*1  220. 

1.63*1  220. 

1.64*1  220. 

1.65*1  220. 

1.66*1  220. 

1.67*1  220. 

1.68*1  220. 

1.69*1  221. 

1.70*1  221. 

1.71*1  221. 

1.72*1  221. 

1.73*1  222. 

1.74*1  222. 

1.75*1  223. 

1.76*1  223. 

1.77*1  224.4 

1.78*1  224.7 

1.79*1  225.3 

1.30*1  226.2 
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DESIGNATION: 


E XX5 


Standard  Method  for 

MEASURING  THE  RELATIVE  1-MeV 
SILICON  EQUIVALENT  FLUENCE 
WITH  FAST -NEUTRON  MONITORS 


1 . SCOPE 

1.1  This  method  describes  the  measurement  of 
^eq/M°nit°r > l~MeV  equivalent-neutron  fluence  per  unit 
monitor  count,  for  a neutron  field  where  the  spectrum  $(E) 

has  been  measured  and  the  1-MeV  equivalent  fluence  for  silicon 
radiation  damage  has  been  calculated. 

1.2  This  method  is  part  of,  and  follows  as  a corollary 

to,  the  four  ASTM  methods  that  are  addressed  to  measuring  the 
neutron  spectrum  (E)  , and  with  it,  calculating  These 

are  as  follows: 

E XXI  Irradiating  a Standard  Set  of  Neutron 

Threshold  Activation  Foils  for  Radiation 
Hardness  Testing 

E XX2  Measuring  Foil  Activities  for  Radiation 
Hardness  Testing 

E XX3  Unfolding  Neutron  Spectra  for  Radiation 
Hardness  Testing 

E XX4  Characterizing  Neutron  Spectra  in  Terms 

of  1-MeV  Equivalent  Fluence  for  Radiation 
Damage  in  Silicon 

2 . SIGNIFICANCE 

2.1  This  method  is  user  oriented,  in  that  it  is 
addressed  to  measuring  /Monitor  for  neutron  irradiations 

ecj 
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subsequent  to  the  one  in  which  4> (E ) was  measured.  It 
provides  4>e^/Monitor  for  a variety  of  options,  depending 
on  (a)  the  access  to  the  results  of  the  SAND  II  code  used 
to  deduce  4> (E)  from  the  threshold  data,  and  (b)  the  avail- 
ability of  the  apparatus  and  facilities  required  for  calibrat- 
ing various  types  of  monitor  foils. 

3 . APPARATUS 

3.1  Nickel  or  iron  foils  and  gamma-ray  detector  as  in 
ASTM  Method  E XX2 . 

3.2  Sulfur  foils,  beta  counter,  aluminum  foil,  Ge(Li) 
or  intrinsic  germanium  detector,  NBS  gamma-ray  sources 
(calibration  standards) , and  access  to  14-MeV  neutron 
generator . 

4.  DETERMINING  * PER  UNIT  MONITOR  COUNT 

eq 

4 . 1 By  Measurement 

4.1.1  During  the  course  of  measuring  <KE)  Co  obtain 
4> eq  * a monitor  foil  may  be  placed  in  the  same  neutron  field 
as  the  set  of  threshold  activation  foils,  or  in  some  more 
convenient  location,  thus  obtaining  ^^/Monitor  directly. 
Subsequent  monitor  counts  will  then  predict  the  absolute 
flux  in  terms  of  if  no  shields  or  nearby  moderators  or 
scatterers  are  changed.  Thus,  a "retractable"  GODIVA-type 
reactor,  calibrated  at,  say  2 meters  above  the  floor,  must 
be  subsequently  relocated  in  height  to  within  0.1  meter  if 
Che  monitor  foils  are  located  at  the  same  distance  from  the 
reactor  as  the  irradiation  samples  (^50  cm) , and  within  a 
few  centimeters  if  the  monitors  are  approximately  a factor 
of  two  nearer  (or  farther) . Nearby  scatterers  must  not  be 
moved  so  as  to  change  the  flux  by  more  than  a few  percent. 

An  upper  limit  on  the  inscattering  effect  of  any  materials, 
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such  as  sample  holders  or  nearby  stands,  can  be  estimated 
(a)  by  using  a rough  average  of  the  total  cross  section  in 
the  vicinity  of  1 MeV  for  reactor  neutrons,  (b)  by  assuming 
isotropic  scattaring,  and  (c)  by  taking  into  account  the  flux 
at  the  scatterer  and  the  solid  angle  subtended  at  the  detector. 
For  thick  materials  such  as  concrete,  a fast  neutron  albedo 
of  10  percent  can  be  used. 

4.1.2  Most  shielding  materials  in  front  of  either 
sample  or  monitor  must  be  kept  constant  in  thickness  within 
a.2  mm  for  less  than  a 5 percent  change  in  flux  {2  mm  for 
polyethylene  and  water  (1)}. 

4.1.3  While  2 or  3 grams  of  additional  moderator  along- 
side a fast  neutron  monitor  will  make  little  difference,  a few 
tenths  of  a gram  of  moderator  adjacent  to  the  cadmium-covered 
zero-threshold  detectors  can  be  important  (i.e.,  1/v  detec- 
tors. and  235U  or  t3^Pu  foils). 

4.1.4  The  monitor  is  generally  chosen  with  a threshold 
high  enough  to  make  it  insensitive  to  neutrons  below  0.01  MeV 

whose  contribution  to  the  radiation  damage  is  negligibly  small. 

32  3° 

Sulfur  S(n,p)  kP  is  usually  selected  because  of  the  high 
threshold  value  for  the  reaction,  the  sensitivity,  the  con- 
venient  half-life  of  the  resultant  “P  radioactivity  (beta  rays)  , 
and  the  ease  of  purifying  sulfur.  Detailed  considerations 
for  using  sulfur  are  described  in  ASTM  Method  E 265-70, 

"Measuring  Fast  Neutron  Flux  by  Radioactivation  of  Sulfur." 
However,  the  3^Fe(n,p) 3Sln  or  3®Ni(n,p)3®Co  reactions  are 
preferable  if  a Ge(Li)  or  intrinsic  germanium  detector  is 
available;  the  absolute  activations  can  then  be  measured 
by  simply  cross  calibrating  against  NBS  standard  gamma-ray 
sources,  providing  the  absolute  value  of  $(E)  if  a spectral- 
shape  measurement  is  available. 
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4.1.5  The  $e(j/Monitor  may  be  measured  at  a different 
time,  if  either  a reliable  secondary  monitor  is  used  (such 
as  a fission  chamber  in  a steady-state  reactor,  or  a ther- 
mometer in  a fast-burst  irradiation).  If  not,  one  of  the 
threshold  foils  used  in  measuring  <J>  (E)  can  be  used  for 
calibration  during  exposures  subsequent  to  measuring  $ (E) . 

4 . 2 By  Calculation 

4.2.1  If  a standard  monitor  foil  such  as  sulfur  is 

not  included  in  the  set  of  foils  irradiated  to  measure  <J>(E)  , 

32 

then  the  S activation  can  be  calculated  with  the  SAND  II 

32 

code.  The  code  is  rerun  with  a trial  value  of  R ( S)  added, 

32  m 

and  R ( S)  is  adjusted  until  A (see  Method  E XX3)  becomes 
m 58  58  0 

zero.  If  the  Ni(n,p)  Co  reaction  had  been  used  in  deter- 
32 

mining  $(E),  the  S activation  can  be  estimated  from  the 
58 

resultant  Co  activity.  The  two  reactions  have  nearly  the 

same  threshold,  and  it  was  found  that  R (^S)  s 2.92  Rm(^®Ni) 

for  three  different  reactor  spectra  (fast  burst  reactor  glory 

hole,  50  cm  from  the  FBR,  and  TRIGA  J-Tube)  . Once  Rffl(^  S) 

is  calculated,  the  sulfur-foil  counter  must  be  calibrated  in 

terms  of  absolute  P beta-ray  activity  resulting  from  the 
32  32 

S(n,p)  p reaction.  One  such  technique  utilizes  14-MeV 

neutrons  to  irradiate  Al  and  S,  the  absolute  14-MeV  neutron 

32 

flux  (and  therefore  the  absolute  P activity)  being  given 

by  the  /Al(n,a)z  wa  reaction.  The  cross  sections  for  both 

reactions  are  accurately  known  at  14  MeV.  Alternatively,  a 
32 

known  amount  of  P (NBS  source)  can  be  mixed  in  with  a 

32 

sulfur  pellet.  This  cross  calibration  of  S is , of  course, 
unnecessary  if  the  ^Ni(n,p)^®Co  reaction,  the  ^Fe(n,p)^Sln 
reaction,  or  some  other  well-known  reaction  leading  to  a gamma- 
ray  emitter  is  used  for  monitoring  <f>  . 
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5. 


PRECISION 


5.1  The  precision  for  measuring  $e(^/Honitor  is  the 
combination  of  the  precision  of  measuring  <Peq  and  of  measur- 
ing the  monitor  count.  The  accuracy  of  measuring  is 
treated  in  ASTM  Method  E XX4,  and  will  not  be  covered  here. 

The  precision  of  measuring  the  ^®Ni(n,p)^®Co  and  ^Fe(n,p)^Sln 
activations  is  about  3 percent  if  counting  statistics  and 
calibration  techniques  allow.  The  foil  location,  the  fluence 
of  each  irradiation,  and  the  counting  times  will  determine 

the  counting  statistics,  and  these  factors  are  largely  subject 
to  the  choice  of  the  user.  The  precision  for  beta  counting 
of  sulfur  foils  is  about  the  same  as  for  the  nickel  and  iron 
monitor  foils,  if  the  sulfur  foils  are  thick  enough  to  make 
the  count-rate  independent  of  thickness.  Many  of  the  sulfur 
counting  considerations  are  discussed  in  ASTM  Method  E 265-70, 
"Measuring  Fast-Neutron  Flux  by  Radioactivation  of  Sulfur." 

6 . ACCURACY 

6.1  If  ^®Ni  and  ^Fe  foils  are  used,  an  accuracy  of 
v5  percent  can  be  achieved  on  the  absolute  activation  if 
counting  statistics  are  not  a consideration  (see  Sec.  5). 

For  sulfur  monitoring,  the  accuracy  is  v5  percent  if  the 
sulfur  is  exposed  simultaneously  with  the  threshold  activa- 
tion foils.  If  not,  and  the  sulfur  activation  must  be 
calibrated  indirectly,  the  accuracy  is  in  this  case  esti- 
mated to  be  10  percent. 
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